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SECTION 1 
SUmARY 

The Lunar Ejecta and Meteorites Experiment (LE4M) was deployed on the 
moon on 12 December 1972. The objectives of the experiment were to me^sure 
the long-term variations in cosmic dust influx rates and the extent and nature 
of the lunar ejecta. While analyzing these characteristics in the data, it 
was discovered that a majority of the events could not be associated with hyper- 
velocity particle impacts of the type usually identified with cosmic dust, but 
could only be correlated with the lunar surface and local sun angle. 

The possibility that charged particles could be incident on the sensors 
led the Principal Investigator (PI) to request that an analysis of the elec- 
tronics be per*formed to determine if such signals could cause the large pulse 
height analysis (PHA) signals. These signals indicate the energy of the hyper- 
velocity particles in the normal mode of operation. 

A qualitative analysis of the PHA circuit showed that an alternative mode 
of operation existed if the input signal were composed of pulses with pulse 
durations very long compared to the durations for which it was designed, by a 
factor of at least 40 to 1. This alternative mode would give large PHA outputs 
even though the actual input amplitudes were small. This revelation led to the 
examination of the sensor and its response to charged particles to determine 
the type of signals that could be expected. 
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A qualitative review of the sensor and application of basic electro- 
static theory indicated that very slow particles, below the normal experiment 
operating range, could produce pulses of the time duration required to excite 
the PHA circuit's anomalous response. 

A grossly simplified model of the sensor was developed on a computer to 
determine the range of particle characteristics to which the sensor would 
respond. This range was then compared with known or expected values for lunar 
dust particles and practical expectations for charge to mass ratios. 

At the same time, the electronics was analyzed using a standard IBM 
analysis program, SCEPTRE. 

The results of the sensor modeling and circuit analysis showed con- 
clusively that charged particles moving at velocities below 1 kilometer per 
second would produce PHA responses of the type observed in the lunai data and 
in addition could cause double accunulator counts, another of the unusual 
events . 

This finding was of such importance to the understanding of lunar sur- 
face dust transport that it was decided to continue the analysis to obtain 
more accurate data on particle mass, charge, and velocity. A theoretical 
calibration of the experiment response to charged particles was required to 
enable a complete analysis of the lunar data to be performed. In addition, 
a practical measurement of the response using the experiment qualification 
model was to be attempted to corroborate the analysis. A complete physical 
calibration was impractical. 
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The analysis was continued on two fronts. A simplified model of the 
electronics was developed because the SCEPTRE simulation was cumbersome and 
costly to use. In parallel with this, a refined model of the sensor was 
developed to remove the limitations of the simple model and provide greater 
accuracy. 

The sensor film, collector grid, and suppressor grid were divided into 
7,360 elements for computational purposes. Using basic electrostatic prin- 
ciples, the charge distributions on each plane were calculated for both the 
applied potentials and the charged particle. The 7,360 simultaneous equations 
that result from the mutual interactions between elements were solved iter- 
atively. The program used a large area of computer memory and was slow to 
converge to a result. No complete results were obtained from this model 
because efforts were made to speed up the convergence and overall running 
time to save future costs. 

Two other programs, which apply the sensor model results to the elec- 
tronics and then analyze the results, were prepared and checked on simulated 
data. Program descriptions are given in the Appendix. 

The conclusion from the analysis to date is that the LEWI experiment 
data contain significant information relative to mechanisms operating at the 
lunar surface. To fully understand and appreciate these mechanisms, the 
lunar events recorded by LEAM must be transposed into parameters of particle 
mass, velocity, and charge and their respective variations in space and time. 
To accomplish this, a calibration of the LEAM in response to charged particles 
must be completed. 
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This report reconmends that the analysis be continued, in conjunction 
with work being performed by the Principal Investigator, to provide a com- 
prehensive picture of the dust environment at the lunar surface. The results 
would be, in addition to characterization of the particles, that unique events 
would be characterized, allowing se^nentation of the measurement range, and 
event types would be correlated with lunar cycles and temporal effects. Hypo- 
theses on dust formation and transport would be refined and opportunities 
would be developed for understanding several unexplained phenomena observed 
on the lunar surface by astronauts and other experimenters. 

A meeting was conducted on 20 July 1976 by the LEAM Principal Investi- 
gator with Dr. W. Quaide and M.J. Smith of NASA Headquarters to discuss the 
present LEAM program status and the importance of continuing both the analysis 
of the experiment response to charged particles and the lunar data analysis. 

A summary of the LEAM study status and the proposed tasks for extended 
study of the charged particle phenomenon is included herein as Appendix B. 
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SECTION 2 
INTRODUCTION 

The study of the LEAN experiment's response to charged particles was 
initiated at the request of the Principal Investigator, when it was observed 
that data over a 2-year period showed an incidence of signals with outputs 
of 6 and 7 PHA counts, far greater than anticipated from data obtained on 
previous space flights. Particles of this energy would normally penetrate 
the front film and provide signals at the rear film, but this was not ob- 
served. There were numerous events which recorded impacts on two film strips 
or collector grid strips, or which recorded two accumulator counts for one 
event. These events could not be explained by the normal experiment response 
to hypervelocity particles. The average event rate of less than 10 particles 
per 3-hour period gave an extremely low probability of two particles being 
incident on the sensor at precisely the same time. The inhibit circuit, which 
was employed to prevent crosstalk between adjacent sensor elements, prevents 
noncoincident events from being recorded in the same time frame. This guar- 
antees that PHA and accumulator data can be identified with the correct event. 
The majority of the events occurred around sunrise and sunset, but thermally 
induced signals were ruled out because the onset of the data occurred up to 
60 hours before sunrise, when the experiment was thermally stable. Normal 
ope'^ation of the experiment was verified by the internal calibration signals, 
which were generated automatically every 15.5 hours. 
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The preliminary analysis was discussed in detail in a Bendix report, 
ASTIR/TM66, prepared 1 August 1975. The electronics analysis using SCEPTRE 
showed that for long input pulses to the PHA peak detector the diode in the 
forward path continued to conduct and maintain the input to the threshold 
detector. This, in turn, allowed the PHA counter to continue incrementing. 

In addition, if the pulse length and amplitude were above certain levels, a 
condition arose which caused double counting of the film accumulator. The 
accumulator increments whenever the PHA threshold detector is triggered. 

Double triggering was caused by the combination of pulse length, amplitude, 
and the circuit time constants. The circuit was designed for pulses of 2 
microseconds maximum length, while the pulses giving the effects discussed 
above were over 80 microseconds in length. 

To determine the type of signal to be expected from the sensor in res- 
ponse to charged particles, a very simple model of the sensor was developed 
which treated the sensor planes as solid conducting sheets rather than 95% 
transparent grids. The model permitted an increased understanding of the 
electrostatic principles involved and allowed determination, within an order 
of magnitude, of the ranges of particle parameters to which the sensor would 
respoi ’. 

The simple sensor results showed that the electrostatic forces involved 
were significant for particles of masses and charges in a range which could 
reasonably be expected to be present on the moon. Also, if the velocities 
were below 1 kilometer per second (km/sec), signal pulse lengths and amplitudes 
could be obtained from the film which would cause the PHA circuit to give the 
observed large values and double accumulator counting. 
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Thus, the simple sensor and SCEPTRE analysis showed that LEAN could 
respond to slowly moving charged particles and give data outputs similar to 
those observed on the moon. The simple model could not give accurate values 
for the mass and charge ranges measurable by the experiment because of Its 
gross simplification of the electric fields. Also, It did not Include any 
modeling of the film strips adjacent to the one being considered, which meant 
that multiple events and Inhibits were Ignored and PHA signal levels were 
generally too small. 

To alleviate the limitations of the simple sensor and to provide an 
electronic model which would provide cost-effective results, a refined sen- 
sor model and a simple electronics model were developed. The refined sensor 
model included a true representation of the grid structures and the inter- 
actions between elements. 
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SECTION 3 

METHOD OF ANALYSIS 

3.1 REVIEW OF LEAN OPERATION 

3.1.1 Sensor Operation 

The sensor (Figure 3-1) normally operates upon impact of a particle 
that causes ionization of film material at the impact site. This ionization 
is collected at the film and collector grid. The negative potential of the 
film attracts the positive ions while the positive potential of the collector 
grid attracts the electrons. These actions cause anal! current flows in the 
film and collector grid circuits, which result in a positive voltage pulse to 
the film amplifier and a negative voltage pulse to the collector grid ampli- 
fier. The film and collector grid areas are divided into 1-inch strips, which 
allow for identification of the impact site. 

A second film and grid assembly is situated behind the first and sep- 
arated from it by 5 centimeters. The operation of this rear assembly is 
similar to that of the front assembly. An analysis of impact locations on 
the two films provides an indication of the direction of travel of the particle, 
while the time taken to traverse the intervening front and rear film space 
provides a measure of particle velocity. 

3.1.2 Electronics Operation 

The typical dual sensor logic is divided into two sections, the first 
ra^k or measurement section, and the rear rank, or buffer section. The meas- 
urement section includes identification pulse storage latches, accumulators. 
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Figure 3-1 The Basic Sensor 
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PHA conversion counters, and TOP conversion counter. The rear rank is a 
parallel-in, serial-out shift register, which shifts data, upon demand, to 
the ALSEP central station in predetermined telemetry frames. The shift 
register is only cleared when new data are to be transferred into it and 
a<"ter the old data have been transmitted to ALSEP at least once. The new 
data are transferred from the front rank storage latches to the shift regis- 
ter, provided that the current frame is not one in which daca are to be 
transferred to ALSEP. If the old data have been transferred to ALSEP once, 
the new data are retained in the front rank storage latches, thus allowing 
data from two events to be retained. Further hits in rapid succession would 
be evidenced by accumulator counts only. The time interval during which 
rapidly occurring events, which exceed the storage capability, would be lost 
varies between 2 milliseconds and 3 seconds, depending upon the position of 
the telemetry sequence in ALSEP. Data have not been observed which approach 
this event frequency. 

The pertinent circuits for this analysis are those associated with 
the front film as shown in Figure 3-2, which shows the elements of one typical 
film channel. These are the circuits which were previously referred to as 
the front rank or measurement section. There are two distinct signal channels 
beyond the film amplifier; (1) the film strip identification channel, or film ID, 
and (2) the pulse height analysis, or PHA channel. Each is discussed separ- 
ately. 
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Figure 3-2 Typical Film Channel 





BSR 4234 


3. 1.2.1 Film Strip ID Channel 

The connion film amplifier provides the -3-volt film bias and a non- 
inverting gain of 3. The output is applied to the first amplifier of the ID 

channel, the PHA amplifier, and the analog inhibit inputs of the three other 
film channels. 

The ID amplifier provides an inverting gain of 5.25 at its nonnal 
input and a gain of 0.49 at each of three noninverting inputs, which receive 
analog inhibit signals from the other film amplifiers. These inhibits cause 
the output of the ID amplifier to renain at or above 0 Vdc if one of more of 
the other films receive a coincident signal which is approximately 10 times 

greater than that on film 1. If film 1 has a signal equal to or greater than 

the other films, an output is applied to the threshold detector. The thres- 
hold detector is designed to apply a logic "1“ to the following NWID gate if 
the input signal at the film amplifier exceeds 1 millivolt (mV). The NAND 
gate sets the following latch circuit, provided that the ID inhibit signal 
from the central electronics is also at logic "1", indicating that no other 
front film latch is set. The latch circuit provides the signal to the output, 
via a buffer, to indicate which film strip has been impacted. 

When the ID latch is set, all the ID signals are inhibited for the 
four front film strips, which has the effect of negating crosstalk and makes 
the ID channels for the front film unresponsive until the measurement cycle 
for this hit is completed. 
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The OR fimction of all the film and collector latches and the micro- 
phone sample one-shot signal starts a measurement cycle. If a collector latch 
only is set during the l-mi11isecond (msec) measurement period, a normal se- 
quence occurs, except that the data transfers and clear are inhibited while 
the clear latch signal is generated. Thus, a collector signal alone will not 
be presented in the data output nor will existing data be changed. 

When the system start occurs, a l-mil1isecond gate signal is gen- 
erated which has three functions: 

1. Provide an enable to the front and rear PHA counters. 

2. Provide a synthetic rear film signa. jo complete the time of 

flight sequence if the normal signal does not occur within 1 
millisecond. 

3. Prevent a premature measurement completion signal . 

The film count accumulator measures PHA signal threshold crossings, 
providing that a film ID latch is set. The ID latches are inhibited for any 
further hits during a measurement cycle, but the accumulator circuit may give 
evidence of later hits. If a second hit occurs within the PHA pulse of the 
first, the PHA is augmented, but no direct evidence of the second hit survives. 

If the second hit is delayed sufficiently to create an independent PHA pulse, 

but lies still within the 1-msec measurement gate, it will cause further PHA 
counting and one additional increment to the film accumulator. If it occurs 
more than 1 msec after the first hit, but before transfer of data into the 
shift register, it will cause an increment of the accumulator only. 
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.2.2 Film Signal Pulse Height Analysis 

The signals from the four film amplifiers are sunned by the PHA 
a>nplifier which, together with Wie film amplifier, gives a gain of -10 from 
film strip to PHA amplifier output. This output is passed to liie Peak Detec- 
circuit, which is a high-gain mnplifier with a closed-loop gain of -*-1.0 
^ negative signals. The detector charges the capacitor C to the peak of the 
irput signal. When the input signal is removed, the diode in the forward path 
prevents discharge of the capacitor C back through the amplifier. 

When transistor 7X is on. capacitor C discharges with a time con- 
stant that is designed to give a 240-niicrosec decay time. The voltage across 
the capacitor is sensed by the PHA threshold detector, which is a high-gain 
operational amplifier. When the voltage across capacitor C is more negative 
than -10 mV, the detector output is clamped at -0.6 V, the "0" level for the 
logic inver-.r of the following stage. When the voltage is more positive than 
-10 mV, a logic "1" (+2.5 V) is presented to the inverter input. 

When the voltage on capacitor C is below threshold, the logic gates 
hold transistor TX on, which causes capacitor C to be in a short time con- 
stant mode, "’en threshold is achieved, transistor TX is turned off via the 
logic u’"’i. the next 25-kHz clock pulse sets the flip-flop. When the flip- 
floo IS set, transistor TX turns on (allowing capacitor C to discharge), the 
rhA counter is enabled, and the accunulator is incremented. When the capa- 
citor discharoes to below threshold level, the threshold detector causes the 
flip- floe to be reset and the PHA counter to be disabled. The length of the 
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pulse from the flip-flop, and thus the length of time the PHA counter is en- 
abled, is proportional to the peak of the input pulse. Thus, the count re- 
corded by the PHA counter is a measure of the pulse height. 

The synchronization of the capacitor discharge with the 25-kHz 
clock reduces the quantizing error. 

3.2 ANALYSIS OF PULSE HEIGHT ANALYSIS (PHA) CIRCUIT 

The description of operation given in Section 3.1 applies to the type of 
particle for which the experiment was designed. That is, a noncharged, hyper- 
velocity particle which would cause a pulse input to the electronics with the 
fol lowing characteri sti cs : 


Amplitude 

1 to 200 mV peak 

Rise Time 

400 nanoseconds (nsec) 

Fall Time 

1,000 nsec 

Width 

600 nsec 


The experiment was tested and qualified for this type of input under all con- 
ditions of lunar environment, and thus shown to meet the design i equirements. 

When considering the effects of charged particles upon the sensor, it 
was realized that, for slow particles, current pulses of much greater length 
than 2 microsec could be obtained. (The sensor dynamics are discussed in 
later sections.) The PHA circuit was then analyzed for the effects of long 
input pulses. 
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A qualitative review of the peak detector circuit shows that, for a 
short pulse, the capacitor C is charged to tto peak of the input signal and 
the decay time of the charge on the capacitor is proportional to this peak 
value. The time constant in this mode is approximately 45 microsec, which 
was chosen to give the maximum count of 7 in 240 microsec. (The PHA output 
indicates at least 1 whenever a threshold is achieved.) When a long pulse 
occurs, the diode in the forward path is held in a conducting state, even 
while capacitor C is being discharged in what is normally called the short 
time constant mode. The effect of the conducting diode is that the signal 
is maintained at the amplifier output. The result at capacitor C is to 
effectively increase the time constant by 200 times, thereby maintaining the 
signal above threshold for a much longer time. The longest pulse which will 
not change the PHA value is theoretically 80 microsec, but the value depends 
upon the time relationship between the start of the pulse and the 25-kHz clock 
and could be less than 80 microsec. 

In addition to the extended count for long pulses, a condition arises 
that causes double accumulator counts. If a pulse of sufficient amplitude 
and length occurs, the falling edge of the pulse causes the input to the peak 
detector to go hard positive, shutting off the diode. The capacitor C now 
discharges normally. The time constants ahead of the peak detector are such 
that its input returns to a negative value, which causes the diode to conduct 
again. If the capacitor C had previously discharged below threshold and the 
signal is large enough (negative) to exceed threshold again, an extra accum- 
ulator count is made and renewed PHA counting occurs. 
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The above analysis Indicates that negative pulses can also give PHA 
thresholds. 

The qualitative analysis was followed by a detailed quantitative anal- 
ysis of the electronics and by laboratory tests on the experiment prototype. 

3.2.1 Circuit Analysis 

The circuit analysis was performed on the typical film channel of 
Figure 3-2 (from the film input to the input of the PHA threshold detector). 
The emphasis was placed on the peak detector portion of the film cNnnel since 
this is the circuit which gives rise to extended counting and multiple accum- 
ulator counts. The remainder of the circuitry was simulated by passive net- 
works and fixed gain terms. 

A detailed simulation was performed using the SCEPTRE* computer pro- 
gram to give a thorough understanding of the circuit operation under all con- 
ditions. This knowledge was then used to develop a simple model of the cir- 
cuits because the SCEPTRE program used an excessive amount of computer time 
for this component configuration. This long run time would make the task very 
expensive for the multiple computations we planned over the ranges of mass, 
charge, and velocity applicable to the problem. 

3.2. 1.1 SCEPTRE Simulation 

The simulation program, SCEPTRE, was developed by IBM for the Air 
Force Weapons Laboratory at Kirtland Air Force Base, New Mexico. The pro- 
gram calculates initial conditions, and transient and steady- state responses 
for large networks. 

*Bowers, J.C. and Sedore, S.R., "SCEPTRE: A Computer Program for Circuit and 
System Analysis," Prentice-Hall, Inc. 1971. 
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The flln and PHA amplifiers were simulated by a simple gain term 
and the transistor TX was assumed to be In the fully conducting state* I.e.* 
ON; thus, the fl1p>fIop and logic control of transistor TX were not simulated. 
The linear transistors were simulated In the nonlinear regions with the best 
data available. The peak detector circuit Is shown In Figure 3-3. 

A typical output from a run is shown In Figure 3-4. The output Is 
the voltage across capacitor C, shown as positive because of the sign conven- 
tion used In the simulation. The output is observed to return negative at 
700 microsec, but on this occasion the amplitude was Insufficient to cause 
further PHA or accumulator counting. 

A summary of the data obtained from several simulations Is shown 
In Table 3-1. All runs were made foi 1-msec duration, which is the measure- 
ment sample time. The times quoted are the length of time the output pulse 
remained above 9 mV, which 1^ the threshold level at the following detector 
circuit. The data show that PHA levels of 7 can be achieved with Inputs of 
30 mV and the multiple pulses do occur. 

The simulation program provides Information on all the Intermediate 
points within the circuit. This Information was used to Identify critical 
components and, thus, enable us to devise a simple model of the circuits. 

Computations were made on Identical data Inputs, using both SCEPTRE 
and the simple model to verify the latter's validity. 
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Figure 3-4 Typical SCEPTRE Output 
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Table 3-1 

SCEPTRE Program Results 


Input 
Anpl Itude 

Output Pulse 
Length (microsec) 

Coments 

50-microsec Pulse 



50 mV 

184 

No subsequent pulses - All normal 

100 

211 


150 

234 


100-microsec Pulse 



50 mV 

213 

No subsequent pulse 

100 

243 

Returned above 9 mV at 730 microsec 
until 890 microsec 

150 

260 

Returned above 9 mV at 670 microsec 
until 1.01 msec 

200-microsec Pulse 



50 mV 

248 

Returned above 9 mV at 720 microsec 

300-microsec Pulse 



10 mV 

189.9 


20 

219.6 


30 

235.0 

Returned above 9 mV at « 780 
microsec 

40 

245.6 

Returned above 9 mV at ^760 
microsec 

50 

254.03 

Returned above 9 mV at a»746 


microsec 

A11 longer than normal 
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3. 2. 1.2 Simplified Peak Detector Model 

Analyzing the data from the SCEPTRE program Identified the Impor- 
tance of the various components within the peak detector, thus allowing us 
to eliminate many of them without affecting the veracity of the result. 

The obvious simplifications are to neglect the transistor ‘'‘'^emal 
capacitances as they are small and the ^tssoclated time constants h ? 
effect on the result. Next, the coupling capacitors 1n the forward one feed- 
back paths are found to have no effect on the length of time the output re- 
mains above threshold or on the cause of the double accumulator counts. 

When the Input signal 1s negative, the circuit behaves as a simple 
amplifier with a gain of 1. When the signal Is positive-going, the diode 
ceases to conduct, allowing capacitor C to discharge. Once the diode ceases 
to conduct, the feedback loop opens and a large back bias is applied due to 
the high open loop gain. The diode will not conduct again until a forward 
bias is applied from the combined effects of the capacitor discharge and 
input level. In the simplified model. Figure 3-5, the diode is replaced by 
a switch, which opens whenever the input increases positively faster than the 
rate at which the voltage across R12 increases. The rate of rise of the 
voltage across R12 is calculated for the switch-open conditions. (The switch 
closes when a forward bias is achieved.) 

The input to the peak detector is an emitter follower with a paral- 
lel capacitor across its load. The effect of the capacitor is to restrict 
the rate at which the emitter can rise towards tho +5-volt supply line. Con- 
sequently, the input transistor cuts off if this input signal rises positively 
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faster than the eaitter load can follow. The ealtter follower just described 
is replaced by a switch whose condition depends upon the direction and rate of 
change of the Input signal. 

The loading of the ewitter follower upon the coupling circuit be- 
tween the PHA anpllfler and the peak detector Is snail, so the coupling cir- 
cuit Is t ' . ated as an Indei^ndent elenent. The signal source VS2 for the 
peak detector Is then the output of the coupling circuit. Siwllarly, the film 
amplifier loading of the coupling circuit between Itself and the film Is small, 
allowing these components to be treated Independently. The signal source V$1 
Is -10 times the voltage across resistor R32 because the film and PHA ampli- 
fiers, together, give an Inverting gain of 10. 

The simple model , Figure 3-5, Is thus comprised of a unity gain 
amplifier, two voitage sources, two switches, and 12 passive components. The 
mcdel has four possible operating conditicms: 

1. Switches A and B closed. 

2. Switch A open, switch B closed. 

3. Switch A closed, switch B open. 

4. Switches A and B open. 

The input signal from the film is divided into many elemental ramp 
functions with known initial value, slope, and time duration. The response 

of the model to such a ramp is calculated (for all four conditions) using 

Laplace transform techniques. The correct response to be applied for any 

particular rmnp element is deternined by first deducing the state of switches 

A and B at the end of the time interval. For »»xamDle, if the switches are 
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Initially closed and a particular ramp input would cause switch B to be open 
at the end of this time interval , the true signal values at the various points 
in the model are calculated using the condition 3 equations. The time incre- 
ments are chosm to be small enough that the errors incurred due to opening 
switch B slightly early are negligible. 

A further complication of the model is that, for large signals, one 
or all of the film, PHA, or peak detector amplifiers can saturate. This cw- 
dition is accounted for using the ramp technique, where the relevant amplifier 
output is treated as a ramp with zero slope. 

3. 2. 1.3 Complete Film Channel Model 

The remainder of the film channel of Figure 3-2 was modeled to 
simulate the correct lEM response to the sensor signals. 

The film and collector grid 10 model accounts for the analog inhibit 
signals from the three sensor elements, at either the film or collector grid, 
respectively, which are not impacted by the particle. A charged particle, 
unlike an uncharged meteorite particle, can induce signals in adjacent sensor 
elements. This affects the charge/velocity characteristics of the particle 
required to achieve threshold, because the inhibit signal from one element 
effectively reduces the signal from an adjacent element. In addition, the 
timing of the element IDs relative to one another and between films and 
collector grids is modeled. The inhibit signals prevent multiple film IDs 
unless they occur within approximately 0.2 microsec of one another. This 
limitation also applies to the collector grids. When a film or collector 
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grid 10 is received, the system starts a measurement sequence with the setting 
of a bracket one-shot udiich lasts for 1 msec. If a collector signal starts 
the sequence, a film ID must be received within this 1-msec period or no data 
transfer takes niace. A film ID alone can cause the system to operate through 
its full measurement sec^ience. 

Vihen a film ID is indicated, the four film signals are summed and 
applied to the peak detector model. Tl« output is recorded for PHA count and 
accumulator count. The accumulator counts PHA threshold crcisings. The PHA 
count is limited to 7 in the LEAK, but in the model it is allowed to reach its 
full value of 26 if a long enough pulse occurs. This is done to obtain more 
information about the response. 

3.2.2 Laboratory Tests 

Measurements were made using the Prototype LEAM Experiment, the experi- 
ment test set, a variable pulse width generator, and a storage oscilloscope. 

The LEAM center support structure was removed from the outer housing and ther- 
mal bag, and the east sensor was removed fron the center support structure. 

This dismantling was required to allow access to the microphone board upon 
which the PHA circuitry resides. The sensor circuitry was now without shield- 
ing, which meant that it was very susceptible to noise, making other than 
qualitative measurements difficult. 

Pulse inputs were injected via the test set calibration adapter box, 
with the input pulse amplitude being measured directly on the film input test 
point. 
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Measurements were made on the A film channels 1 and 2, which gave 
Identical results as follows: 


Input 


Output 


Pulse 

Width 

(micro- 

second 

2 


100 

300 


2 

6 

50 

200 


Pulse 

Amplitude 

(millivolts) 

4.5 

6.5 
28 


30 

30 


-100 
- 28 
- 5 
-1.5 


PHA of 1 registered on test set Imps. 

PHA of 2 registered on test set imps. 

At capacitor C: -250 mV peak pulse; rise 
time 1 microsec; fall time 
to -10 mV, 120 microsec. 

At flip-flop output: 4.5-volt logic pulse 

120-microsec width. 

First noticeable change at flip-flop output. 

Output at flip-flop; logic pulse greater than 
200-microsec width, starting at threshold 
crossing. Second pulse at 950 microsec from 
threshold, greater than 20-microsec width. 
Occassional multiple pulses occurred around 
950 microsec from threshold. 

PHA threshold. 

PHA threshold. 

PHA threshold. 

PHA threshold. 


In summary, the laboratory tests showed that long pulses give large 
PHA counts with the actual value depending upon pulse amplitude and duration. 
Multiple pulses can occur, which add to the PHA count if they occur during 
the 1-msec sample period, and increment the film hit accumulator, giving the 
appearance of multiple film hits. These tests also confirmed that negative 
pulses at the film input can give PHA and accumulator outputs. 
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3.3 REFINED SENSOR MODEL 

A previous report, ASTIR/1N66, detailed the analysis which led to a 
simple model of the sensor. This s1i^)1e model verified that the sensor can 
give valid responses to charged particles with certain mass, charge, and 
velocity characteristics. The model has several limitations which made It 
difficult or, in some cases, impossible to accurately pre«i'ct the response 
to certain particle types, and also gave undetermined inaccuracies in the 
results. 

3.3.1 Simple Model and Its Limitations 

The simple model was based on an analysis that considered the grids 
and film to be infinite plane conducting sheets. This was modified at the 
grids by applying a simple cosine function to the forces on the particle to 
allow the force to go to zero in the grid planes. 

The limitations of the simple model were: 

1. Solid electrodes were used instead of grids with 95% transpar- 
ency. Thus, the grid signals and forces due to induced charges 
were overestimated. 

2. There was no interaction accounted for between the suppressor/ 
collector space and the film/collector space. Thus, the film 
could not see the particle until it passed the collector grid. 

3. Induced charges were calculated by assuming the 1-inch by 4- 
inch strips were circles of equivalent area. 
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4. Only one film strip and collector grid strip <«ere considered, 
whereas a particle will Induce charges In all film strips and 
collector grid strips. This prevents considerations of multiple 
element events at the film or collector grids and gives Inaccurate 
values for particle characteristics which can cause PHA thres- 
holds. 

5. The analysis only considered particle positions between the sup- 
pressor grid and film, with no account being taken of the forces 
on the particle outside the sensor. Thus, all calculations 
assume a particle emerging from the suppressor grid, on the 
film side, with a certain velocity. The true sensor measurement 
range Is not calculated, as the suppressor, due to Its potential, 
will accelerate positive particles and decelerate negative parti- 
cles, while the image forces accelerate all particles. 

To overcome the limitations of the simple model and thus obtain a more 
complete and accurate result, a different approach was utilized to refine the 
model. 

3.3.2 Refined Model 

The sensor is composed of three parallel planes, termed the film, col- 
lector grid, and suppressor grid. The film and collector grid planes are each 
divided into four 1-in. by 4-in. strips and each strip is composed of four 1-in. 
by 1-in. squares. Thus, each plane has 16 1-in. by 1-in. segments. The sup- 
pressor grid is formed by one plane divided into a similar set of 16 segments. 
One of the 1-in. by 1-in. square sections is shown in Figure 3-6. 
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Figure 3-6 LEAM Grid and Film Structure 
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The problem resolves Itself Into two areas, namely the charges Induced 
in the saisor and the potential at the particle. The change In the Induced 
charge as the particle position changes gives a measure of the current Into the 
sensor electronics, while the difference In potential betirfeen successive parti- 
cle positions gives a measure of the work done by the particle and, hence, 
enables calculation of the velocity profile along the path. 

The charges on the sensor elements arise from two sources, the charges 
due to the applied potentials and the charges due to the particle. Both distri- 
butions are require wo determine the potential at the particle, while only the 
latter is required to determine the current flow due to particle movements. 

The potential at the particle is thus seen to be from two sources, the applied 
potential charge and Its own Induced charge. This latter effect Is similar to 
the Image effects used on the simple model. 

The task of modeling the sensor was complicated by several factors. 

The major problem was containing the model within a size that could be handled 
by the computer. The job is equivalent to solving nearly 8,000 simultaneous 
equations. It rapidly became obvious that a compromise had to be reached 
between accuracy and the nunrt)er of elements into which the sensor films and 
grids could be divided. A secondary problem associated with the number of 
elements is that of devising a satisfactory bookkeeping scheme for keeping 
track of which element is influencing which. This task also is affected 
strongly by programming limitations of array dimension sizes and allowable 
DO loop nesting. The final model has 7,360 elements which between them have 
over 27 million interactions. Considerable effort was expended in accommo- 
dating these interactions within 132,701 influence coefficients. The use of 
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this reduced nunber of coefficients required careful bookkeeping and the 
formulation of generalized equations that ^(pressed the relationships of the 
elements to tiie coefficients. 

The coefficients could not all be retained in memory at the same :1me, 
so they were calculated and retained on magnetic tape and called upon when 
required. The most efficient method for operating the sensor model would be 
to have all the coefficients available at once, but as this was not possible, 
a compromise of using two sets of coefficients at a time was used to speed up 
the Iterative process. The two largest coefficients take up 130,000 bytes of 
core. 

The sensor physical shape precludes Its being easily divided into uni- 
formly sized elements. Allied to this Is the task of calculating the Inter- 
actions between the various elements. As the configurations and shapes are 
not found In standard text books, all the Interactions for the potentials pro- 
duced at one element by a charge on another were calculated from elementary 
electrostatic principles. 

The film and grids are divided into 7,360 uniformly charged elements, 
which are 0.125 in. on a side. The charge distributions due to the particle 
and the applied potentials are calculated separately and superposed. 

In either case, the charge on an element is adjusted so that the total 
potential, caused by its own charge and that due to all other element charges 
and the particle if considered, is equal to to the applied potential. The 
charge adjustment is made iteratively by changing the charge on each element 
to the newly determined value at each iteration. The applied potentials are 
set to zero for calculations of the charge due to the particle. 
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The Iterations are continued until the changes In charge distribution 
at each step are less than a specified value, I.e., the calculation has con- 
varged to within an acceptable tolerance of the final value. 

All calculations and results are In terms of a unit coulomb charge 
on the particle. The potential of each element due to all other elements of 
the sensor Is calculated using a set of stored "Influence coefficients." 

These coefficients are the values of potential at an element due to a unit 
charge at another element. To save computer time, they were calculated once 
using first principles of electrostatics and stored for future use. A similar 
set of coefficients is calculated for each particle position, but they are 
determined in real time for each new particle path. 

A computer program was prepared to perform these calculations. Several 
options are made available which are selected by input variables or cards. 

The basic calculations are: (1) to calculate the charge distributions due to 
the applied potentials and store them on tape; these distributions are fixed 
and used often; (2) to calculate the charge distribution due to the particle; 
and (3) to calculate the potential at the particle due to (a) the applied 
potential charge distribution and (b) the particle image charge distribution. 
Items 2 and 3 are repeated for each position of the particle. The charges on 
each film strip and collector grid strip are summed to give the total charge 
on each element at each step. The data relative to a particular particle 
path are stored on tape for future use. 
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3.4 SYSTEM MODEL 

To determine the response of the LEAM to a charged particle, the data 
obtained from the sensor model are used as an Input to the electronics model . 
The sensor model output Is the characteristics of a particular path through 
the sensor calculated using a particle of unit charge. The system mode' 
uses these data In conjunction with the parameters for the particular parti- 
cle In question to derive the actual response to that particle. Thus, the 
profile of the current flow In each film and collector grid strip is deter- 
mined versus time. The profile is then applied to the electronic model as 
discrete ramp inputs for each time interval. 

A program was prepared to accomplish this which performs the following 

tasks: 

1. Reads input cards to determine i^ich of the following options to 
perform: 

a. Selection of sensor, up, east or west and particle path. 

b. Normal or shielded film on east sensor. 

c. Positively or negatively charged particles. 

d. Preselected or random mass and charge values. 

e. Number of particles. 

f. Particle velocities. 

g. Whether output is to be plotted and, if so, the dimensions 
of the axes. 

h. How many of the data points to list on output. 
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2. If a plot Is desired, the plot progron data are generated. 

3. If random particle characteristics are desired, a random nunber 
generator Is employed to derive mass and charge values. 

4. Data relevant to selected particle path read from tape. 

5. Calculates work done on particle between successive steps and 
calculates velocity at each step. 

6. Calculates currents In films and collector grids from rate of 
change of charge. 

7. Determines if film and collector grid IDs occur. 

8. Calculates Input signal to PHA circuit. 

9. If a film ID occurs, the electronics model subroutine Is called 

to calculate PHA and accumulator response to the signal calculated 
In step 8. 

10. Results are listed or plotted as selected by Input cards. All 
results are stored by sensor on tape for future analysis. 

Thus, a single particle path can be analyzed for either positively or 
negatively charged particles at any number of velocities, charges, and masses. 
The stored data for any sensor and any path can then be analyzed by a second 
program, which Is designed to select the particles by type of event or velocity 
and can either plot or list the resulting selection. The types of events that 
can be selected, either singly or In combination, are coincidence, noncoinci- 
dence, multiple accumulator, multiple film or collector grid adjacent or non- 
ad jacent, on any of the sensors or shielded film. 
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The orientations of the film and collector grid strips within the 
LEAM experiment are Identified In Figure 3-7. This Information Is supplied 
so that the analysis data can be readily compared with the lunar data. 
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SECTION 4 

RESULTS AND IMPLICATIONS OF ANALYSIS 

4.1 RESULTS 

An accurate simplified representation of the electronics has been ach- 
ieved in a computer model. This model simulates the inhibit circuits in addi- 
tion to the PHA threshold circuit analyzed previously. 

Vfhen the simplified electronics model Mas completed, it was checked 
out with the simple sensor model. This combined model gave useful results 
because it could be used with the random number generator to generate numer- 
ous particles with differing mass and charge values and calculate the result- 
ing responses very quickly compared with the SCEPTRE program. 

The plots resulting from these runs are shown in Figures 4-1, 4-2, and 
4-3. The PHA values and ('‘Mjble accumulator events appear in bands which 
differ in shape, depending upon the velocity of the particles. The separa- 
tion of events into those with and without double accumulator counts will 
permit a broad classification of the particles observed on the moon. 

The intent with the rt'ined sensor model was that at least one particle 
path would be calculated and analyzed by the end of the contract period end- 
ing on 31 July 1976. 

we have achieved the following towards this goal. A program to calcu- 
late the influence coefficients for the interactions between the 8512 sensor 
elements was prepired, debugged, and 132,701 coefficients committed to 
magnetic tape storage. The sensor program that utilizes these coefficients 
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has been written, debugged, and operated. The main part of this program is 
the iterative loop, which adjusts the element charges to the values needed to 

give the required potentials both in the case of the applied potential distri- 
bution and the distribution due to a particle. Several problems were encoun- 
tered in the implementation of this iterative loop: 

1. The most efficient method of implementation involves holding the 
132,701 coefficients in core while performing the iterations, but 
this takes 530,804 bytes of memory, which is virtually the entire 
capability of the computer. Thus, a method was devised which re- 
quired repeatedly reading the coefficients from tape in blocks. 

2. The calculation of the potential contributions at each element 
due to all the other elements is the most time-consuming portion 
of the iterative loop. The initial implementation of this part 
took almost 30 minutes per iteration to run. Considerable effort 
was expended in reducing the running time until we achieved the 
present t me of approximately 17 minutes, which was done by stream- 
lining each of the 15 subsections of this part and then combining 
them where possible. The number of elements was reduced from 8,512 
to the present number of 7,360 by considering the tops and under- 
sides of the grids as single elements. This potentially impairs 
accuracy, but the difference is insignificant in our model. Finally, 
the whole part was formulated as a subroutine and compiled using 

the FORTRAN H compiler. 
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3. The present probleo is ensuring rapid convergence of the iterative 
loop. Khen originally formulated, the loop was conditionally stable, 
depending upon the magnitude of the changes made in the elemental 
charges at each step. When stable, the convergence was extremely 
slow because of the small size of the changes in charge which were 
permissible. Although time consuming, the present program will 
provide the required data. 

The remaining tasks to achieve the one particle path for one sensor, 
once convergence is achieved, are; 

1. To perform one run of the program to determine the charge distri- 
bution due to the applied potentials. 

2. To perform 10 runs of the program to determine the distributions 
due to the particle. It is assumed that 10 data points will be 
sufficient to allow a good interpolation for the intermediate data 
points. 

3. To perform interpolation to obtain all other required data. 

4. To run sensor and electronics model program. 

4.2 IMPLICATIONS OF ANALYSIS 

The analysis as performed to date indicates that nearly all types of 
events observed on LE/W can be explained and that classification by event type 
will allow more accurate identification of particle mass, charge, and velocity 
characteristics. 
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The hypotheses explaining the events are described below. Mhen the 
model Is made fully operational, the hypotiieses will be verified. 

The coincident film and collector grid events were shown by the simple 
model to be obtained by a positive particle, between the collector grid and 
film, traveling toward the film. 

Noncoincident events can be achieved by a positive particle with a com- 
bination of mass, charge, and velocity that provides sufficient signal at the 
film but not at the collector grid. The collector grid is less sensitive to 
charged particles. Noncoincidence at the collector grid cannot be observed 
because the experiment requires a film ID to allow completion of a measure- 
ment sequence. 

Multiple accumulator events have been observed with the simple model 
and are caused by the electronics response to long duration input signals. 

Multiple adjacent film events are caused by a positive particle having 
a combination of iTwss, charge, and velocity that give a sufficiently large 
signal to achieve threshold on two or more films at once. The same mechanism 
would be expected to result in multiple collector grid events, but conceivably 
it could give only a single one if the signal level were in the right range. 

Multiple nonadjacent film events are of the type where films 1 and 3 
recorded an ID threshold but film 2 did not. This phenomenon can be ex- 
plained by a negatively charged particle traveling toward the film strip 
that does not record an ID threshold, e.g., film 2. It will be remembered 
that the film circuit requires a positive current to produce an ID, which, 
in the case of a positive particle, was achieved by an induced negative 
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charge In the film. This charge was produced by a flow of electrons 
to the film, equivalent to a positive conventional current flow into 
the amplifier. In the case of a negative particle, a positive induced 
charge occurs in the film and, thus, a negative current flows to the ampli- 
fier. This current will not produce an ID, as observed by film 2. Consider 
now films 1 and 3. If the particle has appropriate charge and velocity 
characteristics, it will induce sizable positive charges and, thus, negative 
current flows in them also. As the particle approaches the plane of the film, 
its influence on films 1 and 3 will decrease, falling eventually to zero at 
the film. Note that this is not the case with film 2 whose charge increases 
until impact. Thus, the charges at films 1 and 3 reach a peak positive value 
somewhere before the film and then decrease to zero at impact. When the 
charge starts to fall to zero, there is an electron flow to the film to re- 
place the positive charge; this flow is again the positive conventional cur- 
rent flow into the amplifier. Therefore, if the magnitudes are correct, suf- 
ficient current can flow to produce an ID in films 1 and 3. 

Shielded film events are explained by the fact that the thin dielectric 
virtually has no effect on the particle induced charge in the film except to 
restrict the approach of the particle to it. Thus, the induced signals will 
be identical to the unshielded films for particles in similar positions. 

The following observed cases in the lunar data are less easy to explain 
and require assumptions which cannot yet be proven: 
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1. Multiple film, nonadjacent, events with no collector ID. 

2. Multiple collector, adjacent and nonadjacent, with single film 
ID. 

3. Multiple film and multiple collector, both nonadjacent. 

Analyzing these cases requires further knowledge of the effects of the particle 
on the film when It is In electrode spaces other than the collector grid/film 
space. If the particle can truly induce a signal of threshold amplitude in the 
film when It 1s in these areas, then the remaining cases can probably be ex- 
plained. 

The detailed study of the sensor and electronics has led to a better 
overall understanding of the instrunent responses and has indicated areas that 
affect the LEAM data but which must be left to future analysis. 

Our analysis considers only particles traveling perpendicularly to the 
film. Obviously, particles are likely to be traveling in all directions. 
Particles traveling at the speeds considered here would probably describe 
curved paths in the proximity of the sensor elements, and this has not been 
considered. The implication is that particles, outside the field of view for 
hypervelocity particles, could be electrostatically deflected into the inst-u- 
ment if they have appropriate energy and charge characteristics. 

The verification that the LEAM experiment is measuring charged dust 
particles as well as hypervelocity cosmic dust particles could lead to an 
understanding of phenomena observed by astronauts and other experimenters. 
Observations in this category include several instances of solar light scatter- 
ing over the terminator regions reported by the Apollo crews in lunar orbit. 
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transient lunar events being investigated by experimenters on a worldwide 
basis, and indications at the Apollo 17 site that a substantial amount of 
lunar surface material has been added over the past 1 to 2 million years.* 


Abstracts of Papers Submitted to the Seventh Lunar Science Conference, 
March 15-19, 1976. 


4-10 



BSR 4234 


SECTION 5 

CONCLUSIONS AND RECGmENDATIONS 

There are several conclusions Mhich can be drawn from instrument 
analysis alone, without reference to the lunar data. 

The sensor definitely responds to charged particles that have certain 
ranges of mass, charge, and velocity. The physical dimensions and applied 
potentials of the sensor are such that charged particles incident upon it are 
affected dynamically and some particle selection takes place. Charged parti- 
cles can be attracted into the sensor, thereby increasing its effective field 
of view. In theory, negative particles will cause sensor responses. 

The electronics does not differentiate between signals from hypervelocity 
particles and charged particles, but the circuits are sensitive to pulse shape. 
The pulses from hypervelocity particles, for which the experiment was designed, 
are well defined, both from theory and gun measurements. They are known to be 
of short duration, whereas the sensor analysis has shown that long pulses, 
several hundred microseconds in length, can be produced. The electronics 
analysis has shown that several characteristic responses to long pulses can 
explain certain peculiarities in the LEAM lunar data, namely large PHA counts 
and double accumulator counts. Negative pulses will also give PHA thresholds. 

When comparisons are made between the analyses and the lunar data, it 
can be concluded that different particle types are producing the observed 
events. Some of the events are probably due to particles within a small por- 
tion of the total response range, while some are certainly produced by negative 
particles. 
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The overall conclusion Is that the combined theoretical analysis of the 
electronics and sensor together with the Principal Investigator's analysis of 
LEAH lunar data can provide a comprehensive picture of the dust environment 
at the lunar surface. Therefore, It Is recomnended that the sensor analysis 
be completed in order to allow a thorough analysis and understanding of the 
LEAM lunar data. The achievements to be expected from further study are: 

1. Total ranges of mass, charge, and velocity of particles b hg 
measured by the LEAM Instrimient. 

2. Characterization of particles producing unique events, thus sub- 
dividing total measurement range Into Identifiable segments. 

3. Correlation of particle types Identified In 1 and 2 with lunar 
cycles and temporal effects. 

4. Knowledge gained above will allow refinement of hypotheses on dust 
sources and transport. 

5. Application of icsults to analysis of other lunar surface phen- 
omena observed by astronauts and other experimenters. 

6. Application of results to Pioneer experiment data, allowing addi- 
tional Information to be obtained on deep space particles. 

In accordance with NASA policy, the LEAM experiment data and supporting 
documentation will be archived to make It available for future use by Investi- 
gators anywhere In the world. This report and the results of the Qualification 
model tests constitute essential supporting documentation Invaluable to future 
users of the LEAM experiment data. The bulk of the experiment data Is Incom- 
prehensible without a detailed knowledge of Its response to charged particles. 
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Thus, without this knowledge, the data cannot be applied to investigutions 
of other lunar surface phenomena. Future users of the data could apply the 
results herein to a continued analysis resulting in a comprehensive calibra- 
tion of the instrimient, which would include particles incident anywhere on 
all three sensors. 

A more practical and cost-effective approach would be to require the 
Principal Investigator and the Bendix Project Engineer for the LE/W experi- 
ment to continue the analysis using the extensive knowledge and understanding 
which they have acquired over the past three years. The result would be a 
set of data and documentation with far greater application to other areas of 
scientific research into lunar phenomena than is presently practicable. 
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APPENDIX A 

COMPUTER PROGR-4M DESCRIPTIONS 

The computer programs required for a complete theoretical analysis of 
■3 LE/tfi experiment are described in the following sections. Flow charts 
and listings are included for information purposes. 

The programs complement each other to achieve the fina'. results. The 
numbers given are from the program numbering system for computer data sets, 
used fy the Bendix Corporation Data Center. 

Program P5072CH6 coinputes the path data using subroutines PLFIN. and 
POTCON. Th<» outputs, which are stored on tape, are utilized by P5072o6F to 
determine the experiment response to particular particles. The subrout ne 
used is LES, which itself uses subroutines CONDI, C0ND2, C0N03, and CVOLT. 
Finally, the PHA and accumulator count data for the various particles are 
analyzed or sorted by P5072INT. 

All programs were written ir FORTRAN IV for the IBM-370 sys<-em. The 
plotting routines are those used by the Cal Comp plotting system. 

A.l PROGRAM P5072CHG TO DETERMINE SENSOR CHARACTERISTICS TO CHARGED 
PARTICLES 

A. 1.1 Summary 

The program calculates 

1. Charge distribution on the film, collector grid, and suppressor 
grid due to (a) applied potentials and (b) charged particle. 
These distributions are calculated separately and the one for 
applied potentials is committed to tape fur future use. Those 
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due to the particle are calculated for particle positions* 

Mhich are selected by Input card. 

2. Total charges on each film and grid strip for each particle posi- 
tion. Thus* knowing the particle speed* the current In the film 
and collector grid circuits may be determined. (Ihls calculation 
Is performed In program P5072SGF). 

3. Potential at the particle due to both the applied potentials and 
the particle Image charge. This allows calculation of the work 
done on the particle along the path. 

The program stores position* potentials* and charges on tape so that all param- 
eters for one path are stored for future use. 

A. 1.2 Description 

The calculations center upon determining the charge distributions on 
the films* collector grids, and suppressor grid. The distributions on one 
grid are affected by the distributions on all other films and grids and vice 
versa. Thus, to determine the actual distribution is an iterative process 
which adjusts the individual charge distributions until the calculated poten- 
tial at any element, grid or film matches the applied potentials. When the 
charge distribution due to the particle Is determined, the applied potentials 
are set to zero. 

The films and grids are divided into uniformly charged square elements 
of 3.175 X 10"^ meter on a side. The total number of elements used is 7,360. 
The interactions between elements are determined in a subroutine POTCON using 
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influence coefficients, which have been previously calculated and stored on 
tape. An influence coefficient is the value of potential at one element due 
to a unit charge at another element. 

The resulting charge distribution is used in two ways. The first swns 
the elemental charges on each film and collector grid to give the total charge 
on the respective sensor element at that time. This is done in the particle 
case only and gives the charge due to the particle at each chosen position 
relative to the sensor. The rate of change of charge, caused by particle move- 
ment, determines the sensor output current. The second use for the charge 
distributions is to calculate the potential at the particle caused by both the 
applied potential charge distribution and the distribution due to the particle 
itself. The latter gives rise to the method of images used for calculations 
involving infinite planes. The change in potential along the path through the 
sensor determines the work done on the particle and thus the change .n its 
energy. 

The program has two basic modes of operation: 

1. To calculate the charge distribution due to the applied potentials 
and commit the values to tape. 

2. To calculate the required parameters of potential at the particle 
and total charge on each film and collector grid strip, for each 
selected particle position. 

Other operational options, which are variations and combinations of 
the above two modes, are available and will be discussed later. 
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A.1.2.1 Mode 1 

The mode is selected by guide parameter G1 = 1 on the second input 
card, and guide parameter G2 is set to zero. The initial elemental charges 
are set to half the values estimated for uniformly charged surfaces at the 
potentials of the film, collector grid, and suppressor, and the elemental 
potentials are set to zero. Next, the elemental potentials due to all other 
charges are calculated using the initial charge values and the influence co- 
efficients, which are read from tape. The difference between the potential at 
an element and the applied potential is due to the element's own charge and 
form factor. The charge, thus calculated, is compared with the original charge 
to determine the charge value for the next iteration. 

The comparison includes a check to ensure that the calculated value 
does not lie outside the limits prescribed on an input card. If it is outside 
the limits, the elemental values are scaled to give the limit value for the total 
charge. The charge value for the next iteration is determined by taking a frac- 
tion of the difference between the calculated and original values and adding 
it to the original value. The fraction is selected on the input card, together 
with the number of iterations allowed and the maximum percentage difference 
desired between successive charge values on any element. The maximum percen- 
tage difference detennines the accuracy of the resulting dist. ibution. 

When the program transfers out of the loop, the calculated charge 
distributions are recorded on tape for future use. The transfe** occurs when 
either the iterations allowed are completed or the desired accuracy is 
achieved. 
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A.1.2.2 Mode 2 

This mode Is selected by guide parameters G1 and G2 being set to 
3.0 and 1.0, respectively. 

In this mode, the first step is to calculate the influence co- 
efficients between the particle and the elements of the films and grids and 
vice versa. These coefficients, designated P--Q, are the values of potential 
at an element for a unit charge at the particle and vice versa. The coef- 
ficients are calculated for every particle position that is selected by an 
input card. Subroutine PLEINF is used in the calculation. The charges on 
the films and grids and the potential at the particle are calculated as 
follows, 

1. The charge distribution due to the particle is calculated 
iteratively in an identical manner to that for the applied 
potentials, except that the applied potentials are set to 
zero and the initial element potentials are set to the values 
attributable to the particle (the values of the influence 
coefficients, P— Q). The potential contributions at each 
element due to all other elements are accumulated with the 
P— Q value to give the total potential at each elemeiit. This 
value is compared with the applied potential (now zero) as 
before, and the new elemental charge is determined using the 
same factor. The same criteria are applied as in Mode 1 to 
determine when sufficient iterations have been performed. 
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2. The potential at the par '^'le due to the applied potentials Is 
computed from the Influence coefficients (P— Q) and the charge 
distribution stored on tape in Mode 1. 

3. The potential at the particle due to the charge it induces in 
the films and grids is computed from the influence coefficients 
(P — Q) and the charge distribution calculated for the particle 
alone. 

4. The total charge on each film and collector grid strip is cal- 
culated by summing the respective elemental charges for each 
strip. 

When all the potentials and charges have been computed for a parti- 
cular position, the values are committed to tape as part of a data set which 
is compiled for each path through the sensor. 

The program then reads the next input card for a new particle posi- 
tion. At each position, the program automatically alternates between the 
loop that reads the applied potential charges from tape and the loop that 
iterates to a new charge distribution due to the particle charge. 

A. 1.2. 3 Other Options 

Options are selected by input parameters G1 and G2: 

1. When G1 = 2.0, the program calculates the potential at points 
selected by input cards, in addition to computing and commit- 
ting to tape the charge values of Mode 1. 

2. When G1 = 3.0, the program calculates the potentials of the 
previous option using the charge values recorded on the tape. 
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3. When G1 = 5.0, the charge values recorded on tape in Node 1 
are read in and used as the initial values for the first step 
of the iteration loop. This allows further refinement of the 
charge values without repeating the previous steps. 

4. When G2 = 1.0 and G1 = 0.0, the potentials at the particle 
due to the particle induced charges and the total film and 
collector grid strip charges due to the particle are calculated. 
The potential due to the applied potentials is not calculated. 
This :node has limited use on its own and, if called for, should 
have a dummy card for the JCL card defining FT25F001 to prevent 
erroneous data being stored on a data tape. 

A. 1.2. 4 P5072SIC Program to Calculate Influence Coefficients 

The program to calculate the influence coefficients P5072SIC is used 
once, and the results are stored on magnetic tape. This program calculates 
the coefficients from first principles, based on the physical geometry of the 
elements. The interactions occur many times due to the repetitive nature of 
the physical geometry, but any particular interaction is calculated only once. 
Each interaction is referenced by an index number so that the correct coef- 
ficient can be recalled from tape in program P5072CHG. This program, P5072 lTC, 
determines the correct index number for the particular coordinates of the 
elements under consideration, then calculates the coefficient using subroutine 
INFLCF. 
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All coefficients are stored on tape VOL SER NOS T53344 using the 

following data set names: 

ASD.P067. CFWH 
I CFMSFW 
CFTUFW 
CFEDGU 
CFPFIW 
CFFMFM 
CFTUFM 
CFEDFM 
CFISFM 
CFTUTU 
CFISIS 
CFEDED 
CFEOIS 
I CFEDTU 
ASD.P067. CFTUIS 

A JCL card is required for each data set. 

A. 1.3 Method of Use 

Four input cards are required if full use of the program is to be 
made, including calculations involving particle position. This applies to 
every condition of G1 and G2 except G1 = 1.0 and G2 = 0.0. In this instance, 
the fourth card may be omitted. 

Card 1 controls the iterative process of determining the charge distri 


butions. 


The inputs required, all format code F7.4, are: 
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Columns 1~7; FACTOR* which determines the fraction of old and new charge 
values which are to be used for the value in the next iteration. 

Columns 8-14; PERCEN, specifies the maxinHm percentage difference between 
new and old charge values required before exiting the iteration 
loop. 

Columns 9-21; CYCLES, specifies the maximum number of iterative cycles 
allowed before exiting the loop. 

Card 2 defines the guide numbers 61 and 62 (Format, 2F3.1). 

61 = 0.0 Does nothing with regard to applied potentials. 

61 = 1.0 Charges due to applied potentials are computed and written to tape. 

61 = 2.0 Smne as 61 = 1.0 and also computes the potential at specified point(s) 
from card 4. 

61 = 3.0 Reads charge distribution due to applied potentials from tape and 
computes the potential at specified point(s) from card 4. 

61 = 5.0 Refines charges due to applied potentials. (From Tape). 

62 = 0.0 Does nothing with regard to particle. 

62 = 1.0 Computes charge distribution due to particle. Computes image 
potential at position of particle and total charges on grid and 
film strips due to particle. 

Note: If G1 = 2.0 or 3.0 and/or 62 = 1.0, cards giving XP, YP and ZP must be 
present, where XP, YP and ZP are the coordinates of the particle rela- 
tive to the center of tne film. 
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Card 3 defines the maximum and minimum charge values for each sensor plane 
during the iteration process. These values limit the excursions of 
the charge values to prevent divergence. (Format 6E11.4). 

Card 4 defines the particle path position, the distance of the particle from 
the film and the total number of points (NPTS) to be calculated (par- 
ticle positions). ZP is the distance of the particle from the film 
in meters. XP and YP are the distances from the center of the film 
plane, in meters, as shown below. 


ACl 

AFl 

I AC2 
1 
1 

' ACS 
1 

1 

I 

' AC4 
1 

1 

AF2 

1 

1 

1 

1 

1 

1 

1 XP 

1 

1 

1 

1 

AF3 

1 

1 

1 

1 

1 YP 
1 
1 

1 

1 

1 

AF4 

1 

1 

1 

_J 

1 

1 

1 

J 

1 

1 

1 

1 


A card of this type is required for every particle position or position 
for wi'ich potential due to applied potentials is required. (Format 
3E11.4, 13). XP and YP must have the same respective values on each 
card for each path, i.e., on a particular path only ZP changes. 

Tapes are required for storage of the charges due to applied potentials 
and for the path data which includes potentials and total charges. 
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If the charge distribution due to applied charges is held on tape 
and further refinement of the values is desired, i.e., a smaller value of 
PERCEN, then G1 should be given the value of 5.0. The existing values will 
be read from tape and further iterations performed until the new accuracy is 
achieved. 

Some WRITE statements, that are not shown on the flow chart which 
follows, are included for diagnostic purposes. These print out some of the 
terminal point numbers so that the position in the program can be determined 
and also the potential and charge of selected elements in each plane are 
printed prior to executing terminal points 3508 or 3509. 

A.1.4 Flow Charts and Program Listings 

A flow chart of the program is given in Figure A-1. 

Program listings for P5072CHG, P5072SIC, and subroutines POTCON, PLEINF, 
and INFLCF follow on pages A-13 through A- 53. 
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Figure A-1 P5072CIIG Program Calculates Charge Distributions in Films and Grids 
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PS072CH6 


CARD<S> QIVIMG VftLUES OF 01 « 02 rWST W.M^*BE>RE^NT* 

/ 

Q1»0. e DOES NOTHING WITH RECORD TO RPR.IED POTENTIRLS 
01-1. 0 CHORGES Dl« TO ftPPLIED POTENTI«-S COMPUTED ftt«> WRITTEN TO TAPE 
01-2. 0 SAME AS 01*1. 0 « ALSO COMPUTES POTENTIAL AT SPECIFIED POINT<S> 

01- 3. 0 READS CHARGE DISTRIBUTION DUE TO APPLIED POTENTIALS FROM TAPE A»® 
COMPUTES POTENTIAL AT SPECIFIED POINT<S> 

01=5. 0 REFirCS CHARGES WE TO APPLIED POTENTIW.S. <FROM TAPE > 

02- 0. 0 DOES NOTHING WITH REGARD TO PARTICLE 

02*1. 0 COMPUTES CHARGE DISTRIBUTION DUE TO P«?TICLE. . COMPUTES IMAGE 
POTENTIAL AT POSITION OF PARTICLE. CCW1PUTES TOTAL CHARGES ON 
GRID AND FILM STRIPS DUE TO PARTICLE. 

NOTE IF 01-2. O OR 3. 0 AND/OR 02-1. 0 CARDS GIVING XP, VP ^40 ZP MtST BE 
PRESE T 

NOTE CARDS GIVING VALIKS OF FACTOR, PERCEN AND CVCLES MUST ftt-WAVS BE 
PRE^NT 

IN TKE FOLLOWING ARRAVS THE PREFIX 0 INDICATES T»E TOTAL CHARGE ON fW 
ELEMENT, THE PREFIX P THE TOT«- POTENTIAL AT AN ELEMENT DUE TO ALL OTHER 
CKWRQES ««) THE PREFIX P WITH SirFIX Q THE POTENTIAL AT W ELEI«NT DUE TO 
THE CHARGE ON THE P#«?TICLE ALONE. THE P — Q NUMBERS ARE ALSO THE INFLUENCE 
COEFFICIENTS FOR Tl« EFFECT OF THE ELEMENT CHf^GES UPON THE POTENTIAL 
RT THE PRRT I CLF 

DIMENSION PWQ<£ 2/ 4> A, 7, PTUQ<3. 2, 4/ 4f 2> 8)/ PEDQ<2> 2/ 4f 4, 2, 8> 

DIMENSION PISQ<3:> 4. 4> 2, 2>, PFMSQ<4/ 4^ 8# 8> 

DIMENSION CHG<3;>, SCBLE<3>* QMftX<3>. QMIN<S> 

DIMENSION BLRT^<8> 

COt^iON GFMS<4, 4/ 8# 8 >f PFMS<4. 4, 8> 

COMMON PN<2> 2, 4, 4. 7/ 8>f QIK2, 2> 4, 4/ 7^ 8> 

COMMOf^l QED<2 f 2/ 4/ 4, 2/ 8>f PED<2/ 2* 4, 4, 2> 8> 

COMMOr^l GTU<3, 2, 4. 4, 2, PTUC5> 2, 4. 4, 2/ 8> 

COMMON GIS<3f 4> 4/ 2. 2). PIS<3. 4, 4, 2. 2> 

DRTR BLRNK/8*0. 8/ 

THE RBOVE 15 RRRftVS REQUIRE 25536 WORDS <1. E. 102144 BVTES> 

3500 TO 3514 PROVIDE ROUTING THROUGH THE PROGRRH BLOCKS 
RERD<5. 75:500>FRCTOR. PERCEN. CVCLES 

3500 FORMAT <3F7. 4> 

WRITE <6. 8000 > FACTOR. PERCEN. CVCLES 

8000 FORMAT < 5X. FACTOR = F7. 4. PERCEN = ^ F7. 4. CVCLES « 'fF7. 4> 

PERCENT©. 01>^PFRCEW 

READ<5. 3501>G1. 02 

3501 FORMAT <2F3. 1> 

WRITE<6. 8001>G1. G2 

8001 FORMAT < 5X. ^G1.G2 ^2<F4. 2.2X>> 

RFAD<5. 8020>QMAX. QMJN 

8020 FORMAT <6F11, 4> 

WRITE<6. 8021>QMflX/ QMIN 

8021 FORMAT <5X. ^QMAX M?<Eli. 4/ 5X>/5X/ -'QMIN ^ 3<EJ.l. 4. 5X>> 

IF<G2. GT. 0. 5>LINE«:1 

IFCGl. GT. 0. 5>LINE-0 

JUMP* 0 

INUM«0 

1F<Q1. G1. 4>GD TO 3504 ‘ 

JF<Cil LT. 1. 5. AND. G2. LT. 0. 5>Q0 TO 3504 vlUGlNAL PAGE JB 

OF POOR QUALOX 
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3502 SEM><5. 3583. END»3514»«»/ VP# 2P> riPTS 

3503 F0RHRT<3E11. 4. 13> 

MRITE<6# 8882>XP. VP. ZP 

8802 FORMftKSX. ^XP.VP.ZP « '# 3<E11. 4# !»<>> 

IF<INUn. QT. 1>Q0 TO 5989 
MRITE<25>XP. VP. I8>TS. BLftNK 

GO TO 5999 

C COHPUTES POTENTim.S/irrL. COEFFS DUE TO PPRTICLE RNDRETURNS TO 358^ 

3584 MR1TE<€. 8883> 

8883 F0Rt«T<5X. ^3584' > 

IF<Oi. LT. 2. 5. OR. LINE. EQ. 1>00 TO 2819 
IF<JUTff>. EQ. 1>Q0 TO 2858 
RERD<9>QM. QTU. QED. QIS. QFHS 
REHI^E> 9 

IF<Q1. QT. 4>G^ TO 3585 
60 TO 2858 

C GOES TO 2019 ZOROS RLL CWRGES 6 ^ RETU^4S TO 3585 OR QC^S TO 2858 4 

C CO^FUTES POTENTIAL AT PARTICLE POSITION RETURNING TO 3511 

3585 C0UNT«8. 8 * 

MRITE<E. 8004> 

^!J84 F0RMAT<5X. ^3585-^> 

3588 IF<LINE. EQ. 1>G0 TO 3586 
Vl«-7. 0 
V2«24. 0 
V3«-3. 0 
00 TO 1999 

3506 V1«0. 0 
V2«0. 0 
V3*0. 0 

QO TO 2010 

C GOES TO 1999 A SETS P— TO ZERO. OR GOES TO 2010 ^ SETS P— * P— -Q IN BOTH 
C CASES RETURNING TO 3507 

3507 MRITE<6> S005> 

0005 F0RHAT<5X. ^35070 

IF<COUNT: LT. 0. 5. AND. LINE. EQ. 1>Q0 TO 4999 
GO TO 4499 

C GOES TO ITERATION BLOCK BUT BVPASSES COMPUTATION OF POTENTIAL 

C CONTRIBUTIONS DUE TO ELEMENTS ON FIRST PASS. RETURNS TO 3588 FOR FURTHER 

C ITERATION OR TO 3509 WHEN ITERATION COMPLETED. ALL CHARGES NOW COMPUTED 

3509 WRITE<6. e006> 

8006 F0RMAT<5X. •'35090 

IF<Gi. QT. 4>G0 TO 8850 
IF<Gi. GT. 2 5. OR. LINE. EQ. 1>G0 TO 3510 
8050 WRITE<9>GW. GTU. GED. GIS, GFMS 
REWIND 9 

1F<G1. QT. 4>CiO TO 3514 

IF<G1. LT. :l. 5. AND. Q2. LT. 0. 5>Q0 TO 3514 

JUMP«1 

3510 QO TO 2050 

C COMPUTES POTENTIAL AT PARTICLE. RETURNING TO 3511 

3511 IF<LINE. EQ. 0>G0 TO 3515 
PAR-^SUM 

WRJTE<6. 3517>PAR 
3517 FORMAT <5X. TAR « Ell. 4> 

GO TO 3516 
3515 APP»SUM 

WRITE<6. 3518>APP 
35JA FORMAT < 5X. ’TiPP 4> 

35J6 IF<l.INF. EQ. 0>GO TO *5i3 
GO TO 2108 

C COMPUTES CHARGFS ON FU.M AND GklO STRIP AND RETUPriS TO 3512 
WRITE<6. 3519>AF1. AF2. AES. HF4.. H' 1. AC2> AC3. AC4 
35.1 F0RMHT<5X, •TlL.i CHAkGL »- •*. 4<2X/ Eli. 4>/A1X. GRID CHARGE • ^ 4<2X. E 
2J:I. 4)> 
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WRITE<25>2P. ftPP/ pm, PCI. fK:2. RC3. PC4, ftFl. PF2, BF4 
3513 IF<LI|s^ EO. a RND. G2. GT. 0. 5>L«1 
IF<LINfc. EQ. 1. ftND. Gl. QT. 1. 5>L**0 
IFCLINE. EG. i. ftND. 01. LT. 1. 5>u«l 
lP<Lim. ECl. 0. RND. 02. LT. 0. 5>L»0 
L1NE«L 

IF<G1. GT. 1. 5>G1»3, 8 

IF<L1I^. EQ. 1. ftND. Ql. GT. 2. 5)00 TO 3584 
60 TO 3502 

THIS BLOCK SETS ftLL POTENTIfW. CONTRIBUTIONS FROM ELEMENTS TO ZERO 


1S99 


2000 


2001 


2002 

2003 

2004 

2005 

2006 
200 ? 


DO 2007 K«l>4 
DO 2006 L=l/4 
DO 2005 M«l, 8 
DO 2004 N»l>8 
PFMS<K> L> M. N>^0. 0 
IF<M. GT. 7)G0 TO 2884 
DO 2003 I*=l/3 

IF<M. GT. 2. OR. N. GT. 2)GC TO 200O 

PIS<I.K. L. M^N>a-0. 0 

DO 2002 2 

IF<M. GT. 2)G0 TO 2001 

PTU<I> J. K>L/ M.N>*0. 0 

1F<I. GT. 2)60 TO 2082 

PW<1/ J. K. L,M>N>s^0. 0 

IF<M. GT. 2)60 TO 2002 

PFD<I. J> K. L.M.N)=0. 0 

CONTINUE 

CONTINUE 

CONTINUE 

CONTINUE 

CONTINUE 

CONTINUE 

GO TO j^507 * 

THIS BLOCK SETS RLL ELEMENT CHflRGE£ TO ZERO 




2019 JUr\P=-0 

JFCLINF. EO. 0)60 TO 9008 
/r<TNUM. GT 1>60 TO r<.505 
FVl M^“«0, 0 
GKVr>^0. 0 
Sl*R6=^0. 0 
GO TO 9020 
9000 FVLM^-2. 

6RVD=-3. 0E-1I: 

SUPG^-‘2 0E-13 
9020 DO 2027 Kt^l..4 
DO 2026 L-1.4 
DO 2025 M- J .. 8 
DO 2024 N=l. 8 
GFMS<L.. L, H. N>=^ FYLM 
IF<M. GT. 7) GO TO 282^ 

IF<M. CiT. 2 OR. N. GT. 2>G0 TO 2820 
6IS<1/ K, L. M. N)=SUPG 
GIS<2.. K, L. M. U'r<iRVD 
GJS*'3. K. L. M.. N )-FVLM 

2020 2022 

GN J. K. L. M. N)*.^SUPG 
J.. K* 1, 11. rO^-iRVD 
1F<M. GT 2 >60 TO 2022 
GflKl. .t/K. I .. H. N>^-;i IPG 
6TU<2. J, K. L.. H.. N>^-Gf 

sT.k. i .11 ro^Fvi m 

T, r:.. I .. H N>=i.UFG 
\ .H.N> iiKVD 

r(»NrjNi> 




A-15 



2624 CONTINUE 

2625 CONTINUE 

2626 CONTINUE 

2627 CONTINUE 
GO TO 2505 

C THIS BLOCK SETS POTU CONTRIbUTIONS P— EQURL TO CORRESPONDING P~Q 
C 

2616 DO 2618 K=l> 4 
DO 2017 L«=l*4 
DO 2616 H=:1^8 
DO 2015 N=l>8 

PFI^<K, L, H. N>=PFHSQ<K^ L/ N> 

IF<H. 6T. 7>G0 TO 2015 
DO 2014 I=l>2 

IF<K GT. 2. OR. N, GT. 2>G0 TO 2011 
PIS< I> K> L> ti* N>«PISQ< 1/ K> L. H, N> 

2011 DO 2012 J=d^2 

IF<M. GT. 2>60 TO 2012 

PTU< I, J> K, H/ N>»-PTUO< I, J. K> L> N> 

2012 IF<I, GT. 2>GO TO 2612 

PM< 1, J, L> % N>=PNQ< I, J> 1C H, N> 

IF<ri GT. 2>G0 TO 2612 

PED< I, J. K. L> M* N>=PEDwa> J* K. L, M. N> 

2012 C«^TINL€: 

2014 CONTINl^ 

2015 CONTI^R.lE 

2016 CONTINUE 

2017 CONTINUE 

2018 CONTINl«=: 

GO TO 2587 

C THIS BLOCK COHPUTES POTENTIAL RT PRRTICLE 
C 

2050 Sl-W^^©. 0 ’ 

DO 2058 K==l,4 
DO 2057 L==^: A 
DO 2V156 M=^l. s 
t>0 2055 N==^l>8 

SUM=^SUIHPFMSD<K. L.. M. N>*GFMS<K. L/ M. N> 

IF<M GT. 7>G0 TO 2855 
DO 2054 I=l>2 

IF<M. GT. 2. OR. N. GT. 2>G0 TO 2851 
SUH=^SUM+PISO< I.. K.. L, M, N>^GIS< I. K, L.. M. N> 

2051 r>0 2052 J=l,2 

IF<M. GT. 2>G0 TO 2052 

5UM-SUIH*PTUQ< I . K. L. M, N>^ GTU< I> J, K.. L, I1> N> 

2052 GT. 2>G0 TO 2852 

SUM=SUM^PWQ< I, J. K, L. M. N>^GW< I. J. K. L, M. N> 

IF<M. GT. 2>GO TO 2052 

SUM=SUr*HPEt»Q^ I > J. JC M. N>=»GED< I, J, K- L> M. N> 

2052 CONTINUE 

2054 CONTINUE 

2055 CONTINUE 

2056 CONTINUE 

2057 CONTINUE 

2058 CONTI riUE 
GO TO 2511 

C 

C THIS Bl.OCK COMPUTES CHRRGh.; ul GRID RND FILM STRIPS 
C 

2100 0Cjl-Ci. 0 
BC2-0. €t 
RC3=0. 0 

ftr;4=.0. 0 
RFl wfi 0 
0 

RF2-‘0. 0 
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flF4«^ 0 
DO 2105 KL^l.4 
DO 2104 
DC 2103 N^l^S 
0F1=«=^1+GFHS<KL. 1> N> 
fiF2«0F2^GFHS<KL> 2^ H. N> 
fiF3=0F3^GFMS<KL/ 3^ M, N> 

|Vr4«0F4^MS<KU 4^ W, N> 

IF<K GT. 7>G0 TO 2103 

IF<ri GT. 2 . OR. N. GT. 2 >GO TO 2101 

ftFl=flFl-K5IS<3, KL, 1, rt N> 

RF2«0F2-Mj1S<3, KL, 2> M, N> 

RF3=fiF3+GIS<3> Kl^ 3* H> N> 

AF4=.rF 4^IS<3. KL> 4, M, N> 
ftCl=^flCl^GIS<2, 1. KL. M, N> 
ftC2-ftC2+GIS<2, 2. KL> H> N> 
ftC3=fC3+GIS<2> 3. KL> N> 
fiC4=ftC4+GIS<2. 4, KL> N> 

2101 IF<H. GT. 2>G0 TO 2102 

fiFl=flFl-K5TU<3, 1> KL. 1> H> N>+GTU<3> 2, 1> KL> M, N> 
ftF2=«^2+GTU<3, 1. KL, 2, H. N>+GTU<3> 2> 2, KU H> N> 
ftf^3==ftF3^GTU<3, 1. KL, 3, N>-K5TU<3> 2> 3. KL> N> 

^:F4=^RF4+GTU<3, 1, Kl.> 4, n N>+GTU<3. 2> 4> KL. H. N> 
ftCl=FlCl-KjTU<2, 1. KL> H. N>+GTU<2. 2^ KL. 1> H> N> 
ftC2==fC2+GTU<2, 1/ 2^ KL> H. H>+GTU<2. 2> KL> 2/ H> N> 
ftC3=^ftC3-H5TU<2, 1. 3^ KL, H/ N>-»OTU<2> 2/ KL, 3. M. N> 
ftC4=RC4+GTU<2, 1> 4^ KL> H. N>-KjTU<2> 2^ KLM/ M, N> . . 

2102 0Cl=^fiCl-KiM<2> 1, 1> KL, H, N>^GM<2, 2. KL. 1 . M/ N> 

«:;2^ftC2+GW<2, 1, 2, KL, H, N>^GH<2. 2> KL. 2. M, N> 
ftC3=^3+GM<2. 1^ 3, KL> N>+GH<2, 2. KL. 3, l*t N> 

RC4=ftC4+Gl4<2. 1> 4. KL. M. N>+GW<2, 2, KL, 4, H/ N> 

IF<M. GT. 2>G0 TO 2103 

ftCl^^fiCi+GETKa, 1> 1> KL. M. N>+6ED<2. 2, KU 1. H. N> 
fiC2==ftC2+6ED<2. 1. 2. KL, M. N>+GED<2. 2^ KL. 2. M, N> 
ftC:r<=^fiC3+GEr <2^ 1. 3, KL, M. N>+GED<2, 2, KL, 3. M. N> 
fiC4=ftC4+GED<2. 1. 4. KL, M. N>+6ED<2, 2. KL. 4. H, N> 

2103 CONTINUE 

2104 COT4T1NUE 

2105 CONTINUE 
GO TO 3512 

4499 C«LL POTCON 
C 

C COMPUTE CHRRGE ON EftCH ELEMENT f>L^ TO REPLIED POTENT I RL RNDPOTENTIRL 

C C0f4TRIBUTI0^45 FROM ALL OTHER ELEMENTS RND PARTICLE NOTE THAT IF PARTICLE 

C IS PRESENT THEN ALL AF'^LIED POTENTIALS ARE ZERO. 

4999 CI4G<l>-0. 0 
CHG< 2 )=^. 0 
CHG<3>s=^e. 0 
DO 2510 K=^l/4 
DO 2517 1-1.4 
I'd 2516 M^l. A 
DO 2515 N-1/0 

CHG<3>-CHQ<3>+<V-PFMS<K. L.. M. N>>^0. 11853F-12 
IFCM. GT. 7>GO TO 2515 
V-^Vi 

DC* 2514 1=^1. 3 
IFCI. EQ. 2>V-V2 
IF<I. Fi:^ 3>\ A*" 

IF<M. GT. 2. OR. N GT. 2 >G 0 TO 2511 
CHO'*. I >'-CHG< 1 >^<V-FISO. K. L. N- N 7111799E-‘13 

25il DO 2513 vI-1. 2 

1F<M. GT. 2>G0 TO 2512 

CHG< I >- CHQ< I >^^<V-^PTU< I. J. K. L.. M.- N> >’♦^0. 9322102E-13 

25:12 JF^I.GT. 2>G0 10 2513 

CHCi< 1 I >^CV~PW': I. J. K, L, I1. N> > 40 . 5140147E-13 

IFvM GT. 2>G0 TO 25J 
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CMO< 1 >«=CHQ< 1 >♦< VHPED< 1^ J# lO U H. N> 4137133E-12 

2513 CONTINUE 

2514 CONTINUE 

2515 CONTINUE 

2516 CONTINUE 

2517 CONTINUE 

2518 CONTINUE 

DO 2538 1^,2 
SC8L£<1>«^1. 8 

1F<CHQ<I>. OT. QHRX<1>>SC8LE<I>»QHA><<I>/CHG<I> 
IF<CHQ<1X LT. 0HIN<I»SC8L£<I>»QH1N<I>/CHG<I> 

2538 CONTINUE 

START MIRES 


5888 

5881 

5882 

5883 
580^ 
5885 


CR1T»8. 8 
V*V1 

DO 5885 I»1^2 
iF<I.Ea 2>V«V2 
DO 5884 J=^l#2 
DO 5083 K=l#4 
DO 5082 t=1^4 
DO 5081 M=l>7 
DO 5880 N^l#8 

TEMP1*=<V-FM<I# J# K> L, N»«8. 5148147E-13 
TEMP1*=TEMP14cSC«LE< I > 

TErff>2=TET«>l^ftCT0R^6N< L> M^ N>»<l-FRCTOR> 

TEMP3'-=PERCEN^GM< I. N> 

TEMP4=TEMP1-GW< 1, 1C U M> N> 

TEMP3«f»S<TEMP3> 

TEMP4^«BS<TEMP4> 

IF<TEMP4. QT. TEMP2>CRIT-1, 88 
GH<I> J, K, L. M, N>-TEHP2 
CONTim^ 

CONTIM^ 

CONTIW.IE 

CONTiraJE 

CONTINUE 

CONTINI.IE 


Em HIRES. STf^T GRID EDGES 


V«V1 

DO 581x 1=^1, 2 
IF<1. EQ. 2>V=^V2 
DO 5010 J-1.2 
DO 5C^ K^^l.4 
DO 5808 L=l. 4 
DO 5007 
DO 5806 N=^l/8 

TEMP1^<V-PED<I, J. K. L, M, N>>^0. 4137133E-13 
TEMPl=^TEHPl^SCaLE< 1 ) 

TFMP2=TPMP1*FRC:T0R+GED<I^ J, K, L. M. N>^<l-F0CTfM^> 
TEMP3^PERCEN*GED<I, J. K. L. 1C N> 

TFMP4=^TEriP1 -GED< I, X K, L, M/ N> 

TEMP3=^ftPS<TErr3> 

TFMP<t==^HPS<TEliP4> 

IF<TEMP4 Gf. TEMP3>CRIT=^1. 0 
GED'C I.. J, I , M, N>- TEMP? 

5086 CONUNMF: 

5007 rONTJ NHV 

5008 C(»rniNni' 

500S^ C.C»r4TTNin* 

5010 CONTlNi*E 
501:1 run UNUE 

C FND Gkll* Fl^GFS SIBRl MAIN GRID/FJLM 1 ^ V 
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v«vi 

DO 5ed7 
IF<1. Ea 2>V«V2 
iF<i, Ea ::>v»v3 
DO 5016 J»1^2 
DO 5015 
DO 5014 L«1^4 
DO 5012 H»l#2 
DO 5012 

TEHP1«<VH>TU<1* 10 U It N»^8, S222102E-12 
TEHPl*TEHPl^SCflLE< I > 

TEMP2^ TEHPl«»FftCTOR^QTU< lO L* It M>^<1-FBCT0R> 
TENP3s=PERCEI^TU<l/ 10 L# It N> 

TEHP4=^TEMP1-GTU<I^ lO L> It N> 

TEHP3=«BS<TEHP2> 

TEMR4*=BBS<TEHP4> 

IF<TEMP4, GT. TEW>2>CRIT*^1. O 
GTU<1> J, lO It N>=^TEHP2 
5012 CW4TINUE 
5012 COMTII«E 

5014 CWTINUE 

5015 COMTIICC 

5016 CONTINUE 

5017 CONTII^C 
C 

C EIC^ HAIN GRID/FILM T 4 U. START INTERSECTION SQUARES 
C 

V=:V1 

DO 5022 1=^1^ 2 
1F<I. EQ. 2>V=V2 
IF<I. EQ. 2>V^V2 
DO 5021 K^l>4 
DO 5020 L=-l>4 
DO 501S «=1>2 
DO 5018 N=l>2 

TE«P3 =-<V-PIS< I . K, L.. rt N> 7111739E«12 

TEMPl^ TEMPI ♦SCALE < I > 

TEMP2^TEMP1 FACTOR+GIS< I. K. L. H. N>^<! -FACTOR> 
TEMP2=PERCFN^CtIS< I, K. L.- H- N> 

TEI^4=^TFI*F>l--CiIS<l^ K.. L, M. N> 

TEMP2-ABS<TEr*r2> 

TFl*P4=^ARS<TEHP4> 

IF<TEMP4. 6T. TEMP2>CRIT=:l. 00 
GISa. K. L. M> H>-TEMP2 

5018 CC4>ITINUF 

5019 CONTII^IE 
502€t CONTINLIF 

5021 CONTINflE 

5022 CONTJNiiE 
C 

C END INTERS! CT ION SQilf^ES. START FILM MAIN SQUARES 
C 

v-vr< 

DO 5026 K=^l,4 
DO 5025 L^l,4 
DO 5»-i24 H-J .8 
DO 502? N^t.8 

IFMP1=.*:V-PFHS 1,1 »M. W> '4A. :4 185?F“12 
TFMP I - fEHP.l4 Si:.rH 
TFMFv: TFMPHPMrT 
TFHP? . PFROrU » >FN 
rFMP4^ IFHPt-MHt-. 

TFMF U Af:S. | j- hk* > 

TFMP4 APS<TFHP4‘» 

I F •: Tf HP4 GT. 1 1- Mf ’I >CP I T^ ;l fir i 


F*'?:> 

nR+GFMS<l5 I . M. N>’* <:1-4“ ACT0R> 
S* K . ) . \h M> 

*:F. I..IL M> 
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sa23 

sa24 

5825 

5826 

8188 
. 8181 
8182 

8183 

8184 
C 

C 

C 

C 

C 

C 

C 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

5999 


6888 


688:1 

6082 

6884 

6885 
6806 
6087 

C 

C 

c 

c 

c 


0FHS<I0 U K> N>«TFMP2 

CONTINUE 

CONTINUE 

CONTINUE 

CONTINUE 

HR1TE<6^ 8ieeXPIS<I^ 2. 2* 2/2>^ OIS<I, 2, 2> 2^ 2>» I«l, 3> 

F0RHftT<3<5X. -^PIS « ^ Ell. 4> 5X> 'GIS « -'>E11. 4/>> 

NR1TE<6.. 181><PTU<I> 2. 2> 2, 2» GTU<I^ 2^ 2^ 2^ 2, 8>^ I«l, 3> 

FORHftTO 5>C ^PTU = Ell. 4> 5X^ -^GTU = -^,E11. 4/>> 

NRITE<6> cfl82XPED<I> 2« 2> 2. 2> 8>. GED<I, 2. 2, 2* 2> 8>^ 1=^1, 2> 

F0RHftT<2<5X> 'PED = Ell. 4, 5X, -'GED = Ell. 4/>> 

HRITE<6. 8183XP«<I^ 2> 2. 2, 7, 8>^ GMCI^ 2, 2^ 2, 7> 8>> 1*1. 2> 

F0RHftT<2<5X. 'PW * Ell. 4. 5X> 'GW * -'.E11.4/>> 

WRITE<6. 8184>PFMS<2. 2. 4. 4>. QFMS<2> 2. 4. 4> 

F0RHBT<5X. 'PFMS ^ '.Ell. 4. 5X. 'GFHS * '.Eli. 4/> 

END FILH HRIN SQUARES 

COHPUTATION OF NEW CHARGE DISTRIKITION OWLETED 
IF CRIT IS NOT ZERO THEN CH^GE OF CHARGE ON AT LEAST ONE ELEMENT 
EXCEEDED SPECIFIED PERCENTAGES. IF CRIT IS ZERO THEN ITERATION HAS 
REACHED REQUIRED ACJClff^V 

COUNT*COUNT-H 

IFCCOiJHT. LT. CVCLES--0. 1. AND. CRIT. GT. 8. 5>G0 TO 3588 
CONDITIONAL RETURN TO START OF ITERATION PROCESS 

WHEN ITFRATIOtil COI^ETED PERMITTED UUmER OF CVCLES REACI«D 

COMPUTE POTENTIAl. AT POSITIOf^ OF PARTICLE W>E TO APPLIED POTENTIALS 
AND PARTICLE IMACiE Ct«^S 
GO TO 3589 

CCWPUTE POTENTIW. AT ELEMENTS OkfE TO ONE CCHJLOI« AT *=^TICLE 
START PARTICLE ON WIRES 


2=:2P-a 8097828 
DO 6887 

IF<1. EQ. 2>2^2P-8. 8€t67564 

DO 6086 J=^1.2 

DO 6085 K=^1.4 

DO 6084 L=-l. 4 

DO 6083 M=l. 7 

DO 6002 N=l. 8 

IF<J. EO. 2>G0 TO 6800 

A= 0. 0 

B=.0. 083175 

X^0. 02921^K+0. 003175^M~0. 085725-XP 
V=«0. 02921 *L+8. 003175>*=N“0. 087313-VP 
GO TO 6801 
A^e. 003J7T> 

B^8 0 

X--8. 02921 ♦L+8. 8031 75:*^N-0. 887.^1 3-XP 
V= 8. 0292.1 4:K+a C^317r.:fH“0. u85?25~»’*P 
CALL R FTNFO:.. V.::. A, A. PDTI > 

PWQ<J. .T,KM ,M.N>=^POri 

CONTINUF 

CW4TJNUF 

CONTINUE 

CONTINI IF 

CON NUF 

CON I NUF 

END PARTICl.E ON WIRFS 
START PARTICLE ON GRID EDGES 
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DO 6015 1*^/2 
ir<l. EG. 2>Z»ZF-8. 0067564 
DO 6014 J«1^2 
DO 6013 K»l#4 
DO 6012 LsdU4 
DO 6011 H=1^2 
. DO 6010 H«l>8 
IF<J. £0. 2>Q0 TO 6008 
RbO. O 
B»a 003175 

K*=0. 02921>HC^. 0254««-0. 111125-XP 
V=0. 02921*i«L-^e. ee3175^*tl-0. 087313-VP 
GO TO 6009 

6008 003175 
0 

X=0. B292±*L^. 003175^N-e. 087313-XP 
V«8. 02921»K^a 0254^M«O. 111125-VP 

6009 COLL PLEINF<X/ V> 2, fi, POTL> 

PEDQ<I^ L, n. N>~POTL 

6010 CONTIM^ 

6011 COHTIKUE 

6012 CONTimiE 
€013 CONTir^C 

6014 CONTINUE 

6015 CONTINUE 
C 

C END PORTICLE ON GRID EDGES 
C 

C STftRT PARTICLE ON HRIH GRID/FILM T It U 
C 

2=ZP-0. 0097028 

DO 6023 1=1.3 

IF<I. EQ. 2>Z=ZP-0. ^67564 

IF<I. EQ. 3>2=2P 

DO 6022 J=1.2 

DO 6021 K=l>4 

DO 6020 L=l>4 

DO 6019 M=1.2 

DO 6018 N=1.8 

IF<J. EQ. 2>60 TO 6016 

A=0. 001905 

B=0. 003175 

X=-0. 02921’«^K+0. 027305*:M-0. 1139825-XP 
V^0. 0292i^L+0. ^3175^N-0. 8873125-VP 
GO TO 6017 

6016 A=0. 003175 
B=0. 001905 

X^0. 02921^L-^0. 003175>i44-0. 0873125->^ 

V=0. 0292j>i*K+0. 027305^H-0. ..139825-VP 

6017 Cftt.L Pl..Eir#^<X. V Z. A.B.POTL> 

PTUQ< I. J. K. L. M. M>-POTL 

6018 cot^Tir^r 

6019 CONTINUE 

6020 CONTINUE 

6021 CONTINUE 

6022 CONTINUE 

6023 CONTINUE 
C 

C END Ff^TlCLE ON MAIN C>RID/FILI1 T It U 

C 

C START PARTia.F: ON GRID/FllM INTERSECTION SQUARES 
w 

0097020 
DO 6020 1-1.3 
IF<I.FO 2>r -7P-(t 0067564 


A-21 



B$R 4234 


1F<I. EG. 3)Z«ZF> 

PC eeS7 K«i«4 
PC 6026 Lk1«4 
PC 6025 2 

DO 6024 N*l-2 
(M. OOIPOS 
B«e. 001905 

X*=0. 02921>*4C+0. 027305'*4*-0. 1139825-XP 
V*fO. 02921*L-*O. 0273e®-*M-O. 1139825-VP 
CflIX PLE1HF<X> V# 2/ ft, B, POTL> 

PISG< I, K, L, H, N>“POTL 

6024 CONTINtC 

6025 CONTIWJE 

6026 CONTimJE 

6027 CONTI»AJC 

6028 CONTINIE 

END P«?Tia.E ON GRID/FILH INTERSECTION SOUftRES 

STftRT Pffl?TlCLE FILM »«IN SQU«^ ELEMBITS 
Z“ZP 

DO 6022 K*=l,4 
DO 6021 L==l>4 
DO 6020 M=l, 8 
DO 6029 N=l,8 
ft*0. 002175 
B=0. 002175 

Xa0. 02921>H<'*-0. e©2175*H-0. ^72125-XP 
V==0. 02921*L-^0. 002175*N-0. 087212S-VP 
CfiLL PLEINF<X, V,S, fLB,ROTl.> 

PFHSQCK, L, H, N>=POTL 
-0^ CONTINUE 
. s)20 C0NT1NI« 

6021 CONTINUE 
60~2 CONTINUE 


C 

C 

C 

C 


C 

C END PftRTlCl.F ON FILM NOIN J^OUftRI-S ELEMENTS 
C 

C fill. POTENTIM-S IH»=^ TO CI¥iFGE CiN P«?TICLE COMPLETED. NOTE THftT PMQ, PEDQ, 

C PTUO,PISO #! PFMSQ ftRE HLSO TW! INFLUENCE COFFFICIENTS Ff« COMPUTING 

C POTENTIfH. AT PARTICLE DI.IF TO CHARGE DISTRIBUITCMT ON ELECTROl-FS 
C 

fiO TO 2504 
2514 CONTINUE 
END 
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SUBROUTINE POTCON 

DIHENSION CI^<9€7^>>CFHSFU<5194>,CFTUFN<21364>>CFeDQM<lM82 
DIMENSION CFPFIM<1097S>> CFFHFH<28eS>> a=TUFH<89B4>. CFEDFH<5S%> 
DIMENSION CFISFM<94e8>. CFTUTU<18595>^ CFISIS<€05># CFEDED<7436> 
DIMENSION CFEDIS<3136>^ CFEDry<15512>> CFTU1S<24S4> 

COMMON GFMS<4> 4, 8# 8>> PF.tS<4/ 4, S. 8> 

COMMON PM<2> 2^ 4> 7# 8># QN<2# 2. 4. 4« 7# 8> 

COMMON Q£D<2« 2# 4> 4« 2> 8>. FED<2^ 2/ 4^ 4/ 2# 8> 

COMMON 0TU<2^ 2^ 4^ 4> 2> 8>^ FTU<3^ 2/ 4M. 2^ 8> 

COMMON QIS<3^4^4^2«2>>PIS<3#4.4^2^2> 

E€ailVRLENCE<CFTUFM<l>, CFTUTU<±>^ OnUFH<l>^ CFPF1M<1># CFm<±>, CFISIS 
2<1>> CFMSFM<1)> CFFMFM<1>. CFISFM<1>. CFEDIS<1>^ CFEDTU<1>> CFTU1S<1>>. < 
2CFEDGN<1>> C^^EDFM<1>> CFEDED<1>> 

C 

C ST8RT INTERSECTION SQLiftRES ON INTERSECTION SQU^S 
C 

R£BD<28>CFISIS 
DO 4S48 LL»1.4 
DO 4S47 NN^1>2 
II^=“3<«LL+W4 
DO 4646 MM»1.2 
DO 4645 KK=1.4 
IMP2=^2HtKK-^W1 
DO 4644 M=1^2 
IND2=^M-*IMP2 
DO 4643 K=l,4 
IND2=I>«>2+3 
ir^3=^If»S<ir«)2> 

IND4=^-IW^1 
DO 4642 L=l>4 
IND4=ir«>4^3 
IND6=^ll>t: If»S< IND4+1> 

IND7=11^IRB5< IND4+2> 

INXl=ir«>6+IND3 

JNX2=^IND?+IND3 

r>0 4641 11=^1. 3 

TEMP13==0. 0 

DO 4640 1=1,3 

IF<I. EQ. II>GO TO 4637 

IF<I. EO. 3. OR. II. EO. 3>GO TO 4638 

INDl^-243 

GO TO 4639 

4637 IND1=1 

GO TO 4639 

4638 IF<I-m. EQ. 4>IND1==364 
IF<I + II. EQ, 5>1ND1=485 

4639 ir4DFXl=IN3^:i+INi »:1 
INDEX2=INX2-HND1 

THI^rO^TEMPO-i-CFISIS* 1NDEX1>^QIS< I, K, L, K l>^CFISIS<ir4DEX2>4:GIS< 1/ K 
?, L, M, 2> 

4640 CONTI^4UF 

PIS< 1 1, KK.. LL, HM. NN>=PIS< 1 1, KK, LL. MM, NN>+1 FI1P13 

4641 CONTI t4UE 

4642 CONTINUE 

4643 CONTINUE 

4644 continue 

4645 CONTINUE 

4646 CCiNTINUE 
47 COimNUE 

4648 f ONTir4UE 
C 

C F.HO INTERSECTION SOUftRES ON INTERSECTION SQUGReC 
C 
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C 

C stmt HIKES ON MIRES 

C 

f^IND 20 
REM><i6>CFMH 
DO 4516 NN»1>8 
DO 4515 LL»1#4 
1ND1-15*LL<H«4 
TEMP5«^44iLL-iNN-0. 5 
DO 4514 MM»1#7 
DO 4513 
IND2-i3>N<K<H1M 
TEMP6s=14»KK4NrH0. 5 
HOLD1«0. 0 
H0LD2^. 0 
M0LD3««. 0 
H0LD4^. e 
TEHP17= IND2*-TEW^ 

DO 4512 N=l>8 
ir«>3s:4f-lNDl 
TEMP7*=N-TEHP6 
DO 4511 L«1.4 
lND3=ir«>3-H5 
1ND5^464>IBBS< IND3>^1 
TCMP7=TEMP7^14 
TEf«>ls=ft8S<TEMP7>+0. 7 
IND7“ INT<TEMP1>*M876 
DO 4510 M“l,7 
lND4s^H-ir«>2 
DO 4509 K=1^4 
IND4-ir6>4^13 
irc»£Xi-i0BS< ir«>4>*«^ira>5 
TEMP8=I ND4^TEf^l7-^K 
TEMP8=fBS<TEJ^8>-0. 3 
II«>8=49^1NT<TEr«P8> 
lH0EX3^im7^1HT>B 
I NDEX2= I WDEXl+2438 
I NDEX4= 1 NDEX3+2401 
TEMP±2^CFWW < I NDEXl > 

TEHP9:^-CFNN< INDEX2> 

TEMf>0=-CFWM< INDEX3> 

TEMP10=:CFI^^<: INDEX4> 

TEMP21=GN<1, 1, K. L. M. N) 

TEMF22=6W<1. 2, K, L. M> 

TEMP23=^GW<2/ 1. K. L. M. N> 

TEI1P24==GW<2/ 2/ K. L, M/ N> 

H0!.Dl-H0LDl-»-TEMP9^TEMP23‘fTEI1P12^TEI1P21-i'TEMP10*TEI1P24^-TEHPi3^TEr*r22 
H0LD2- H0LD2-MFNP9*^ fEMP2 J +TF.MP:1 2^T EhP2:^+ "T EHP10*^TFMP22'»-TEI1F‘i2^T£t1P24 
MOt.D^^HOL r>3-MF11P9f»'TEf1P24+TFHP:1 2^1 EMP22+TEMPi 04^TEMP22>TEMPiI:^TEMP21 
D4- H0LD4->'TFMP9’t'TEMP22^-TEMP12*1 EMP24+TEMPlCi=t 7 EMP21-M EMP12:’*^TEHP23 

4509 CONTINUE 

4510 CONTINUE 

4511 CONTINUE 

4512 CONTINUE 

PW< J / 1. KK, LL. MM. WN>=-PW<1.. 1. KK. LL. MM. NN>^HOLDi 
PW<2. 1. KK. LL. MM. NN>--«PU<2. 1. KK. LL. MM. NN>-t*H0LD2 
PH<1. 2. KK. LL. MM. NN>“PN<1. 2. KK. LL. MM. NN>-i*HOLr):< 

PW<2. 2. KK/ LL. MM. Nr^>^PN<2. 2. KK. LL. MM. mn>+H0LU4 

4512 CONTINUE 

4514 CONTINUE 

4515 CONTINUE 

451n CONTINUE 

C: 

C HIRES ON WIRES FINISHED 

C 

C STRRT F1U1 MHJN SQURRE ON WIRES ftND VICE VERSR 
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RBHIND le 
REM><U>CFn5FM 
DO 4531 NN«1>8 
DO 4539 LL-1.4 
lNC«»15>t>LI.«NN 
DO 4^9 HH-1.7 
. DO 4528 KK«1>4 
TEHPlel4*ia<^HH+0. 5 
DO 4527 JJel>2 
TEHP12«GH<1> J J< IOC LL> tttl. 

TEMP14cGH<2> J J. KIC LL« MM> »M> 

TEMP13=0. 0 
TEHPie^. 0 

IF<JJ. EQ. 2>e0 TO 4521 
TEHP4=-TEHP1 
DO 4520 K-1,4 
TEW»4*14+TEMP4 
DO 4519 M=i>8 
TEHP2=TEMP4+M 
TEI«l>2=fiBS<TEHP2>+0. 7 
IND2--INT<TEMP2) 

IND3--ir®6 
DO 4518 L=l,4 
ll«)3^IND3+15 
DO 4517 N=i,8 
1M>4“IND3+N 

lNDEXl==^49*If«!S< IM>4)-*-IND2 
INDEK2=INDEXl+2597 
TEHP3=CFHSFW< I Nt«Xl > 

TEMP5=CFHSFM< INDEX2> 

TEh»»13=TEHP13+GFMS<K^ L, H, N>>•=TE^tf>3 
TEHP16=TEI«>16+GFHS<K> L, M, N>>«'TEHP5 

PFMS<K> L/ H/ N>»=PFHS<K, L, M. W>+TEMP3»TEW12+TEMP5>»=TEI*'14 

4517 CONTirflJE 

4518 CONTINUE 

4519 CONTINUE 

4520 CONTINUE 
60 TO 4526 

4521 TEI-»>4=-TEI«»l 
DO 4525 L=l, 4 
TEMP4=TEMP4+14 
DO 4524 N=l,8 
TEMP2=TEt1P4+N 
TEMP2-RBS<TEtlP2>+e. 7 
IND2-INT<TEMP2> 
lh«:>3-- -IND6 

DO 4523 K-1,4 
IND3=IND3+15 
DO 4522 M- 4,8 
lND4=rIND3+M 

lNDEXl=-49»IfifeS< IND4>+IND2 
INDEK2- INr4EXl+2597 
TEMP3-CFMSFW< INDEXl) 

TENP5=^CFMSFW< TNr<EX2> 

TEMP;1 3-TEHPi:;.+GFMS<IC L, M, N>*TFMP3 
TEMP16«TFMP16+tiFM?:<.K, L, t1, N>^=1 FMP5 

PFMS<K, L, M, N>= PFMS<K, L, M, M>+TEMP3^TEMPJ.24-TEMP5*TEMP14 

4522 CONTINI.E 

4523 CfiNTINUE 

4524 CONTINUE 

4525 CONTINUE 

4526 PIKI, JJ, KK, LI.., MM, NNT=^PW*:i, JJ, KK, l.L, MM, NN>+TEMPJ 3 
PW<2, JJ, KK, LL, MM, NN>»=PN<2, JJ, KK, l.L, MM, NN>+TEMP16 

4527 (ONTINUF 

4528 r ONTINUF 
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4S2S CONTlNliE 

4530 CONTlf^ 

4531 CONTIS 

END FILH HNIN SQUARES ON MIRES AND VICE 
START QRID/FILM MAIN T 4 U CM HIRES VICE VERSA 


STARl LATE EDGES ON MIRES AND VIC^ VERSA 

REMir«> 11 
READ<12>CFTiFM 
REMIND 12 
REf^<13>CFEDGM 
13 

DO 4548 W4=1.8 
DO 4547 LL«1#4 
IW>1«15^LL+NN 
TEW>1^16>M.L+2*NN-7, 5 
DO 4546 MH-1,7 
DO 4545 KK-1,4 
TEW>2*^14^KK+2*»‘riM-6. 5 
TEMP3=14»H<K+WH0. 5 
TEMP13=^GI4<1/ 1, KK^ LL, m, NN> 

TEMP14=:GW<1. 2/ KK# LU MM. NN> 

TEI*P18=^6W<2/ 1> KK. LL, MM, 

TEMP19=GW<2, 2, KK> LL, MM/ NN> 

TEMP12^0. 0 
TEMP15=0, 0 
TEHP16=^0. 0 
TEMP17==0. 0 
HODl^O. 0 
HOD2-^0. 0 
HOD3=0. 0 
MOD4=^0. 0 
DO 4544 N=l/8 
IND2^N-IND1 
TEMP4=N-TEMP3 
t>0 4543 L==l/4 
IND2=IND2i-15 
IND3'-=49^IABS'C IND2> 

TEMP4:^-TEMP4^14 
TEMP5r.ftBS<TFMP4> 

INt>4=riNT<TEMP5+0. 7> 

DO 4542 M==l,2 
TEMP6-4‘1-TEMP2 
TEMP7==M-TEHP1 
DO 4541 K:-l/4 
TEMP6^‘1 EMP6+14 
TEMP7==TEMP7+16 
TEMP8=-RPS<TEMP7>-0. 2 
I ND8-49* I NT < TEMP8 > 

TEMP9*^ ABS<TEMP6>+0. 7 
IND9==^ir4TaFMP9> 

INX1>:fMD3+IND9 
INX2-^INC*«+IND4 
7000 INDEX J/ INX1 

1 NDEX2s^ I NX2+:l 03BB 
JNDEX:i«^ INXJt-2597 
JNDEX4^ 3NX2+i:i.j:<2 
INI)FX'>“ UM2+5194 
INDEX6= INX2+79.S0 
TEMP3B^ ( f FD(jH< index I > 

1 EMP37^ *: FH.DCiN< irn.»rx5> 

T EMP3B^ < F F C 'CiN < I Nt 'feXIi > 
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TEHP39»CFEDGM< INDEX€> 

TEHP30>6ED<1/ 1* L, N> 

TEIif>31»QED<l> K. L> ft N> 

TEfff>32>>QED<2> 1/ 1C U ft N> 

TEfff>33>QED<2» 2> tC L» ft N> 

HODl»HODl+TEffl»36*TEHP30+TEf1P37«i‘TEHP3l4TEt1P38*TEf1P32‘«-TEt1P39>i*TEf1P33 
HOD2-HOD2+TEfB»3«%TEHP31*TEHP37*TEMPS0+TEMP38*TEf«»33+TEf«»39*TEf«»32 
H0D3*H0O3*TEHP36>KTE«P32+TEf«»37*TEMP33+TEHP38*TEMi»3e+TB1P39*TEMP3i 
HOD4«HOD4+TEHP36»«‘TEf1P33+TEf^37#TEMP32+TEHP38*TEW>31'HEMI»39*TEf1P30 
PED< jU ±. K, L/ ft N>«PED<1. K/ U M> N>-*-TEMP36fTEMP13+TEHP37.*TEI1Pi4-*-TEfff> 
238*TEf»>18+TEHP39*TEHP19 

PED<1> 2/ K, L, M, N>»PED<1, 2. K, L, It N>^TEHP38*TEMP14+TEHP37*TEI1P13+TEMP 
338*TEMP19+TEMP38*TEMP18 

PED<2, 1/ lO L, M/ N>»PED<2, ±, K. L, Mi N>+TEMP36*TEftf»18+TEMP37'«TEMP18+TEfff» 
438*TEMP13+TEW»39t«TEMP14 

PED(2i 2i Ki L> Mi N>«PO><2i 2i Ki Li Mi N>-^TEMP36#TEMP19-*-7EW»37*TEW^18+TEI«» 
538*TEMP14+TEMP39*TEf»»13 
DO 4539 I=li3 
1F<I. EQ. 2>60 TO 7810 
1F<I. EQ. 3>Q0 TO 7828 
80 TO 7838 

7810 INDEXl=lNXl+2597 
INDEX2«INX2+13132 
Jf®EX3»lNXl 
I f4DEX4s I MX2-I-18388 
80 TO 7838 

7020 INDEX1=INX1+5194 
I NDEX2» I NX2fl5876 
lf«>EX3*=INXl+7791 
1NDEX4-I NX2418628 
7030 TEI1P20=CFTUFW<1NDEX3> 

TEMP21«CFTUFN< 1NDEX4) 

TEMP10-CFTUFM < I NDEXl ) 

TEMP11<^CFTUFN< I NDEX2 > 

TEMP22=QTU< Mi Ki Li Mi N> 

TEMP23^ QTU< li 2i Ki Li Mi N> 

TEMP12uTEMP12+TEMP10*TEMP22+TEMPll*TEMP23 
TEMP15®TEMF t 5+TEMP10*TEMP23+TEMPll*TEMP22 
TEMP16“TEMP16+TEMP20*TEMP22+TEMP21*TEMP23 
TEMP17=TEMPi.7+TEMP20*TEMr23+TEHP21*TEMP22 

PTU< li li Ki Li Mi N>=PTU< li li Ki Li Mi N>+TEMP10*TEMP13+TEMP11#TEMP14+TEMP 
220*TEMP1S+TEMP2.1*TEMP19 

PTU< li 2i Ki Li Mi N>«PTU< li 2i Ki Li Mi N>+TEMPll*TEMP13+TEMP18*TEMPi4+TEMP 
321*TEMP18'*-TEMP2e*TEMP19 
4539 CONTINUE 

4541 CONTINUE 

4542 CONTINUE 

4543 CONTINUE 

4544 CONTINUE 

PW<li li KKi LLi MMi NN>=PW<li li KKi LL. .'tt NN>+TEMP12+H0D1 
PW<li 2i KKi LLi MMi NN>«PW':ii 2i KKi LLi lilt NN>+TF.MP15+H002 
PN<2i li KKi LLi MMi NN)-PW<2i li KKi LLi MMi NN>+lEMPlt.+Ha(>3 
PN<2i 2i KKi LLi MMi NN>-PW<2i 2i KKi LLi MMi NN>+TEMP17+H0D4 

4545 CONTINUE 

4546 CONTINUE 

4547 CONTINUE 

4548 CONTINUE 

END GRID/FILH T & U ON WIRES ftND VICE VERSA 

END PLATE EDGES ON WIRES AND VICE VERSA 

START PLATE/'FILM INTERSECTION ON WIRES AND VICEVEP.SA 

RFAD<14>CFPFIN 
DO 4575 NN-liS 
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D0 4574 LL-1.4 
TEHPl=16^^2%NN-7. 5 
DO 4573 MK»1^7 
DO 4572 H3<>^l/4 
TEHP2«14^K:K^2»^MH-6. 5 

1/ KK> LL# «H. flN> 

TEMP12»0M<1# 2. KK# Ll^ MM> NN> 

iE¥fp±3^m<2 , too LL# rat rai> 

TEHP14*«GW<2/ 2^ KK. LL/ rai. HN> 

TE«P15»0, e 
TEHP16ea e 
TE^ff»17«0. 0 
TEra^8=0. 0 
DO 4571 M»l#2 
TEra»3«M-TEMP2 
DO 4570 K*=1^4 
TEHP3»TEHP3-i-14 
TE«P4s^ftBS<TEMP3>^0. 7 
ira>4^==iwT<TEra»4> 

DO 4569 N=l>2 
TEt1P5*^N-TEra>l 
DO 4568 L«1.4 
TEt1P5=TEMf=^-t*16 
TEMP6=RBS<TEMP5>-0. 2 
I ND6==49*«‘ I NT < TEHP6 > 

INX«IND4-^IND6 
DO 4567 1=1/3 
1F<I. EQ. 1>G0 TO 4566 
IF<I. EQ. 2>G0 TO 4576 
IF<I. EQ. 2>GO TO 4577 

4566 INDEXi:=INX 
INDEX2=INX+2744 
60 TO 7040 

4576 INDEX1=-INX^*2744 
INDEX2=-INX 

GO TO 7040 

4577 INDEXl lNX-»-5488 
INDEX2-1NX-^S232 

7040 1 EriP7s-CFPF I W < 1 NOEXi ) 

TEMPOS CFPF IW< INDFX2) 

TEMP15^.TEt1P15+GIS< 1/ K. L/ M/ N>*TEMP7 
TEMP16-TEMP16+GIS<I. L. K/ N. M>^TEMP7 
TEMP17~TEI1P17-i-6IS< ) . K. L. K N>’*1 FMPS 
TEMPie=^TEHPi8+GI5<l/ L. K/ th ri>>*TEMP8 

PIS< 1 / K. L. M/ N>=PISO/ K/ L/ M. N>+TEMP7^TEMPll+TEriP8%TEMP13 
PIS< 1. 1 . K/ N/ M>=PIS^ I . L/ K N. M>+TEMP7>^TEriP12+TF.HP8*TEMP14 

4567 CONTINUE 

4568 CONTINUE 

4569 CONTINUE 

4570 CnNTTNUE 
457 J CONTINUE 

PN<1/ 1/ 1:>/ LL. MM/ NN>=PW<1/ 1/ KK/ LL/ MM/ NNf-MEMPiS 
PlKl/ 2> Kk, LL/ MM/ NN >=PW< j / ?/ LK/ LL/ MM, NN>+TF.MPi6 
PN<2/ 1/ KL/ LL/ MM/ NN *^PN<2/ 1/ KK/ LL/ MM/ NN>+TFMP17 
PN<2/ 2/ KK. LL/ MM r4N.*^PN<2/ 2/ KK/ LL/ MM/ NN>-»-TEMPlB 

4572 CONTINUE 

4573 CONTINUE 

4574 CONTINUE 

4575 CONIINUK 
C 

C FI4D PI.rtTE.TILM INTERSECTION ON HIRES PND VICE VERSft 
C 

C STRkT Fit II MHIN SOUBRE ON f“ II M MB IN SQUBPE 
C 


kFWH4f> 14 
M-.BImLIGm EFMKM 
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DO 4S85 KK>1<4 
DO 4584 MH-1/8 
tf«>l»=15*KK+HM 
DO 4583 LL»1>4 
DO 4582 NN*l/8 
IND2al5*LL+NN 
TEHP13»0. 0 
DO 4581 M=«l»8 
1ND3«H-1ND1 
DO 4588 K«l>4 
IND3«^^IND3+15 
INr;4=IftBS<lND3)+l 
DO 457S» N="i 8 
Ihl>0>=N-IN02 
DO 4578 L-1/4 
1ND5«IND5+1S 

INDEX=^53>«'IflBS< IND5>+IND4 

TEM? L . .*FMFM<lNDEX>>»=GrMS<K, L. M. N> 

TEM^ l - ''tPO+TEMPl 
4578 COirrNt- 
4578 f * "^I' -'E 

4580 CG1.\” 'U'^ 

4581 COtriNUE 

PFI1S<KK, LL m, NN>“PFi1S<KK. LL, MM, NN>-rTF.MP13 

4582 CONTINUE 

4583 CONTINUE 

4584 CONTINUE 

4585 CONTINUE 

END FILM M81N SOUftEE ON FILM MR IN SOURRE 

STRRT OR ID/FILM MR IN T & U ON FILM HR IN SOURRE RND VICE VERSR 


STRRT GRID EDOFS ON FILM MRIU RCURRE ELEMENTS RND VICE VERSR 


REWIND 15 
RERD<16>CFTUFM 
REMIND 16 
RERrM'17'^CFEr.X-M 
REWIND 1? 

DO 4586 KK=1, 4 
DO 4595 MM=^:l, 8 
TEMPI J 6>t KK+P'+MM-T'. 5 
DO 4594 l.L- .1,4 
DO 4r-9.? NU=l.-8 
INr)l=-15+LL+NN 
TEMPI 2-n'=' 0 

TFMPi:i--GFMS<KK, LI.., MM, NH> 
TFMP;t4-fi. 0 

TEMPI.G^GFMS'iLL, KL. NN, MM> 
DO 45 V N -1,8 
lNL-2- N-INIM 
DO 45:-' I I.: J, 4 
]ND?= lUI<,'+lt. 

INt'"- L.*-.*lRH^.vTNrv’.' 

DO 45:-'it M-.1.,' 

TFMP;- 11-TFMH 
DO 45iV K -.L. 4 
I FMP2-- n-.HP; '^'6 
TFMP'.' RP^'<TFMP.';:.’-H\ 7 
iNL"i=- iN'i •:TFHp:r> 

1NX-/Nr.“; + IND4 
ItH >KXJ INX 
JNI'FX2- in: ♦, 

judi-xt.'^ih:' - 6 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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1EMP4«CFTUFM< 1NDEX2> 

TENP«sCFEDrM< INDEX1> 

TEMP7««ra>FM< INDEX2> 

TEHP12<:TEHP12^TEHP6>^iED<l, K. L# K» N>«TEHP7«GED<2« ju K. L* n» N> 
TEHPl^sTEMP14+TEHP<5«QED<l« 2, K L« tt N>«TEHP7«GED<2* 2* 1C U fC N> 

PED<1« 1, K.- L> H. N>sPED<l^ K, U H, N>-i-TEHP6«TEHP13 
PED<2» 1^ IC U H, N>>#>ED<2« 1. K> L# tC ll>^TEHP7*TEHP13 
PED<1, 2, 1C L« tt N>«f>EI><l« 2« 1C U 1C «>-^TEHP«*TEIff>15 
F*EI><2> 2. 1C L, K- N>«PED<2, ICI^ H, N><«-TEHP7*TEHP15 
DO 4588 I»1^3 

1F<1. EQ. 2>TEHP4«CFTUFH<INDEX3> 

IF<I. ECL 3>TEHP4sCFTUFn<IM)EXl> 

TEMPISMSTUO, 1C L> H. ri> 

TEMP1?*GTU< I, 2> 1C L. 1C N> 

TE«P12=TEi •i.2+TEMP4»7EHri6 
TEMP14*TEHP14+TEHP4»TE«F17 

FTU< I> 1C U 1C N>sPTU< 1, 1> K, »C N>+TEHP4*TEHP13 
PTU<1, 2, 1C 1C N)=PTU<1, 2> 1C L> »C N>+TEMP4fTEMP15 

4588 CONTINUE 

4589 CONTI»«JE 
4598 CCNTINUE 

4591 CONTINUE 

4592 CONTI^^e 

PEMS<KK. LL, HTC NN>=PF;JS<KIC LI - MH- NN>-*-TTM»12 
PFHS<LL . KK. '4C MM>-' ' <LL, K;., NN, HH)+TEN>14 

4593 CONTINIE 

4594 CnNTir«.lE 

4595 CONTII^IE 
459t C<»4T INLIE 

C 

C END GRiD/FlLM H8IN T 4c U ON FILM HfilN SQUARE fiW> VICE VERSA 
C 

C END GRID EDGEi. ON FILM MAIN SQLlf#?E ELEMENTS ATO VICE VERSA 
C 

C START intersection SQUARES ON FILM MAIN SQUARE ELEMENTS AND VICE VERS 
C 

REAr-:i8>CFTSFM 
DO 46it. KK=1, 4 
TEMPI- lS^IK-7. r- 
DO 4615 MM=1,8 
TEMPI = TEMPI +?• 

DO 4614 LL-1,4 
TEMP2=-16’H-L-7. 5 
O 4613 NN=:l, 8 
fEMP2-TEMP2*2 
TEMP13-=P 0 

TEMPJ.2“GFMS<KK, I .L.. MM, I#I> 

I'D 46J2 M-1,2 
TEttf 3-- M-TEMP1 
DO 4611 K- 3 , •! 

TEMR-3= TFMPl+1.6 
1 EMP'i ^APS T FMF .Of O. 7 
IN»j4-^INI<TH1P4> 

DO 4610 N=-:l ,2 
Tfc, .'-•r.--N-iFMP2 
DO 4609 1-34 
TFMPI. TFMP‘-.f,16 
IFMP6 AS-' IFNI '.', .-t., < 

IND6 ''.6« INI • ri 
1 NI 4 iNr*4+iiii . 

ini.f:;i- .i 

ini'f::'-- iNi'F::ifT:.'7.-- 
IHH-;- I ■ f F'lSFh', INflViJ . 

IF HP.-. -rFI-.FN-ONl'F;:.... 

u-.iif:.'- pisfh' indfx.'> 

TniP.I 1 H'- ’ ;• fIFi IF 3 I f01S<. ■ M. N +1F NPl- AI'm. I • I ... L- M, FCOflFF PI?. 
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>«IS<2.IGL*H»N> 

PIS<1, 10 L* N. N>s^IS<JO lO U H. N>vrENP:<2*TEMP12 
PIS<2, 10 lb N>sPIS<2» lO U H. N>-«-TEI«»33*TEHP12 
PIS<3i# 10 U n. H>sPIS<3* lO L> Ml N>-»TENP31«TEMP12 
4€a» COHTINUE 
CONTINUC 
A€±± CONTINUE 
4«12 GONTXNUE 

PPHS<tao UL# Mb NM>*PFnS<KIO LU Wb NN>-«TEMP13 

4613 CONTINUE 

4614 CONTINUE 

4615 CONTINUE 

4616 CONTINl- 

END INTBRSEICTION SQURRES ON FILH HRIN SQUARES fM> VICE VERSA 
START PLATE/FILn T 4 U ON PUITEi^ILH T 4 U 


START GRID EDGES ON GRID EDGES 

REMIND 18 
READ<19>CFTUTU 
REMIND 19 
RfcAD<21 >CFEDED 
REMIND 21 
DO 4636 

TEMP1:^-<M»*7. 5> 

DO 4635 KK= 1 , 4 

TB«L=TEMP1+16 

IND 1 = 3 »KK-»MM 

DO 4634 NN ^>8 

TEMP 2 =- 2 *NM- 7 . 5 

DO 4633 LL=J .4 

TENP 2 =^TEMP 2+16 

IND 2 = 15 HL-MN 

TEHP 13 =^. e 

TEMP 25 =©. e 

TEMP 33 =^e. 0 

TEM» 43 = 0 . 0 

DO 4632 H=l ,2 

TEMP 3 = 41 -TEMP 2 

IND 3 =M-IND 1 

DO 4631 K=l ,4 

TEMR 3 =TEHP 3+16 

TF ■' »=rtBS<TEIE 3 >- 0 . 3 

IN 1 ^ 10 :^ 56 *INT(TEMP 4>+2915 

IND 3 =^INC» 3+3 . 

IM> 4=1 ABS< I»ir> 3 >-M 
DO 4630 N“l, 8 
TEMP 5 =-<TEMP 1 + 2 »N> 

IMt» 5 ^N-ir «>2 
DO 4629 L--l,4 
TEMP5^TEMP5+I6 
TEMF6^h8S<TEHF-5>+0 7 
IW‘ 6 - INT<TEHF>.> 

ir4D7 - I J 1 Nl/ ;> 

lNI«8-lND4+jrfl'V 

INI*EX1-IND8 
TNDEX2 IWD9+1.166 
INDEX’ - IttttEXl +':.83 
T ftt'E X4 1 r4r^EX.'+’.<l i 6 
TFMPV=-fiFI.*< J .. A , I- . L, M, N> 
TFMEf:-OED<l, 2, K, I./ M. N’* 
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TEHP9MiF <2#i#IOL^H#N> 

TEHPlfMiED<2^ 2^ IG L# N> 

TEKi»CFEDEI>< IND£X1> 

TEX2«CFCDED< INDE>C2> 

TE)a^CFEDED< IND£X3> 

TCX4>CFEDED< INDEX4> 

TEHP13»^TEMP13+TEXl^TEHP7^TEX2<tTEMPS^TEX3*TEMP^TEX4^TEHPie 

TEMP22«TEHP22^TEXl*TEMP8^TEX2*TEHP7+TEX3>itTEHPie^TEX4»TEHP9 

TCMP33*-TEMP33^TEXl^TEHP9^TEX2^TEMPie<‘TEX3^TEMP7^TEX4^TEHP8 

TEHP43=TEJ1P43^TEXl^TEMPie^TEX2#TEMP9^TEX3*TEHPe^TEX4»TEMP7 

DO 4628 

TEHP53=^0. 8 

TE:1P14=^. 8 

DO 4€27 1-1,1 

1F<I. ECL II>GO TO 4€17 

IF<I. ECt 1. 02. il. EQ. 1>G0 TO 4618 

111«^ 

GO TO 4619 

4617 iii=e 

GO TO 4619 

4618 Ir<l^ll. EO. 4>II1^1 
IF<I+II, EQ. 5>IIIM 

4619 IMDEK1=^IND8+583>»III 
INDEK2=IHD9^31364tl 1 1 

TEI8>52=TEJ8>53-K^FTUTU< INDEXx>«*^TU< I, 1, K, L, H, M>+CFTUTU< INDEX2>*GTU< I 

2, 2, 1C L, H, N> 

TEHP14=TEMP14^CFTUTU<INDEX2>*aTU<I, Ir K, L, K, N>^CFTUTU<1NDEX1>^TU<I 

3, 2> 1C U H, M> 

4627 CONTINUE 

PTU< 1 1, 1, KK, LL, HH, NN>=PTU< 1 1, 1, KK, LL, MM, NN>^TEMP53 
PTLK 1 1, 2, KK, LL, MM, NN>=PTU< 1 1, 2, KK, LL, MM, NN>^TEMP14 

4628 CONTINUE 

4629 CONTIIsR-lE 

4630 CONTINLff^ 

4631 CONTI^^IE 

4632 COr^TINl^ 

PEIXI, 1, KK, LL, MM, NN>-PED<1, 1, KK, LL, MM, NN>+TEMP13 
PED<±, 2, KK, LL, m, NN>=^^tXl, 2, KK, LL, MM, r>W>+TEMP23 
PED<2, 1, KK, LL, MM, NN>=PED<2, 1, KK, LL, MH, NN>+TEMP33 
PECX2, 2, KK, LL, MM, m>=Pc£D<2, 2:. KK, LL, HM, r^^>+TF.MP43 

4633 C0NTir4l.lE 

4634 C0r4TJNl-IE 

4635 CONTINiJE 

4636 CONTINUE 

Er«> PL0TE/FILM T 1 U ON PL0TE/FILM T ^ U 
El^ GRID EDGES ON GRID EDGES 

ST0RT EDGES ON INTERSECTION SOURRES RND /ICE VERSfi 

REFlD<22>CFEDIS 
DO 4684 LL=i, 4 
r>0 4683 NN^1,2 
TEMPX^ NN“16’*'LL+7. 5 
TEMP2=^:1 5-4^LL-2^r4N 
DO 4682 KK«1, 

DO 4681 MM=^1, . 

TFMP3^ i 5-4’^KK-2:nr: 

TEMP I-MM«^6*KK+7. 5 
TEMPi2 CiIS<:1, KK, 1 1.* HI!, NN> 

TEHP34^ GIS<2, KK., LL, MM, r4N> 

TEHP0^GIS<3, KK, LL, MM, NN> 

TPIS:12: 0. 0 
fPTS34M?i. O 
IPIS'V-0 M 
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TEH P 6 «TEMP1 

DO 4€80 N»l«8 

TEMP6-TEMP6-2 

TEMP4-TEMP4-2 

TEHP7«TB«>2 

T EW P8 *TE MP 3 

DO 4679 K«l>4 

TEMP7*=TE«P7*4 

T E MP 0uTE>r0 »4 

DO 4878 L«l«4 

TE MP »- T E MP 6 »X8*L 

TEMPie=TEHP4«-16>(4. 

TEHP9»8BS<TEHP9>-e. 3 
IMM.el4«INT<TEMP9> 

TEMP18=flBS<TEMP18>-8. 3 
IND2>14*INrT(TEMP18> 

DO 4877 H»l>2 
TEMP11*'TEMP7+H 
TEMP14>TEHP8«M 
TEHPll=feS<TEMPll>-*«L 7 
IND3»1NT<TEMP11> 

TEMP14=flBS<TEMP14?+e. 7 
IMD4*1HT<'TEMP14> 

INDEXl<-'INDli-IND4 

INDEX2=lNDEXli-784 

INDEX3sIND2^I»«>3 

INDEX4»INDEX3-^784 

INDEX5»INDEXl-*d588 

INDEX6sI.OEXl-»2352 

INDEX7==INDEX3+1588 

IHDEX8^INDEX3<^2352 

TEX1bCFEDIS< INDEX±> 

TEX2«CFEDIS< INDEX2> 

TEX3=K,1^DIS< INDEX3> 

TEX4«CFED1S< INDEX4> 

TEX5«CFEDIS< TNDEX5> 

TEX6=CFEDIS< INDEXO 
TEX7=CFEDIS< INDEX?) 

TEX8sCFEDIS< INDEX8> 

TGED1=:QED<1, 1, lO L, M, N> 

TGED2«GED<2, K, L, H. N> 

TGED3=GED<1, 2, lO U M, N> 

T1ED4=GED<2- 2* K, Li. H, N> 

TPIS12=TPIS12+TEXl*TGEM.+TEX2i<TGED2+TEXS*TQED3+TEX4*rr.;ED4 
TPIS34-TPlS34+TEX2*TCiEDi+TEXl*TQED<i'*-TEX4*TQED3+TEX3'*>TGED4 
TPIS5--TPIS5+TEX5*TGEDl+TEX6*TGED2+TEX7*TGED3+TEX8*TCiED4 
PED<1, 1, K, L, Mi N>=PED<li 1, K> L, M, N)+TEX1*TEW>12+TEX2>^TEMP24+TEX5*TEM 
2P5 

PED<2i 1, K> Li Mi N>»PED<2i li Ki ti Mi N>+TEX2»TEr^l2+TEXl*TEMP24+TEX6*TEI1 


3P5 

PED<li 2i Ki Li Mi N>=-PED<li 2i Ki U Mi N>+TEX3»TEMPl2+TEX4t'TEMP34+TEX7*TEM 


4P5 

PED<2i 2i Ki Li Mi N)=PED<2i 2i Ki Li Mi N)+TEX4*TEMP12+TEX3>*‘TEMP24+TEX8*TEM 
5P5 

4677 CONTINUE 

4678 CONTINUE 

4679 COt>ITm.lE 

4680 CONTINUE 

PIS<li KKi LLi MMi NN>*PIS<li KKi LLi MMi NN)+TPIS12 
PIS<2i KKi I Li MMi NN)*.PIS<2i KKi LLi MMi NN)+TPIS24 
PIS<3i KKi LLi MMi NN>«PlS<3i KK, LLi MM. NN>+TP1S5 

4681 CONTINUE 

4682 CONTINUE 

4683 CONTINUE 

4684 CONTINUE 
C 


r s — 
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EM> EDGES ON INTERSECTION SQUARES AND VICE VERSA 

START EDGES ON HAIN QRIO/FIIM T A U AND VICE VERSA 

REMIND 22 
READ<23>CFEDTU 
DO 4783 LL=1>4 
TEHP1»7. 5-16<«LL 
DO 4782 NN=d^8 
IND1*-<15#LL^NN> 

TEHP2«TEMP1-2>*NN 
DO 4781 HN=1>2 
TEMP3*=Mrt+7. 5 
TEHP4=1. 5-2i«HM 
DO 4788 KK*1^4 
TEMP3=TEMP3-16 
TEMP4=TEMP4-4 
DO 4699 
IND2=IND1+N 
TEMP5==TEMP2-2^N 
DO 4698 L=l,4 
IND2^IND2^15 
TEHP5=TEI1P5-M6 
1ND2=14^IABS< 1ND2> 

TEMP6^f®S<TEW^>+8. 7 
lND6=INT<TEr«>6>^2968 
DO 4697 H^l>2 
TEW^7=^TEHP2-i*M 
TEMP8-TErr4^H 
DO 4696 K=:1..4 
TEI*»=-7-^TEMP7+16 
TE^tf>8=^TEMP8-►4 
TEMP9=-f«$S<TEHR7>-e. 3 
I r«>9=56^ I NT < TEHP9 > 

TEHP18=ABS<TEI^8>^^0. 7 
I M>ie=- 1 NT < TEMPI© > 

I rOEXl=^ IND3 + 1 rJDl© 

INDEX2^ ir>IDEXl+742 
I r«>EX3- 1 r«>EXl+1484 
INDEX4^ lNDEXl+2226 
IND€X5==IND6-t^IND9 
INDEX6= lNDEX5^31i6 
I NDEX7=^ I NDEX5+6272 
I NDEX8^ 1 NDEXC»+94©8 
TEX1.<FEDTU< INDEX1> 

TEX2^=CFEDTU< INDEX2> 

TEX3==CFEDTU< I NDEX3 > 

TEX4^ CFEDTU *C I NDEX4 > 

TEXG^CFEDTLK INDEX5> 

TEX6-CFKDTLK INDEX6> 

TbX?=CFFDTLK IN[*EX7> 

TEX8^-C.FFDTU< INDEX8> 

DO 4690 1 1=^1. 3 

TEMP12-CjTLK 1 1. 1, KK. LL. MM. NN> 

tempi:- 0. h 

TEMPI 4=:*Gnu 11.2. KK.. L L. MM.. NN ^ 

TFMPIO^ M. M 
DO 4694 

IF<1. ECf. ; >G0 TO ’< .89 
IFvII. FO. 2>G0 10 4680 

iFai. Hc». nycio to 4686 

TFMPlIi- lEMPJ :<+TFXl =^»;a:D< I> 1. K. 1... H.. N.>-fTEXO=*GED< I. 2, K. L.. M. N> 
UMPIO- TE.MPJD+lFXr»^r,FD< 1. 1. K.. 1.. M. N>+TEX:1 I. 2. K. 1.. M. N> 

PFfi< !. 1.. K. L. M. N>*PFD< I. 1.. K- L. N. N>+TFXl>^aENF J,2+TEX5*^TFMKi4 
PFD<J.. 2. K. i... M.. N>-PFD<I/ 2. K.. I . \h N>+TFXO^ TEMPI 2^TF.Xi4^TFMP14 
(if! 10 46*=*4 
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4€&5 ig l> h. h>^vek&*ge[><u 2 ^ ig u it n> 

TEHP15bTEHP15-i>TEX 6*GED<I, U H. N>^TEX2*0EDCI. 2^ IG U IG N> 
f>ED<X^ ±*KsUH^ N>»PED<I^ IG U It N>^TEX2<t>TEHP12^TEX6*TEHP14 
P€D<X^ 2# IG U ft N>i4>ED<l^ 2> IG U It N>^TEX64^TEHP124>TCX2«>TEHP14 
GO TO 46^ 

4686 TEI6>13*TEI1P13^TEX3«i«QED<I, IG L> It N>-i^T£>C7«GED<I^ 2> IG U It N> 
TEHP15«TEHP15^TEX7#GED<I> K> It M>^TE>a^QED<l^ 2^ IG 1^ It N> 
f>£D<I/ 1# IG L> It N>«PED<1^ 1, L> It M>^TEX3^TEWP12^TEX7^TEMP14 
PED<T# 2. L/ It N>>4>ED<I# 2^ IG U It N>^TEX7^TEMP12^TEX3*TEMP14 
GO ^ ^ 4694 

4689 1F< Ea 1>G0 TO 4698 
IF<iI. EQ. 3>G0 TO 4691 

TEMP13«=TEI«>13+TEX1*GED< I, 1, K, It N>+TEX5<'GED<1. 2^ IG It N> 
TEHP15=TEHP15^TEX5>i<jED<I, 1, IG U «, N>+TEX1^GED<I, 2. IG U It N> 

F£D< 1* L, H.. N>*=PED< 1 , K, H, N>+TEX1#TEMP12^TEX5^TEHP14 
FED<I^ 2^ IG It ll>*FED<I^ 2, K> L, It N>+TEX5*TEHP12+TEX1^TEMP14 
GO TO 4694 

4690 TEHP13=TEHP13+TEX2’i<3ED< I, 1, K, L, N>^TEX6*GED< 1/ 2^ IG U It N> 
TEHP15-TEHP15+TEX6i«GED< I, 1> K, L, It N>+TEX2H5ED< 2^ IG U It N> 

PED<1> IG L, N>*PED<I, 1> IG L, It N>+TEX2<^TEMP12+TEX6*TEMP14 
PED<I# 2. IG U H, N>=PED<I> 2^ IG L> M, N>+TEX6^TEW12^TEX2*TEHP14 
GO TO 4694 

4691 TRMP12=TEHP13^TEX4>tcGED< I> 1^ K, L. It N>+TEX8^MjED< 2> IG U It N> 
TEMPlG=TEHP15+TEX8<4iED< 1, 1^ K. U N>^TEX4^GED< I, 2, U It N) 

PED< 1, 1* K. U It N>=*FED< 1. 1, K> L. It N>+TEX4^TEIi=^12+TEX8^TEMP14 
PED< 2* K> L> M, N>*PED< I, 2, L. h. N>+TEX3*TEMP12^TEX4^TEHP14 

4694 CONTINUE 

PTU< 1 1, 1^ KK, LU MH, NN>=PTU< 1 1, KK, LL, MM. »4>+TEMP13 
PTU< 1 1. 2. KK. LL. MM. NN>=PTU< 1 1. 2. KK. LL. Wt r#4>+TEMP15 

4695 cor>rrir^ 

4696 CONTINUE 

4697 COr^TINUE 

4698 CONTir^JE 

4699 CONTINUE 

4700 CONTINUE 

4701 CONTINUE 

4702 CONTINUE 

4703 CONTINUE 

END EDGES ON MRIN GRID/FILM T U PND VICE VERSft 

STRRT MAIN GRID/FILM T ^ 1 ON INTERSECTION SQURRES RND VICE VERSR 

RFNirO 23 
RFRD<24>CFTUIS 
DO 4721 LL=-1. 4 
TEMP1=^7. 5-16'^LL 
DO 4720 NN^l. 2 
TEMP2=-TFMP1+NN 
IND1^-C<^LL+NN> 

DO 4719 KK>1. 4 * 

TEMP3-7. 5“16»fKK 
DO 4718 MM^l. 2 
TFMP4-Mr-HTEMP3 
IND2^^-<3*KK+MM> 

DO 4717 N=:l/8 
TEMPr.~TEMP2-?>N 
TFHP6«TEMP4-2^N 
DO 4716 4 

IFMP?:~TFMP5+16 
TFMP6==TFMP6-M6 

TKHP7»^ RPS< rr.Hpr. >-o. i 

I ND7^1i4c I NT TEMP7 > +1 
TFMP8-RbS<TEMP6>-0. 7; 

J 1 4D8=1J I NT •: 1 FHP8 > +1 
DO 471C. rt-1.2 
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1ND3»IND1<I^ 

DO 4714 K«l>4 
1ND3«IHD2^3 
IND4«IND4-^3 
IM>5=»IABS<IND2> 

IND6^If»S<IHD4> 

IND9»IND<S-MND7 
INDie=IHD5^IND8 
DO 4713 II*1>3 
TEMP12«G1S< 1 KIC LL^ HH^ NN> 

TErtf>13=0, e 
DO 4712 1=^1^ 3 

IF<I. EQ. 3. OR. II. EQ. 3>G0 TO 4784 
INDll-616 

IF<I. EO. II>INDll3^ 

GO TO 4705 

4704 IND11»1848 

IF<I. Ea II>ir«>ll=:0 
IF<I+i:. EQ. 4>IND11=1232 

4705 ir«>EXl=^lND9^INDll 
I r«>EX2= I ND10+ 1 NDli 
TEX1=CFTUIS< INfCXl> 

TEX2=CFTUIS< INDEX2> 

TEW^13==TEI#^13+TEX1^GTU< 1, 1, K, L. N>+TEX2*GTU< I, 2> K, L. H, N> 

TLKI, 1, K. L, K. N>=^PTU<I, 1. L> M> N>+TFXl*TErff^l2 
TU< I, 2. K, L. H, N>=^PTU< I, 2, K. L. M, N>+TEX2^TEMP12 

4712 COr4TINUE 

PIS<II, KK. LL, m. m>^PlS<U. KK. LL^ mh NN>+TEI^13 

4713 CONTINUE 
471‘1 CONTINUE 
47ir. CONTINLIE 
471j CONTINUE 
471. CONTINUE 
4718 CONTINUE 
47d55i CONTINUE 

4720 CONTINUE 

4721 CONTINUE 

END HfilIN GRID/FILM T ^ U ON INTERSECTION SQURRES AND VICE VERSA 

REWIND 2‘' 

RETURN 
END 
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SUBROOTINE F>LEINF<X« Z. ^ B. f>OTL> 

TEHPZ^X^X-i-V^V-i-Z^i^Z 

1F<T»«>2. GT. e. 16E-2>G0 TO €;8e2 

IHGX^ie 

SUH=G 

lF<a LT. e. 1F-2>IMGX=^1 
1F<B. LT G. lF-2> JMftX*^! 

TFW^i=^ V^e, 55^E 
TEWP^^X^G, 55*ft 
TEW^==TEHP3-0. 

DO 6001 i^i> max 
DO 6000 J^l, JMRX 
VEMP6^TFMPl-0. 

TEHP7‘-Tf£MP6^TEMP6-HEMP5HeTEMP5^-Z^2 

R=S€«T<TEW^7> 

SUH-SUH+l^R 
6©©e CWTimE 
6001 CONTINUE 

POTL=<SUM^0. 89877E10>/<imX^JMaX> 
GO TO £mr< 

€m?: TEI^2=S0PT<TEMP2> 

POTL^O. 89877E10/TEHP2 
6002^ RETliRN 
END 


page Q 

^OOil QUAxjjJ 


A-37 



1 oooooo 


PS072SIC 


FC 


PROQRfM TO Cf«.CULATE LEAH CXPERIMEHT SENSOR INFLUENCE COEFFICIENTS 


DIHBI5ION CFm<9S7B>*CF¥ISfU<^±S^>,CFl\JFH<2l2€4>*CF^^ 

DIHEKSION CITFIH<l©976>FCFfn«^<280S>FamrH<8904>*CFB>FH<5S36> 
DIMENSION CFISFM<9408>^ CrrUTU< 18595 >. CFISIS<€:85>^CFEI>ED<7438> 
DIMENSION CFEDIS0128>, CFB>TU< 15512 ># CFTUIS<2464> 

EQUIV^ENCE <CFWM<1># CFMS|-N<1>^ CFTUFWC1>f CFEDGW<-^ >> CFPFIW<1># CFFMF 
2M<1>. CFTUFM<1>, CFEDFH<1>, CFISFM<1>. CFTUTU<1>, CFISIS<1>. CFEDED<1>^ C 

3FEDIS :X>.^EDTU<1>* CFTUIS<1>> 

DO 2005 K“l>4 
T1«=0. 0 

IF<IC HE. 1>T1«^ 

TEMPl^^e. 02921«H<-^. 051435 

IND1*=13>H:-3515 

MMl^:=l 

IF<K. EQ. 1>MM1N=^7 
DO 2094 M*=MltlN>13 
T2=^T1 

IF<M HE. 7>T2«1 
X=TEMP1^0. e^3175^M 
lND2=^INTri-H1 
DO 2083 L«l>4 
T3=^T2 

IF<L. HE. 1>T3==1 

TEMF2«=0. 02921^L-e. 95461 

IND3-690^L+IWi2 


?000 



^084 


NM:IH=:1 

IF<L. EQ. 1>W1N=8 
DO 2002 N=NM1N.15 
T4-T3 

IF<N. NE. 8>T4=^1 
V«TEMP2-»0, 093175*N 
IN03+46’»^N 
DO 1.001 I«J>2 

2-0. 0029464 h:^:I- 1> 

I NDEX== I ND4+2438HC I 
IF<T4-»'I. LT, 1. 5>Q0 TO 2900 

CRLL INFLCF<X/ V. Z> O. 8, 9. 003173. 6. 003175. 8 0. COEFF) 
CFWW< INDFX>«^COEFF 
GO TO 2001 
CFWW<INDEX> 0. 0 
CGNTUftJE 
CONTINUE 
CONTINUE 
CONTINUE 


?f003 


CONTI NUF 


WIRE ON I'ilRE WITH J=JJ COMPLETED NUMBER OF COEFFS « 4876 

ST0RT NFXT ON HIRE ON WIRE WITH J NOT EQURL 10 JJ 

DO 20:10 4 

TEMPI ^ H. 025^21 ♦1-0. 0330223 

INDL*1376+14^L 

NMJ N - i 

IF a EQ. J>NHIN. 8 
DO 20U5< N==^NMIN. 14 
X- TEMPI ♦O. 003.< 75^14 
IND2‘ INDIAN 
UO 2000 Ks«.l,4 

TFHP2^0. 02921 ♦K-O. 0339223 
INrCss-6p.6^K^ir4D2 
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IF<K. Ea 1>MM1N«8 
DO 20®7 

V*=TEHP2*0. 003175*M 
1 ND 4 * 1 N 03 :^ 49 ^H 
DO 2006 1*1/ 2 

Z»0. 0020464*<1-1> 

». e, «. eOJlTS, 0. «. 6, »317S. COEFF) 

CFHW< 1HDEX>«C0EFF 
2eee CONTINUE 

2007 CONTINUE 

2008 CONTINUE 

2009 CONTINUE 

2010 CONTI»«JE 

I«1TE<6* <CFMH< I >F I»3f 9€78f 129> 

WRITE<12> CFI*I 

end file 12 

flLt HIRE ON HIRE INFLUENCE COEFFS. CFHH< ), COttPLETEO. TOTBL NUilBER - 
<r»T HRICT OH FILH HMH SOORRE OH HIRE CFHCSHC >. 0R«> 1. FIRST. BRIO 


9€78 

2. 


DO 2015 K=l/4 

TEHP1=0- 02921*K~0. 8^38225 
IW>l“14*K-3745 


miN=i 

IF<K. EQ. 1>MM1N=8 
DO 2014 M=MM1N/ 14 
X-TEMPl+O- 003175»H 
IHD 2 =I^^Di+N 
DO 2013 L=-1r4 

TEHP2=0. 02921*L-0. 05461 
I rfl>3= I ND2+735*t 


NHlH-1 

IF<L. EO. 1>NM1N=8 
DO 2012 N»NM1N> 15 
V=TEHP2+0. 003175>*'N 
1ND4=IND3+49*N 


DO 2011 I"l' 2 

Z=0. 0126492-0. 0029464*I 
I NDEX-- 1 ND4 +2597 + 1 
CftLl. INFLC.F<Xr Vr 7., 0. 0R 0 
CFHSFW C 1 NDEX> =-COEFF 


003175.. 0. 003175r 0. 003175r COEFF) 


2011 
2012 
20.13 
20 J 4 
2015 


C 

C 

C 

C 

c 

c 


CONTINUE 
CONTINUE 
COWTINUF 
CONTI HI .IE 
CCrt’ITINUF. 

MR1TF<6.. 3100> 

WRnE'^E... 3000> Cf:FW5FW* 
MftnF<l.'i> (FHSFW 
END FIl-E 12 


>R If-2r5194r59> 


fill. 


M0IH SOURRF FUN 


ON NIFFS COMPLETED. TOTRL NUMBER = 


;trrt next on grid support STRU. 
■n.V\ F.LEHENTS.R OH HlPf^ fPlUFHv 


IUPF..TOP RNf' UNDERSIDES 

eftSF EIRS'f 


5194 

flHD EQUIVRLENT 


DO 2021 1>.l H , 

TEMR:1 C> 02921 •■♦ K— 0. 05 1 4t. f .y 
•'EMP2-J'i^K.-“ 5 
DO 2020 MM-1 r7 

.if |.'R;l-0. 00/ 1?’5+MM 
- IF’ : 
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20jL6 

2017 

2018 

2019 

2020 
2021 


2022 

2h2? 

2024 


IF<K. EQ. 1>«M8X«1 
DO 2019 H»1.HH8X 
X«TEWP3^0. 027305*^ 

TEMP5»TEHP4^M 
TEHP5«=ftBS<TEW5>^0, 7 
INDl«INT<TEHP5>-3724 
DO 2018 L»1^4 

TEMP5=0. 02921^L-0 05461 

IND2*=^1N>1^735^L 

NHIN«=1 

IF<L. EQ. l>miN«8 
DO 2017 Ns=r#IIN>15 
V*=TEHP5^0. 003175»*44 
ir«>3=INf>2-Ml9^N 
DO 2015 I«l>4 
2=0. 0 

IF<I. EQ. 2>2*0. 0029454 
IFU. EQ. 3>2=0. 0097028 
IF<I. EQ. 4>2=0. ©067564 
INDEX« IND34*2597^I 

CftLL II^LCF<X. V. Z. 0. 0. 0. 003175> 0. 003175.. 0, 001905/ COEFF> 

CFTUFW< INDEX>=COEFF 

cormr^ 

cmumE 

CONTINUE 

CONTINUE 

CONTINUE 

CONTir4UE 

GRID SUPPORT STRUCTURE/ TOP ftND UNDERSIDES/ PND EQUIVBLENT FILM ELEMENTS 0f*4 
MIRES WITH J=JJ/COMPLE ED. CFTUFM< >. TOTAL NUMBER =10388 

START NEXT ON CFTUFW< >FOR J NOT EQUAL TO JJ 


DO 2027 L«l/4 

TEMP1=0. 02921HCL-0. 0530225 

ir>»Dl=14x‘L+7598 

NMIN=i 

IF<L. EQ. 1>NH1N^8 
DO 2026 N=NHIN/14 
X=TEMPl+0. 003175’^N 
IND2^ INDl-i-N 


DO 2025 K> l/4 

TFMP2=.0. 02921^ K--0. 05588 
TEMP3=16’^K-'8. 5 
MMAX-2 

IF^K. FQ. 1>MMAX=1 
DO 2024 M=-i/MMAX 

TFMP4=TFWP2+0. 0273O5:4^M 
TEMP5^TEMP3+M 
DO 2023 N14-1, 8 

V=^TFMP4-0. 003175^NN 
TE MP6- TF MP5-2^ NN 
TFMPt>- 49’*' ABS • “"F HP6 > +0. 7 
IND3=INr>24‘IN'! <TEHP6 ) 

DO 2022 1=1/4 
Z=0. 0 

IF<T. FO. 2>7»0 0029464 
IF < T. F»;^. 3 >7' 0 009*; 028 
JF<I. FQ. 4 *7=0. OHF./064 
INI'FX^-jr4M -^2V44t I 
CAU 1 44F I ,r.F •: X/ V/ 7/ ©, 0.. ©. 
f;F1UFW< TNDFX>’*r i ih'FF 
cowr:fiUF 
comiNJ »F 
(0NT1NI.IP 



:<175/ 0. 00X905/ 0. 0CCJ75/ C.orKF V 
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2025 CONTINUE 

2026 CONTINUE 

2027 CONTI»^ 

M%ITE<6, 3100> 

l«%ITE<€> 3000><CFTIFM<I>, I*=4, 21364. ^9) 
NRITE<12> CFTUFW 
EW> FILE 12 


0LL CFTUFH< ) COMPLETED. N»>®ER IN THIS BLOCK » 10976. TOTftL ■ 21364 


START NEXT ON PLATE EDGES ON MIRES WITH 3«JJ. CFEDQH< > 

« 


2028 


2029 


2020 


2021 


2032 


202S 


DO 2033 

02921^K-0. 05401 
TEI1P2=14#K-7. 5 
DO 2032 

TEHP3-TEMP1-0. e«3175>^Hri 

TEMP4«T0tP2-2*HWl 

HM0Xe2 

IF<K. EQ. l>MHfiX«l 
DO 2031 M=i.MM0X 
X«TBiP3«i-0. 0254^H 
TEMP5=TEHP4+H 
TEI^5=^0BS<TEMF5>-*‘0. 7 
I MD1« I r 4T < TEHP5 > -3724 
DO 2030 L-l>4 

TEMP5-0. 02921‘»^L-0. 05461 
IND2=IND1+7554cL 

IF<L. EQ. 1>NHIN=^8 
DO 2029 H«HHIN>15 
V^TEMP5+0. 003175>**N 
IND3=^U4D2+49>*'N 
DO 2020 1=^1. 2 

2*=0 »029464^<I-1> 

I NDEX^ I ND3-^2597^ I 

C0LL INFLCF<X. V, 2> 0. 0. 0. 003175. 0. 083 175. 0. 0> CO£FF> 
CFEDGW < I NDEX > =^COEFF 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 




PLATE EDGES ON WIRES WITH J=^JJ COMPLETED. NUMBER OF COEFFS « 5194 

START NEX’- ' .4 PLATE EDGES ON WIRES WITH J NOT EQUAL TO JJ. CFEDGW<5195 0N> 


DO 20:s9 L==1.4 

TEMPING. 02921 )fL-0. 0530225 

TEMP;>-14>»L-'21. 5 

NHIN-l 

IF<L. FQ. d >NHIN--A 
DO 2039 N^^NHIN. :14 
X« TEMP 1+0. A03l?5%N 
TFMP3- TEHP;:+N 
TFMP3- APS* IFMP < >+2426. 2 
DO 2037 

TEHPF. 0. O292i.4K.-0. «i3W;‘S 

DO 20: NI4' ,1.8 

TFHK I I MP8-0. 00 Vi 7V.+NN 
TEMH>:.. ThHP'l-^+NN 
UMAX- 2 

1F<K EO. il >MMAX^;t 


Qf POOE 
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2024 

2025 

2026 
202> 
20^ 
2029 


^C\4? 

pcuO 


DO 2025 

V»TEMP5^0. 0254HCM 
TEW>7«TEHP6-^H 
TEI*^7«49«*4^S<TEHP7>-MENP2 
ir«>l^INT<TEHP7> t 

DO 2034 l«l/ 2 

2=0. 0029464**‘<I-1) 

1 r«>EX= I NDi+2744*** I 

CmJL INFLCF<X> V, 2. 0. 0/ 0. 003175. 8. 0. 0. 803175/ COEFF> 
CFEDQW< INDGX)=COEFF 
CONTINttE 
CONTINUE 
CWTINUE 
CWTINUE 
CONTINl^ 

CONTIhH.lE 
i4RITE<6/ 3i00> 

W^ITE<6. 3000><CFEDGW<I>. 1=2. 10682/ 89> 

WRITE<12> CFEXXM 
END FILE i;. 

RLL PL0TE EDGES ON WIRES COMPLETED. TOT«L OF CFEDa4< > IS 10682 
ST8RT NEXT ON PLfiTE/FlLM INTERSECTION SQUHRES ON WM8ES. CFFF1M< > 


DO 2046 K=l/4 

TEMPi=e. 02921*K-0. 0574675 
TEI-r2=14’*'K-7. 5 

IF<K. EQ. i >MM8X=l 
DO 2045 M^i^MMRX 

TEI1P3«=TEMPl-^0. €27305*M 
TEI1P4=TEI*lP2-*-M 
DO 2044 MM=l/7 

X»TEMP3-0. 00317^5♦^U1 
TEMPr>«^TEMP4- 2^MM 
TEMP5‘--ftBS< TFMP5> 

DO 2043 L-1/4 

TEMP6^^0- 02921*^L-0. 05588 
TEMP?^=16>kI ~8. 5 
NMft>:=^2 

IF<L. EQ. i >Nfl8X^l 
DO 2042 Nt=:l, NMFiX 

TFMP8=-TF.MP6+0. 0?7305^N 
TEMP9:^TFMP7+r4 
DO 2041 NM=^1>8 

V= TEMPS-0. 00?:1 75H NN 
TEMPI 0^ 1 FMP^-2^N?I 
1 FMPJLO^a?*^ BPS < TEMP .1 0 > -*-TF:MP5-^0. 2 
I NDl--^ I NT < TFMPit 8 > -2768 
DO 20*10 I--=l>4 
0 

IF<I. EQ. 2>Z-0. 0029-*t.M 
IF<T. EQ. 3>Z=0. 0097028 
IF<I. F.a 4>Z= 0. 0067564 
I NDEX« I NDJ.+2744 4^ I 

CBLL INPLCF'.;X/ V. 7, 0. 0. 0 003'i75/ 0. 001905. 0. 001905. r,OPFF> 
L.FPFIWV ) r^DFX>«^COEFF 
CONTINUE 
CONTINUF 
ruin INUF 
f .OMT 1 N* *F 

r,.,UT JNUr 

O'UI INMF 
‘ ..m INIIh 
Kp • n* '.6‘ .» 
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:< 00 e><cFPFiw<i>/ i»ii. 10976 . 129 ) 

CFPFIW 
END FILE 12 
C 

C RLL PLRTE/FILM INTERSECT I Or4S ON WIRES COMPLETED. TOTAL OF CFPF1W< >^10976 
C 

C ALL CONTRIBUTIONS TO WIRES HAVE BEEN COMPLETED. NUMBER OF COEFFS .894 
C 

C START NEXT Oi4 CONTRIBUTIONS TO FILM MAIN SQARE ELEMENTS. NOTE THAT CONTRIBU 

C TIDNS FROM WIRES ARv£ THE SAME AS FILM ELEMENTS TO WIRES 

C 

C START OF FILM MAIN SQUARE ON FILM MAIN SQUARE ‘ CFFMFM* > 

C 

DO 2051 K-1.4 
T1:==0. 0 

IF<K. NE. 1>T3==1 
TEMPl-0. 0292^'^K-0. 054 Si 
IND1=15»»*K“1241 
MMIN=-1 

!F<K. EQ. 1>MMIN=8 
t>0 M=MMIN. 15 

T2^T1 

IF<M. NE. 8>T2=i 

x==ti:mpi+c. 00:<i75hvM 

IND2-=IND1+M 
DO 2049 L«i. 4 
*; 2=T2 

IF<!.. NE. ±>JZ^± 

TEMP2«:0. 0292i^L.--0. 05461 

INDri;-IND24.795^L 

NHIWa^l 

IF<L. EQ. 1>NMIN==8 
DO 2048 N“-=K*'1IN. 15 
T4^T3 

IF<N- NE. 8>T4=1 
V-TEMP2+0. 003175 ♦^N 
IWr>EX“.J|»D3+53*N 
IF<T4. LT. 0. 5 >»j0 TO 2047 
TEMP3=^0. 001175 

CALL I WFLCF<X. V. 0. Cl, TEMPI, TEMPI, TEMP3, 1 FMP3, COEFF> 

CFFMFM V I WDFX > ^COEFF 
GO TO 2048 

2047 CFFMFIK J Nr.'FX>«0. 0 

2048 COr*TlNUF 

2049 CONTINUE 

2050 CONTINUF 

2051 CONTINUE 
WRITEC6, 110»n> 

WRITEC6, 1000* ac^FMFMv ! i=:j , 2809, I<9> 

WRITEi.:L2> ( FFM^M 
END FILE 12 
C 

C FII..M MAIN SQUARE ON FIl M MAIN SQUARE COM TED. T^TAL CFFHFH'! > ^ 2809 
C 

C START NEXT ON GRID SUPPORT STRUr.TljR{r. ^Qp UNDERSIDES. AND FWI'-^ALENT 

C FILM EUHENTS 0 4 FILM MAIN S»XIARE FLEr S. CFTUFM< ^ 

C 

DO 205 r . =^1, 4 

TE;1P:1 -0. 02921%K-0. 05588 
temp;.* : •A64-K-8. 5 
MMAX^ 2 

1F<K EQ. 1>MMAX-1 
DO C056 M^^1,MHAX 

TF.HFl^- TFMP r^-0. 027105>t'M 
TFHP4-TFHP2+M 
DO 205t. 
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2052 

2053 

2054 

2055 

2056 
2^7 


7055* 


^<=TEMP3-©. 00317S«»tt 
TEHP5=TEHP4-2*MH 
TEHP5=flBS<TEMPS)'*«. 7 
INM:=INT<TEnP5>-42S£ 

DO 2054 L=l,4 

TEMP5«=e. 02921>*t.-0. 054€1 

IM>2»IM>1-»^840>*L 

NHIN=1 

1F<L. EQ. 1>NHIN=^ 

DO 2053 N=4t1IH, 15 
<^TEMP5^0. *.^a175^ 

. 'D3«IND2- 5 j44 
DO 2052 1=1,2 

i=0. 0 

IF<1. EQ. 2)Z=0. 0007028 
1F<I. ECL 2>2=0. 0067564 
INDEX=IMD2+2068»I 
TEMP6=0. 003175 

CflLL 1NFLCF<X, V, Z, TEMP6, TEMP6, TEHP6, 0. 001905, CCCFF> 

CFTUFH< INDEX)*eOEFF 
CtfNTINUE 
CONTINUE 
CONTINUE 
CONTirflJE 
• NTINUE 
CONTIM.IE 
HRITE<6, 3100> 

MRlTE<6,3000><CF'nFH<I>, 1=4, 8904, 89 > 

HRITE<12> CFTLIFM 
END FILE 12 

GRID SLE>PORT STRUCTURE, TOP RI4D UNDERSIDES fiND EQUIVRLENT FILM ELEMENTS ON 
FILM MAIN SQUARE ELEMENTS COMPLETED. TOTAL CFTUFM< > = 8904 

START I«XT a<:ID ET^GES ON FILM MAIN SQUARE ELEMENTS. CFEDFM< ) 


DO 2063 K=l,4 

TEHPi=0. 02921»K-0. 0530225 

TEMP2=16*K-8. 5 

MMAX=2 


IF<K. EQ. 1>MMA.X=1 
DO 2062 M=l, MMAX 

TEMP3=TFHPl-*-0. 0254»H 
TEMP4=TEHP2+M 
DO 2061 MM=1, 0 

X=TEMP3-0. 003175»MM 
TEMP5=TEMP4-2kMM 
TEMP5=ABS< TEMPO >+0. 7 
I ND1= I NT < TEMPS >-4256 
DO 2060 L=l, 4 


TEMP5«0. 0.'?921>^L-0. 05461 

IND2-IND1+S40H- 

NMIN=“1 


IF<L. EQ. 1>NMIN=8 
DO 2059 N=NHIN, 15 
V=lFMP5+0. 003175fN 
IMD3=-IND.^+56*N 
DO 2058 I =1,2 
Z-0. 00970.20 


oeigwaljageb 

OF. POOR quality 


IF<I. \ Q. 2>7--0. 0067564 
I NDr X=- 1 nrv.+2968 *■ 1 
TEMP6 0. 00 A 175 

CALI. U4H.r.K.;X., V,. 7. ieHP6, TEMP6, 1FMP6, tv A, COF.FP ) 
CFhDiMK ir^H :;>• i mkfk 
tom iNui- 
CONTINUE 
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2060 

2061 

2062 

2061 


C 

C 

C 

C 

C 


206<1 

206r» 

2066 

2067 

2068 

2069 

2070 


C 

C 

C 

c 

c 
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CONTIHUE 

COHTII4UE 

COHTINP^ 

CONTINUE 
WRITE C6# 1100^ 

HRITE<6,ie00><CFEDFH<I># I«ll, 5936^-79> 

WRITE<12> CFEDFM 
END FILE 12 

GRID EDGES ON FILH HAIN SQURRE ELEHENTS COKPLETED. TOTRL CFEDFH< >«5936 

STRR NEXT INTERSECTION SQURRES ON FILH HRIN SQURRE ELEMENTS. CF1SFH< > 

DO 2070 K=1^4 

TEMPl^. 02921^K-0. 05588 
TEMP2=16<4<-8. 5 
MMRK^2 

IF<IC EQ. l>^ttiRX=l 
DO 2069 H=1^MI-«X 

TEMPl=TEI^l+0. 0272fet5^ 

TEW>4=:TENP2-^M 
DO 2068 MM=^l.-8 

X=TEMPl--0. 001i75^MM 
TEMP5=^TEHP4-2»*'MM 
TEMP5=^RBS<TEMP5>-*^. 7 
1W>1=^INT<TEMP5> 

DO 2067 L=^l>4 

TEMPS - 0. 02921*L^0. 05588 
TEMP6=16^l.-:-: 5 
NMRX=^2 

IF<L. EQ. l>r^fiX==l 
DO 2066 N=1,NMRX 

TENP7=TEMP5+0. 027305^N 
TEMP8=TEMP6+N 
DO 2065 NN=--1>8 

V=TEMP7-0. 002175’»NN 
TEt1P9-TEMP8“2^NW 
TFMP9=^RBS<TEMP9>-«0. 1 
I ND2^ I NDi+56>»: I NT < TEMP9 > 

DO 2064 1=1.. 1 
Z- 0. 0 

IF<I. EO. 2OZ=0. 0097028 
IF<I. FQ l>Z=-0. 0067564 
1 NDFX= I •: I -1 > 

TEMP9^0. 0O:<175 
TFMPJ0=0. 00:1905 

CRLL INFLCF^X. V.. Z.. TEHP9.. TFMC>9. TEMPiO. TEMPI©. COFFF> 

CF I SFM< I NDK:> ^COFFK 
CONTINUE 
lONTINUF 
CONTINUE 
CONTINUE 
CONTINUF 
CONTINUE 
CONTINUE 
WPITF'>“. .I:t00> 

WklTF* 6. :i000*^M:;FIS*^M<r>.. 9408. 47> 

WRiTF':.l2> ChJSFM 
END FILE 12 

INfEPSECTION SOUftRES »»N FILM MR IN SOUftRES COMPLETFD. TOTRL CFISFM< ,*=9408 

STRPT NEXT Or4 PLfiTE/FJLM T BND U ON PLRTF/FIl.M T RND U. CFTUTUC >FOR J=JJ 

r*0 2076 K-i. 4 
Ti^ 0. 0 


A-45 



r> o o o o o 


BSR 4234 


2071 

2072 
2073: 

2074 

2075 

2076 


IF<IC ME. 1>T1»1 
TEMPl^. e2921»K-e. 08382 
lMWL=3^K-840 

IF<K. Ea 1>HMIM»2 
DO 2075 n=^41N>3 
T2«T1 

1F<K ME. 2>T2»1 
X«TEMP1-*^. 027305*+l 
IND2^1ND1^ 

DO 2074 L»l,4 
T3=T2 

1F<L. ME. 1>T3*1 
TEHP2=0. 02921*itL-0. 05461 
IND3=lND2<i>l€5>i4. 

IFCL. Ea 1>I«IN=8 
DO 2073 N:=NMIM,15 
T4=T3 

IF<N. NE. 8>T4»1 
V=TEHP2+0. e03175*N 
IMD4=IND34.11.|cN 
DO 2072 1=1.5 
Z=0. 6 

IF<I. EQ. 2>Z=0. 0001524 
IF<I. EQ. 3>Z=e. 0029464 
IF<I. EQ. 4>Z=0. 0097020 
IF<I. EQ. 5>Z=0. 0067564 
INDEX=IND4^583*I 
IF<T4+I. LT. 1. 5>G0 TO 2071 
TEMP3:^0. 001905 
TEMP4=0. 003175 

CRLL INFLCF<X. V. 2. TEHF3. TEMP4. TEMP4. TEMP3. COEFF> 

CFTUTU a NDFX > =COEFF 
GO TO 2072 
CFTUTUaNDEX>=0. 0 
COTiTirWE 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 

01J. J=JJ PLRTE/FILM T RND U ON PLRTE/FJLM T RND U COMPLETED. NUMBER = 2915 


STRRT NEXT ON SRME FOR J NOT EQURL TO JJ. CFTUTU< > 


DO 2083 L . 4 

TEMPI m3. 0292J^L-0. 00254 
TEMP2=.16’»L“8. 5 
MHMRX=2 

IF<L. EQ. 1>MMNRX=J. 

DO 2082 MM^i.. MMMRK 

TFMP:<^ lKMPi-0. 0?7305*=MM 
TF:MF4^ TEMP2+MM 
DO 2081 U‘ ‘3.. 8 

X=TEMP 003 j 75^ f4 

TFMP5~TFNP4•-?'^N 
TFMP5«R8?;<TEMP5>+0. 7 
ir4Dt.- J NT ». TEMPO .'-221 
DO 2080 K-M .. 4 

7FliPr,^v?i. O;::92i^K-0. 65588 
TKHP6‘-’16^K-8. 5 
MI‘IRX 2 
IF<K. EQ 

DO 207^ M.^1.*MMRX 
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2077 

2078 

2079 

2080 
2081 
2082 
2082 


TEHP8»TEHP6^M 
DO 2078 NN^8 
V»TEMP?-0, 003: 

TEnP9»TEHP8-2#NN 
TEHP9«A8S<TEMP9>-0. 3 ► 

IND2»INDl*i'58*INT<TEMP9> 

DO 207? 1*^1^ 5 

Z*=0. 0 

1F<I. EO 2>Z:=0. 0001524 
IF<I. ECL 3>Z=0. 00294€4 
IF<I. Ea 4>Z=0. 0097028 
IF<I, E€l 5>Z=0. 0067564 
1NDEX-IND2^3136%I 
TEHP9=e. 001905 
TEHP10*0. 003175 

CftLL IHFLCFOi, V, Z. TEMP9, TE«P10> TEMPS, TEHP18, COEFF> 

CFTIJm< INDEX>=CC»=F 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 
MRITE<6, 3100> 

MR1TE<6, 3000><CFTUTU<I>, 1=3, 18595, 83> 

HRITE<12> CFTUTU 
END FILE 12 

«.L PL8TE/F1LM T 8ND U ON PLfiTE/FILM T RND U COMPLETED. FOR J NOT EQU8L TO 
JJ NUMBER IS 15680. WITH 2915 FOR J=JJ T0T8L CFTUTU< > =18595 

ST8RT NEXT ON INTERSECTJW SQUfiRES ON INTERSECTION SQUARES. CFISISC > 

DO 2089 K=l, 4 
T1=0. 0 

IF<K. NE. 1>T1=1 

TEI-IP1=0. 02924 ♦K“0. 08382 

IND1=3%K-180 

MM1N=1 

IF<K. EQ. 1>MMIN=2 
DO 2088 i1=MMIN.. 3 
T2=T1 

1F<M. NF. 2>T2=1 
X=TEMPl+0. 027305>*^M 
IND2^INr*i+M 
DO 2087 L=l .. 4 
T3=^T2 

IF<L. NE. 1>T3=1 
THriP2=0. 02921 >*^L-0. ©8382 
IND3- JND2^33>^L* 

NMIN" 1 

IF^T.. FQ. 1>NHIN:^2 
DO 2086 r4=-NI1irL3 
14^13 

IP<r4. NF. 2>T4^:I 
V=^TFMP2-^0. Ci.‘ V7 01V+ N 
INIH JI4r»3^.J I *44 
DO 2085 J~J .. 5 

a 0 

IF<I. EQ. 2 >7=0. ‘0001524 
IF<.J FO. ?> 7=0. 0029464 
IF a. FQ. 4 >7=0. 0097028 
IF<I. FQ. 5 >7=0. 006756^1 
INDFX=IND4+J2i:+:I 
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2085 

2086 
208? 
2088 
2089 


C 

C 

c 

c 

c 


2890 

2092 

209> 

209-1 

209r 

C 

c. 

c 

c 

c. 
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1F<T4^I. LT, 1. 5>G0 TO 2884 
TEHP3=e. 001905 

CftLL 1NFLCF<X> V, 2. TEMPI, TEMP3, TEMP3, 1 *~MP3, COEFF> 

CFISIS< INt€X>=COEFF 
GO TO 2085 
CFISIS<ir«>EX>=0. 0 
CONTINUE 
COr>ITINUE 
CONTINUE 
CONTINl^ 

COI4T INLIE 
MRITE<6, 3106> 

MRITE<6, 2000><C4=^IS1S<1>, I=5>605, 12> 

MRITE<12> CFISIS 
END FILE 12 

INTERSECTION SOURCES ON INTERSECT I ON SQUARES COMPLETED. TOTRL CFISIS< >*685 
STRRT r^XT ON GRID EDGES ON GRID EDGES. CFEDED< >. J=JJ FIRST 


DO 2095 K=^l,4 
T1=0. 0 

IF<K- NE. 1>T1-1 
TEMPING. 02921’«'K-e. 08801 
II«>1=33^K--840 
MMIN=^1 

IF<K. EQ. 1>MMIN=2 
DO 2094 M=:MMIN, 3 
T2=T1 

IF<M- NE. 2)T2=1 
X«TEWl+0. 0254>^M 
ir02=^INDl-»*M 
DO 2093 L^l.*4 
T3^T2 

IF<L. NE. 1 >T3^1 
TEMP2=^0. 0292J^L-0 05461 
ir«>2==INr>2+165^L 
NM1N=:1 

IF<L. EQ. 1>NMIN=8 
DO 2092 N=NMIN» 15 
T4=T3 

IF<N. NE 8>T4^ 1 
V-TFMP2+0. 0031 75^N 
IND4“ I t4D3+J 
DO 2€i91 1=1,2 
2=0. 0 

IF<I. EO. 2>7=0. 0029464 
INDEX=^IND4+5S3*I 


IF<T4+I. LT. 1. 5>ri0 TO 2890 

CRI L JNFLCFO:, V, Z, 0 O, 0. 003:175, O. 003175, ©. 0, COFFF> 
CFEDED < 1 H\ ‘EX > =rOFFF 
GO TO 209.1 . 


C.FFDFrM:iNi.'FX>=^0. 0 
CCfNTJNME 
rONTINl IF 
CONTINUE 
CONI iraiE 
CONIINUF 



GRID Fr»GFS (*N GRIP' FD0FS, vFH 'ErM; >,wnH J -.U CCiHPLElED. r4UHBt>- 1166 

STHRT NbX1 0?4 GRID EDGES 044 UklfA EDGFS WITH J N01 EOURL TO JJ 

DO 3002 L=:l,4 

TFHP.UA. 02:92 J ^1-0. O053975 
TFHP2=^l6>»l.-ft 5 
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2096 

2097 

2098 

2099 

2100 
3001 
2002 


t 

C. 


HHHRX-2 
IF<L- EQ. 

DO 3001 

TEHP3=^TEI^l-0. 0254^HH 
TEHP4^*TEMF2-ittt1 
DO 2100 N^l>8 

X=TEMP3^0. 003175^N 
TEHP5=TEMP4-2<^ 

TEHP5=ftBS<TEI^5>-^0. 7 
INDl=INT<TEHP5>-197e 
DO 2099 K=l,4 

TEHF5=a 02921>M<-0. 0530225 

TEHF>6=16^K-8. 5 

HMRX=2 

1F<K. EQ. l>MiftX«l 
DO 2098 H=1.M*-WX 

TEr^7=TEHF5+e. 0254^M 
TEI^=TEMf>6+H 
DO 2097 NN=1.8 

V=TEI^7-0. 002175>t^Nr4 
TEHF9=TEMP8-2^NN 
TEMP9^RBS<TEMP9>-0. 3 
I ND2^ I NDl+56^ j NT < TEMP9 > 

DO 2096 1=^1^ 2 
Z=0. 0 

IF<I. EQ. 2>Z=0. 0029464 
I NDEK= 1 ND2+3136’*' I 

CftLL U4FLCF<X> V, Z.* 0. 0. 0. 083175. 0. 0. 6. 063175. COEFF> 
CFEDED< INDEX>=COEFF 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 
COI4TINUE 
14R!ITF<6. r5l00> 

HR1TFC6. 2000>*;CFEDED<I >. 1^ 2. 7438. 52> 

NRITE<:t2> CFEUFO 
EWP FILE 12 

OLL EDGES 0^4 EDGES COMPLETED. TOTRL = 1166+6272=^7438 


STRRT NFXT Or-4 Ff^GES ON IWTERSFXTIOr4 SQUARES. CFEDISC > 

DO 3009 K=^:l. -^1 

TEMPI =0. 0292*J+f>0. 0263525 
TEMP2==‘4^=K-2. 5 
MHA>;- 2 

I Fa.: FO J >MMhX=- 1 
DO 2E»05r: H=- I ■ HMAX 

temp:* - TEMP \ +0. ri2r.4+H 
TFMK 4 - TH Ih\-+H 
D(» "0»tV MH'J»2 

X aPMP •'0. 027. Gr.tMM 
TEMP' .•-UHP4“2'»MH 
TEMP^.^-hPSMPHK.. »+m 7 
ir4r*j-.£i4r oFMPo:*-7P4 
DO :0?i06 1. 1.4 

TFHPO O. 0292J +L-0. 00254 
TFMPGt J 6+1 -8. 5 
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2002 

2004 

2005 

2006 

2007 

2008 
2009 


C 

C 

C 

C 

C 

C 

C 


IF<L. EQ. l>NNMftX»l 
DO 2005 NN=l.mMRX 

TEMP7*TEMP5-e. 02720544W 
TEHP0=TEI1P6^NN 
DO 2004 N-1.8 

VS5TEW7+0. 002175 N 
TEMP9«TEW>8-2^ 

TEHP95^0BS<TE^tf>9^>-e. 2 
IND2=INDl+14^INT<TEr«>9> 

DO 2002 1=1.4 
Z»0. 0006762 
IF<I. FQ. 2>2»0. 0029464 
IF<1. EQ, 2>2=0. 0097028 
IF<I, EQ. 4>2=0. 0067564 
INDEX= ir«>2+784^ I 
TEHP9=0. 001905 

CnJL IHFLCFOi, V. Z. TE^tf>9. TEHP9. 0. 002175. 0. 0. COEP"> 
CFED I S< I HbEX > :^COEFF 
CONTINUE 
CONTINUE 
CONTINUE 
COHTINLIE 
COr^lTINUE 
COI^T INLIE 
CONTINUE 
WRITE<6. 2100> 

WRITE<6. 2000><CFEC>IS<I>. 1=4. 2126. 87> 

WRITE<12> CFFDIS 
END FILE 12 

EDGES ON INTERSECTION SQURRES COMPLETED. TOTRL CFEDIS< > =2126 


STftRT NEXT ON EDGES ON SUPPORT STRUCTURE/FILM T 4 U CFEDTU< > 
FOR J=-JJ 

DO 2015 K=^l. 4 

TEWPl-0. 02921^K-0. 0262525 
TEMP2=4^K-2. 5 
MMRX=-? 

IFCK. FQ. 1>M!1AX=J 
DO 20i4 M=l.HMfiX 

TFHP2-TFMP:1-»-0. 0254>^M 
TFMP4‘-1EMP2+M 
DO 2CvJ:< MH-1.2 

X=TF.MP2-0. 027205^MM 
TFHP5=TEHP4-2^MM 
TEMP5=:^ttBS<TFMP5>+0. 7 
I NDri - 1 NT •: TFI*IP5 > - A 064 
DO 7S0A2 L-i.. 4 

TEHP*s>-0. 0292:1 '♦A -0. 05461 
IND2- INDA4*2A0^L 
NMJN-1 

IF<L. EQ. 1>NMIN=8 
DO 2011 N=NMIN.. J5 
V=TFMP5+0. 002i754<N 
IND2^IND2+14=+=N 
DO 2010 1=1.. 4 
Z=0. 0000762 
IFCI. FO. 2 >7=0. 0029464 
IFU. FCs!. 2 >7=0. 0097028 
IF<1 FO. 4>7=a 0067564 
ira>EX- iNro+74:4i 
TKHP6^0. 002:1 : 5 
CPU.. JNFIXL vX.. V. 7.. 0. 001905. TFHP6, 1FMP6, 0. 0. r.riEFi*> 
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3810 

3811 
3012 

3813 

3814 
3015 

C 

C 

C 

C 

C 


3016 
30 J 7 
3018 
3015# 
30;>0 
302:1 
3022 


C 

C 

C 

C 

r 

( 


CF€DTU< 1HDEK>*^C0I£FF 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 

EDGES ON SUPPORT STRUCTURE/FILH T 1 U FOR J«JJ COMPLETED. NUMBER » 2968 

STRRT NEXT ON SRME FOR J NOT EOURL TO JJ. CFEDTU<2969 ON > 

DO 3822 L=^l>4 

TEW>1«0. 02921^L-8. 88254 
TEMP2»16*N.-8. 5 
MMI-IRX=2 

IF<L. EQ. l>W«1flX=l 
DO 3021 m-l^NMMRX 

TEMP3=^TEMPl-0. 827305*MH 
TEMP4=TEMP2-H4H 
DO 3020 N=^1>0 

X=TE^F3♦a 0031753|»N 
TEMP5=TEriF4-2’«‘N 
TEMP5=BBS<TEMP5>+0. 7 
INDd - INT<TFMP5>-168 
DO 3019 K=1.4 

TEMP5=0. 02921>^K“0. 0529225 

TENP6---16HCK-8. 5 

MMBX=2 

IF<K. EQ. 1>MM«X=1 
DO 3018 M=l.MMfiX 

TEI1P7==TEMP5-i*e. 025^’»^M 
TEHPe^^TEMPE-^M 
DO 3017 m^±.s 

V=lEMP7-0. 003175^NN 
TEMP9=--TEMP8-2^NN 
TEMP9=^RBS<TEMP9>-0. 3 
I ND2= I NDl+56^ I NT < TEMPS > 

DO 3016 1^1,4 
Z=0. 0000762 
IF<I. EQ. 2>2=0 0029464 
IF<I. EQ. 3>2^0 0097028 
IF<I. EQ. 4>7=0. 0067564 
I NDFX«= I ND2+3136^ I 
TEMP9=0. 003175 

CBI.L INFLCF<X> V. Z. 0. 901905/ TEMP9/ 0. 0/ TEMP9/ COEFF) 

CFEDTIK INDEX>=^COEFF 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 
WRITE<6/ 3100> 

WRITF.<6/ 3000><CFEDTU<I)/ I«=2/ 15512/ 235> 

WRITE Cl 2> CFFDTU 
END FILE 12 

ALL EDGES ON SUPPORT STRUCT!. IRE/FILM T AND U COMPLETED. TOTAL CFEDTUO-15512 

START NEXT ON SUPPORT STRUCTURE /FILM T AND U ON INTERSECTION SQUARES. 
CFTUISC > 

DO 302ft K= l/-1 

TFHP'Uft. 0292 !^K-0. 98382 
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INDl*=3*M<-620 

NHIN»i 

IF<IC ECL i>HHIN«2 
DO 3027 HNiHlN>3 
X«TEMTi-»«. 027205 *kH 

DO 3026 L=l>4 

TEMP2*^0. 02921*L-0. 08254 
TEMP3»^16*L~8. 5 
N^«4RX>^2 

1F<L. EQ. l>NWIfiX«l 
DO '3025 NN^l>NNmx 
TEMF4*TEr^2--^ 

TEHf^^TE^IP^ m 
DO 3024 N .08 

V-TEfF' 002175*N 
TEWH&= ^E^tf>5-2♦N 
TEMP6^ ftBS<TEI^6>-0. 3 
If«>3= ND2^-11^INT<TEMP6) 

DO 30 >3 I“l,4 
Z=0 0 

IF<1 EQ. 2 >Z=a 0029484 
1F<I. 3>Z*=0. 0097028 

IF<I. EQ. 4>Z*0. 0067564 
INDEK=IND3+616*I 
TEI^7=0. 001905 

CRLL INFLCF<X> V/ Z, TEHP7/ TEf*E>7/ 0. 003175> TEHP7, COEFF> 
CFTUIS< INDEX>=COEFF 

3023 COI^TINUE 

3024 CONTINUE 

3025 CONTINUE 

3026 COt^TINUE 

3027 CONTINUE 

3028 CONTINUE 
WRITE<6. 3100> 

WRITE<6. 3000><CFTUIS<I>, 1=4, 2464, 41> 

WRITE<12> CFTUIS 
END FILF 12 

RL SUPPORT STRUCTURE/FILH T mO U ON INTERSECTION SQURRES. COMPLETED. 
TOTRL CFTUIS< >=2464 

RU- INFLUENCE COEFFICIENTS CRLCUI.RTED. TOTRL NUMBER=132701 
3000 FORMRT < IK 2X, E10. 3 > > 

3400 FORHRT< '0'% -^NEW DRTR SETV> 

END 


ORIGINAL PAGE IS 
0? PCX)B QUALITY 
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SUBROUTINE 1NFLCF<X. Z/ B> S# T. COEFF> 
TEM>2«Z^Z 

TEHP1*X»X^V»V-MEMP2 

1F<TEMP1. 6T. 0. 16E-3>G0 TO 1004 

IMAXt^ie 

JMRX-10 

KMRX^ie 

SUH«0. 0 

iF<a LT. a iE-2>irw<«i 
IF<B. LT. 0. lE-2> JMftX«l 
IF<T. LT. 0. lE-2>KfW<«l 
TEW>8=T 
TEHF9=0. 0 

IF<S. GT. 0. 1E-2>Q0 TO 1000 
TEW^8=0. 0 . 

TEHP^aT 

KMftX«l 

1000 TEHPl=V-^0. 55^B-0. 

TEMP2«X^0. 55>fR-0. 55*TEW»8^-0. Sh^TEHFB 
DO 1003 K=1.KI1AX 

TEHP4=TEMP2+0. l^i^K^TEHPS 
DO 1002 I«l. IMRX 

TEHP5=<TEMP4-0. 1^R^I> 

DO 1001 J=^l. JMRX 

TEMP6==TEriPl-0. l^B^J 

TE«P7=^TEMP6^TENP6+TEMP5^TEMP5^TEMP2 

R1=SQRT<TEMP7> 

TEMP6«TEMP6+S 

TEMP7=TEMP5-TEHP0 

TEMP7*TEMP7^TEMP7+TEMP6^TEMP6^TEMP2 

R2^^SQRT<TEHP7> 

TEMP7»S-^TEHP9 

TEMPS* < R1-^R2+TEMP7 > / < Rl-»-R2-TEMP7 > 

SUM==SU 1+RL 0G<TEMP6> 

1001 CONTINUE 

1002 CONTINUE 

1003 CONTINUE 
TEMP8*S 

IF<S. LT. 0. 1E-2>TEMP8*T 

COEFF^ < SUM*©. 89877E10 > / < TEMP8* I MAX*** JNRX>**KMRX > 
GO TO 1005 

1004 TEMP1*SQRT<TEMP1> 

COEFF*0. 89877EJ.0/TEMP1 

1005 RETURN 
Er^D 
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A.2 PROGRAM PS072SGF TO DETERMINE SENSOR SIGNALS DUE TO CHARGED PARTICLES 
AND RESULTING aECTRONICS RESPONSE 

A.2.1 Sunmary 

This program uses the particle path characteristics calculated by 
program P5072CHG to generate sensor signals for selected particles and then 
determines the response of the electronics to these signals. The particle 
parameters of mass and charge are selected either by Input card or by a ran- 
dom number generator. The output from the electronics model Is stored on 
tape for future analysis and may be plotted. If desired, by selecting the 
proper code on an Input card. 

A. 2. 2 Description 

A. 2. 2.1 Determination of Velocity 

The starting velocity Is one of the parameters supplied by Input 
card. The new velocity of the particle at the end of each Incremental step 
Is determined from the new particle energy and the particle mass, using the 
relationship that energy equals half the product of mass and velocity squared. 
The new energy Is determined by subtracting from the starting energy the work 
done In traversing the step distance. The work done Is calculated from the 
potentials. Program P5072CHG provides two potentials at each step. One Is 
the potential due to the applied potentials EPOT(J) and the other Is the po- 
tential, per unit charge, due to the particle charge CPOT(J). The work done 
between two points Is equivalent to the product of the potential difference 
between the points and the charge. Thus, the work done Is determined from 
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Q*EPOT(J) + Q*Q*CPOT(J) calculated for the two steps* J and (J-1). The tine 
taken for the step Is determined fron the average velocity over the step and 
Its distance. 

A.2.2.2 Determination of Grid and Film Currents 

The currents are equivalent to the time rate of change of charge. 
Program P5072CHG provides the total charge on each film and collector grid 
element at each step. The current Is found from the difference in charge at 
the beginning and end of the step, divided by the time interval, calculated 
above. 

A. 2. 2. 3 Film and Grid ID Thresholds 

The two systems are Identical except for signal polarity, so only 
the grid ID will be described. The Input to the linear amplifier after the 
Input circuit Is calculated for each collector grid strip using the ramp 
function response equations arrived at by Laplace transform. The Input to a 
threshold detector Is the sum of the amplified signal from the Impacted film 
plus the factored Inputs from the other films applied as analog Inhibit sig- 
nals. 

The time to reach threshold Is determined by performing a linear 
Interpolation between the present and most recent steps. Thresholds at other 
collector grids are only permitted If they occur within 0.2 microsec of the 
first ID. 

The time of the first ID, either film or grid. Is used as the start 
time for the PHA measurement in the electronics subroutine LES. 


A-55 



BSR 4234 


A.2.2.4 PHA Anpllfler Signals 

The Inputs to the four film amplifiers are sunned to be used as 
Inputs for subroutine LES. To save computing time, data points are not ac- 
cumulated until the Input reaches one-tenth of the input threshold level. 

A.2.2.5 Program Flow 

The first Input card is read to retrieve the parameters which select 
the various options. The number of particles to be analyzed, the velocity of 
the particles, the mass and charge, if random numbers are not used, the num- 
ber of the Input data set for path data, and the number of the output data set 
and codes for printout selections are retrieved from this card. The second 
card gives plot axes dimension Information. 

The first step Is to define the physical stopping point within the 
sensor as either the film or east sensor shield followed by the Initialization 
and setup of the Cal Comp plotter. This setup can be bypassed if no plotting 
is desired. Next, the particle path data are read in from tape as the potentials 
and charges at each part* "e position, plus a header record which defines the impact 
position on the sensor relative to the center and the total number of data points. 

If random particles are to be selected, the first mass and charge 
values are calculated. The random number generator scales the values deve- 
loped so that they fall within a range specified by the axes dimensioi. infor- 
mation given on the second input card. All random numbers generated are <!sed 
to save computing time over the method which uses all numbers for deriving 
masses and charges and then rejects those which do not fit the problem. The 
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variables are Initialized, and the particle position at nhlch Measurable sig- 
nals cm be detected Is then detnnined ffxr use as the starting point for the 
rcnalnlng calculations. Ihe steps start at 10 aeters froM the sensor and step 
In rapidly until either the sl^l reaches one-tenth of tl e threshold of a 
collector grid or flUi circuit or a point 0.4 centlaeter fron the suppressor 
grid Is reached. This Is done to give a starting point for the potential 
■easureaents and the calculations of work done on the particle, since absolute 
potentials are Measured relative to infinity or a point of zero potential . 

The sensor currents at this position are written out If the selection 
code deaands thea, followed by the calculatlois of work done md the Magnitude 
of the reaalnlng particle energy. Providing that the renalning energy Is posi- 
tive, the new velocity md the tine Increaent are detenilned, followed by a 
calculation of the new flla and collector grid currents. 

The flla and collector grid threshold ID status Is then deterained 
together with the value of the PHA amplifier input signal. 

This sequence Is continued until all particle positions have been 
analyzed or the remaining energy reaches zero, indicating that the sensor 
forces have stopped the particle. 

When the sequence Is complete and if a film 10 has occurred alone, 
before a collector ID or less than 1 msec after a collector grid ID, the data 
are passed to subroutine LES, which calculates the electronics response. 

The results of the electronics analysis, nanely the PHA, film and 
collector gtid ID and accumulator counts, together with the particle charge, 
mass, and velocity are stored for future analysis and. If required, the points 
are plotted. 
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The prograi then returns to the start to read the next selection. 
A.2.2.6 Subroutine LES 

Section 3. 2.1. 2 describes the operation of the stq>1e aodel of the 
electronics. The Modeling Is accoHpllshed In subroutine LES. 

The prograM uses the data points passed to It by dunay argunents 
and sinllerly returns values for the PHA count and accunulator count. 

The output signal fron the sensor Is a pulse uhose length and ampli- 
tude are determined from the path characteristics and particle characteristics 
in the MAIN program of P5072SGF. The signal Is In the form of discrete ampli- 
tudes at discrete times. This subroutine treats the signal as a series of ramp 
functions by developing straight-line equations for the signal between adjacent 
values. 

The program then evaluates the slope of the ramps to determine the 
status of the two switches. The result of this evaluation determines which of 
three subroutines will be used to calculate the value of the output signal. 

The subroutines called are COND 1, COND 2, and COND 3, which calculate the 
responses using predetermined equations that were arrived at by using the 
Laplace transform technique. A fourth subroutine, CVOLT, Is used to calculate 
the voltages across the capacitors at the end of each step, as these are re- 
quired as Initial conditions for the next ramp function. 

A.2.3 Method of Use 

All references to Job Control Cards (JCL) are for the IBM-370 system. 


A- 58 



BSR 4234 


Operation of the progran requires a alnlNn of two Input cards, plus 
one data tape produced by P5072CHG giving particle path data for the path and 
sensor to be analyzed. The nlnlaun Input allows, on the one path, either: 

(1) analysis of one particle with Its wass, charge, and velocity selected by 
input card, or (2) analysis of any nuaber of randomly selected particles, all 
at one selected velocity, up to a maximum of 999 particles. The results will 
be printed and optionally plotted. If more than one particle is desired in (1), 
different random numbers in (2), or different velocities or different paths 
are desired, then additional sets of cards must be added with the new codes 
and the appropriate path data sets must be available on tape. 

The information required on the cards is as follows: 

Card 1 

Column Requirements 

1-3 A number from 1 to 999, format 13, representing tte number of 

particles to be generated. If discrete particles are selected, 
the value should be 001. 

4-6 A ninnber, format 13, which determines the rate at which tbe 

element charge values will be written out, e.g., if the value is 5, 
every 5th step will be printed out during analysis. 

7-13 A number representing the particle/s initial velocity. Format F7.2. 

14-23 Not used. 

24-29 An odd number used to start the random number generator. Format 
16. 
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Card 1 (Cont.) 

Coliwn Requirtwents 

30-32 The nunber of the particle, e.g., 5, for the 5th particle of 

the total set of randoai particle'* generated for which the charge 
data are desired. If zero, all particle data will be selected on 
the basis of the nunber In colmns 4 through 6. 

33-42 The particle charge, format E10.4. If random numbers are selected, 
this value may be blank. 

43-52 The particle mass, format E10.4. If random numbers are selected, 
this value may be blank. 

53-55 A number which If greater than 10 will cause random particles to 
be produced and a plot of the results generated. If greater than 
1 but less than 10, random particles will be generated. If less 
than 1, discrete values must be put In columns 33 through 52. 

56-58 The Input data set number which matches the JCL card, e.g., 
FT12F001. 

59-61 The output data set number. 22 and 25 must be used for shielded 

film data sets. 

Card 2 

Card 2 provides Information required by the CalComp plotter to set up 

the axes and by the random number generator to set up mass and charge values. 

The axes charge and mass Information Is developed as follows, with references 

being made to the following figure. 
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Consider tte charge values: because of the range of values, the logar- 
ithm of the charge Is plotted on a log scale. The distance along the y-axis Is 
given by QC0N*1og (log Q* where QCON Is a constant and Q Is the charge. 

The random number generator develops numbers between 0 and 1.0. The 
value of 0.5 Is subtracted «.o give a range of -0.5 to ■•'0.5, which Is then 
multiplied by QRG, t\ . desired range of log Q values. We now have the cor- 
rect range cente id about the origin. The mean value of the desired range, 

QMN, in Inct s from the origin, Is added to the generated value to place 
the rang^ In the required :.rea. The value obtained (y) Is the position along 
the log Q axis. In inches, of the desired log Q. 
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Therefore, scale value ■ exp (y/QCON). 

To obtain log Q, we now nultlply by the scale factor, QFR, then 

Q = exp (log Q). 

An Identical procedure Is followed for the nass N. 

The values required by Card 2 are: 


Column 

1-4 


5-8 

9-16 

17-26 

27-30 

31-34 

35-42 

43-52 

53-56 

57-60 


Requirement 

QRG, range In Inches of required log Q values on the plot (Format 
F4.2). 


QMN, mid-point of range In inches from origin (Format F4.2). 

QFR, scale factor (value of log Q axis at origin). Ignore any 
minus sign (Format B.2). 

QC(W, axis constant for size of axes to be plotted (Format F10.8), 
i.e., axis length for one cycle In inches (CYC) = QCON* In 10. 



AXLEM length 


the same as the equivalent log Q definitions, 

M. 


of the log Q-axIs In Inches (Format F4.2). 
of the log M axis In Inches (Format F4.2). 


for 


If a plot Is not required, a card must be submitted but It may 
be blank. 
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A.2.4 Flow Charts and Program Listings 

A flow chart for program K072SGF Is shown In Figure A-2. 

Program listings are provided for program P5072SGF and subroutines 
LES> CONDI. C0ND2. C0ND3, and CVOLT on pages A>67 through A-77. 
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P5072S6F 


MF IS DflTft SET REF. NWBER FOR PftRTlCLE PATH DATA 
NDOUT IS DATA SET REF. NUI«©ER FOR OUTPUT DATA SETS AS FOULOHS 
POSITIVE CHARGE : UP 4 EAST NORMfH. PATH 20 

l«ST SENSOR 21 

EAST SHIELDED FILM 22 

NEGATIVE CHARGE : UP 4 EAST NORMAL PATH 23 

I^ST SENSOR 24 

EAST SHIELDED FILM 25 

NOTE : REF. NOS. 22 4 25 MUST BE USED FOR SHIELDED FILM DATA SETS 


DIMENSION DRTA<2, 100©>, IBUF<1000>, ^ANK<8>, CC<4, 500>, CF<4, 500> 
DIMENSION IG<4), IFM<4>, C1G<4>, C.2G<4>, C1F<4>, C2F<4), VCG<4>, VCF<4> 
DIMENSION VIDQ<4>, TIMF<4>, TIMG<4>, S'R1F<4>, VRZ VR1G<4), VR2G<7> 
DIMENSION VIDF<4>, CPOT<500>, EPOT<500>, DIST<500> 


\ 


ORIGINAL PAGE IS 
OF POOR QUAm 
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Rem. N 
INTEGER OUT 

^<6^«se>«0. tm «FR. 

IF<M>OUT. EO. 22>SHLD-5. 4864E-3 
IF<NDOUT. EQ. 25>SM.D>^. 4864E-3 
IF<KODE. LT. ie>GO TO £03 
DEFINE PLOT f»<ES £ TITLE 
CBLL PLOTS<IBUF,10OO,0> 

Cm.L PL0T<1. 0, 1. 0, -3) 

CVC**2. 30259<<<MC0N 
DTV*^/'CVC 

CRLL LQflXIS<0. 0, 0. 0, £H'-LN 
CVC-2. 30259<<<QC0N 
DTV=1/X:VC 

CftLL LQfiXIS<0. 0,0. 0,6H-LN 
Vm<=RXLEQ+0. 5 

SlL 'ooo 15HVEL0CITV <ti/S> 

CflLL NUMBER<.999. , 99:. . , . 14, VEL, O 0, -li 
FINISH AXES DFJ^INITIG,< vtu, e. U, 1> 


.'1, “€, ftXLEM, 0. 0, MFR, DTV> 
Q, 6/ RXLEQ, 90. , QFR, rrv> 


- 0. 0, 15> 


CHARACTERISTICS 


C 
C 

C PARTICLE PATH 
603 REMIND NF 

RFAD<NF>X, V, NPTS, BLANK 

685 SS?[!Sf lO' '-“I- ■<>. <CF<L, K>, L-1. 4> 

60S* DO 4000 KP«1, NUN 
MRITF<6, 6035)KP 
IF<KODE. LT. DGO TO 625 
610 IV-IX»65539 

iF< IS’>615, 6:16, 616 

615 lVs=IV+2147482647+l 

616 VFI.<^IV 

VFL^VFl.<«. 4656613E-09 
IX-IV 

VFL*-VFL-0. 5 
RANDQ<-VFL>«QRG+C!MN 
QSV^FXP < RANDCVQCONJ 

QVAL^^QSV^ < -QFR > 

Q-EXP<OVAL> 

616 IV=-IX<«65539 

IF K TV^61 9, 620, 620 
61.0 IV-IV+214748.3647-<-l ' 

620 VFI.- IV 

VFL«VFL>«, 465661SE-09 
IX<-IV 

VFl -VFL-0. 5 
RAI’4i'.<M- VFl.;«iaRC+MMN 
MSV- FXP RANDIl,''MC(tN > 

MVAI 'H'^sV+C-HKR:* 
n=‘Fxp<;wvrl;.' 

INJTImI J2ATT0N 
625 TIMF< 0. 0 
Nr<i<^0 

N( I 
NFlf« 0 

Nt<(il< 009 
HI>FTs '='99 
VFl I' VFL 
I P -0 

DO 700 M,4 


OBIGINAI! PAGE IS 
OP POOR QUAury 
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IFH<1>«0 
IQ<1>«0 
CiQ<l>^. e 

C1F<I>«0. 0 
VCG<1>«0. 0 
VCF<I>«0. 0 
VUXKI>«>0. 0 
VIDF<I>*0. 0 
C2G<I>«0, 0 
C2F<1>«0. 0 
700 cmriHLM 
ENRGV=0. 

c ccTERHiwE srmiim coNDiiums 

DO 726 J«2.12 

IF<D1ST<J>. LE. 1. 369£-2>GO TO 730 
» JJ^J 

TINC-<DlST<J-i>-DIST<J»/VEL 
DO 72S K«;L,4 

C2G<K>«Q#<CC<K/ J>-CC<IO J-1>>/TINC 
C1G<K>«C2G<K> . 

C2F<K>«0^<CF<K. J>-CF<K. J-«1>>/TINC 
ClF<K>s^F<K) 

1F<C2Q<K>, LT, -0. 4E-9. OR. C2F<K>. QT. 0. 4E-9>G0 TO 730 

725 CONTINUE 

726 CONTINUE 

730 POT1 -Q**^EPOT<JJ-l»Q^QHCPOT<JJ-l> 

735 Js^JJ+ND 
LP»LP+i 

IF<IPt4UI1. EG. 0>GO TO 736 
IF<KP. NE. IPNUM>GO TO 738 

736 IF<LP. LT. OUT>GO TO 738 

737 IP^0 

MRITE<6. 6040>C2Q. C2F/ TIME/ DIST< J-l> 

738 IF<DIST< J>. GT. SHLD>60 TO 739 
DfiT0<i/ Jr>«TIME 
D0TH<2/ J>=0. 0 
GO TO 2500 

739 ND=^ND+1 

POT2 =^Q^EPOT < J > Q^CPOT < J > 

W0RK=:^P0T2-P0T1 
REM-ENRGV-WORK 
IFCREM. GT. 0. 0>K J TO 745 

740 IF<DIST< J>. GE. . 67564E-2>G0 TO 741 
WRITE<6/ 6050> J 
GO TO 2500 

741 IF<DI?:T<J>. GE. . S7Cv8E-2>G0 TO 742 
WRTTE<6/6060>.; 

GO TO 2500 

742 NRITE<6/60?0>J 
GO TO 2r>00 
VEl 2-S0Rl<2^PKM/M> 

RVEI,‘- < VEUl +VHL2 > /2 
TINO 1-1 >-DJST<J>>/HVFL 

TlME-^-TJHF+TTNC 
VEL:U VFL2 

C CP! GUI RTF 
DO 8fi0 
G2G<K> 

C2FCK> 

800 COUTINUE 
C CCiU.FGTOF: GRID JD 
DO nOOO K^i/4 
SL P « < C 2 G < K > -- C 1 G < K > > / T I NC 

SI. Mil ^ 0. 5^ EXP < -2. 27272727:iK3^^T ) NC > < 2. 2E5^C1G < K > VCCi < K > > 
fUH:^^<:l-EXP<-2. 272727273E3=4=TJNC>>%<2. 2E5^C:tG<K>-48. 44;SLP> 
HIK.:.. 2. 2Er»rf.| P=fTTNr 


NFN GRir* ftWD FILM CURRENTS 
K^ l/ 4 

-(MiCCXK. J>-CC<K/ J-i>>/TINC 
^ 04 <CF<K.^ J>"CF<K/ J-1> >/T1NC 
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. RlQ<K>*^SUMl-*-SUM2-^SUM3 

SUH4*=0. 5>«‘EXP<-2. 272727273E3*«‘T1^K;>%<2. 2E5^CiG<K>-VC0<K>> 
SW.3«SLP»*'48. 4*»*<l-EXP<-2. 272727273E3*^TINC>> 
VR2G<K>«SUM4-^SUM5 
C1B<K>=C2G<K> 

VCG<K>=VR1Q<K>-VR2G<K) 

leee continue 

1F<ND. LE. NDGT+DGO TO 1016 
IF<TirflE. QT. TIMQ<1>+2E-7>Q0 TO 1106 
1610 1F<NQID. EQ. 4>G0 TO 1100 
VR2G<5>«VR2Q<1> 

VR2C<6>«VR2G<2> 

VR2G<7)«VR2Q<3> 

DO 1100 K«l,4 
IF<IG<K>. EQ. 1>G“ TO 1100 

VID«15. 75^W2G<K>-1. 47**^<VR2Q<K+l)-fVR2Q<K+2>-^VR2G<K^2>) 
IF<VID. LE. 12E-DG0 TO 1650 
VIDCi<K>:=^^VlD 
GO TO 1100 

1050 IFCNGID. GT. ©>GO TO 1660 
NDGTs=ND 
IG<K>^1 
NOID:^=l 

TIMG<1>=<. 12E-1-VIDQ<K>>^TINC/<VID-VIDG<K>> 

T1MG<1>«TING<1>+TIHE-TINC 

IF<NFID. QT. 0)GO TO 1160 

TIHID-TII1Q<1> 

GO TO 1100 
1060 NGID=-NGID+1 

TJHG<NGID>«^<. i2E-l~VIDQ<K>>^TINC/<VID-‘VIDQ<K>> 
TIMG<NGID>=T1MQ<NGID>-I^TIME-T1NC 
IF<TIHG<NGID>-TIMQ<1>. LT. 2E-7>IG<K)«1 
1100 CONTINUE 
C FILM ID 

DO 2000 K«l/4 
SLP^.<C2F<K>--C:lF<K>^/TINC 

SUHl-0. 5^EXP<-2. 27272727:NE3HcTINC>^i'2. 2F.54:C1F<K>+VCF<K> 
SUM2=^<i--EXP<-2, 272727272E^<♦TIr4C>>'^^<2. 2E5+ClF<K>--4©. 4^SLP> 
SUMr^=2. 2E5^SLP^TINC 
VRiF < K > =^SUMi+SUt12+SUM3 

SUI‘14^^0. 5’*^EXP<“2. 2727272V2E:i^jmcy^<.2. 2F54:C1F<K>-VCF<K> > 
SUMCr- SLP>^4a 4=* <l-EXP<-2. 272727273E3^TINC> > 

VR2F < K > ==SUM4+SUM5 
C1F<K>«C2F<K> 

VCF < K > -VPJ F < K > - VR2F < K > 

2000 CONTINUE 

1F<ND. l.E. NDFT+l>GO TO. 20:10 
JF<TiriE. GT. TIMF<1>+2E-7>G0 TO 210© 

2010 ir<Nr EG. 4>G0 TO 2100 
VR2F<5>=V <2F<:1> 

VR?F<6>=-VR2F<2> 

VP2F<7;«=:VR2F< ^ > 

DO 2J00 4 

1F<IFM<K>. FO. 1> GO TO 2:100 

7'.>VR2F<K *>+:!. 47-+ <VR:>F<K+1>+VK.2F<H>2>+VR:4’i:1 
IKVID. LF. :1;H-0 GO TO 2050 
VTDF'.J>^-VID 
GO TO 2:tM0 

20r.0 1K<NF JD. GT. Ci>GO TO 
Nl »F T * ND 
JRKIO^ *J 
NFIDvl 

T 1 HF < 1 > 4 < . 1 2H -i -V 1 \ *F < K .> > ^T 1 NC/ < V I D- V I DF < f O > 

1 1 I F* < J. > 1 HF i. .1 > +T I MF-T I NC 

IFvNGTD. GT. Ci>GO TO 2100 
T1HID> ^T;HF«a> 
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W TO iiM 
. 2060 NFIIMnO^i 

TIHF<WFl1^>-C lStt-l-VlCF<K»*TJHC/<Vll>-VH>F<IC> > 
Tiir<trn»<*TfHForii»^Tnc>TtNc 
IF<TIir<NFIO>''TItr<l>. LT. 2E-7>1FN<0«<L 
2100 CONTINUE 
C CALCULAYE PHR ««.UES 

- SVMLF«ViaiF<i:MVMJ^<2>-»VRftF<3>«VRlF<4> 

IFCNDT. onr. TO 21SO ' 

IF<SVR1F. UT. . 1E>3>GD TO 735 
MITfKl* 1>«T1ME-TIND 
pOTO<2.;i>>€l. • 

DATf»l*2>sTItC 

MmK2*2>sSVMF - ' - 

N0T»2 ^ J- 

00 TO ^ 

2150 NDT«4i>T-i4. 0 

DOTR<l«NDT>«TinE 
DATA<2* NDT>»SVIfclF 
GO TO 735 

2500 IF<MFID. EO 0>G0 TO 2700 ^ 

IF<NG1D. EQ. e>00 TO TT'O 
1F<TIHB<1>1-. IE-2 LT. iHF<l>>GO TO 2700 
2510 ETIHE=TIHH>*. lE-2 


3000 COLL LES<DAT(LNDT,ETIHE. Oljrr.- KP.. iPHUH<NPHfLHnCC> . 
IF<KODE. LT. 10>Q0 TO 3560 
C PREPflRE TO PLOT PHR VRLUE 
PHR=M PH R ' 

C PLOT VfltUE . ' ' 

CRLL SVHBOL <RRM>I1.. RffNDQ, . 01^ 75. 0. 0» -1> 

CRLL NUMBERO90. > 999. , . 07. PHR. 0. O, -1> 

lF<t«CC LT. 2>G0 TO 3500 

CALL SVMB0L<999. . 999. , . 03. 125. 0. U. -1> 

GO TO 3500 
2700 HRITE<6. 7000> 

35^ fiCM«S<0) 

RQL»ALOG<RQ> 

RML-RL00<H> 

I0T0T-1G<1>+1G<2>*2+IG<3>*4+.Q<4>*8 
1'-T0T«IFH<1 IFM<2>*2+IFM<3>»4+IFH<4)*8 

HRITE<NT>0UT >VEL. O. M, IGTOT.. IFTOT, NPWb NAOC 
HRITE<6. 6O80>O.H. IG, IFM 
HR1TE<6. 7010>fteL, RHl. 

4000 CONTINUE 

IF<KODE. I.T. 10>6O TO 4500 

CRLL SVMBOL<20. , 0. 0, . 14. 11, 0. 0. -1> 

CRIJl SVMBul.<20. . J0. , . 14, 11, 0. 0. -2> 

CRLL PL0T<22. ,-l. , -3> 

4500 GO TO 600 
5000 COf'ITIHl.lF 

IF<KCC€. LT. 10>GO TO 510© 

CRLL fn.OT<e. , 0. , 999> 

5100 CONTINUE 

6000 F0RMRT<2I3.. F7. 2, E10. 4, 16, 13, 2E10. 4. 313> 

6010 F0RHRT<2<2F4. 2, E8. 2, FIB. 8>, 2F4. 2> 

60.^0 F0RMRT<1.X, 2<2F5. 2, E8. 2, Fll. S>, 2F5. 2> 

6030 «^0RMRT<5X, 'VELOCITY ■= E10. 4, 2X, 'X ORO = ', F8. 3, 2X, 


“k 




OF POOR 


ORD “ %F8. 


NO. OF STFFS = ^ 7 1 4.. 2:X> RRNDU = M6/ 2X. -'FfiTH DRTh SET « ^ 

ninpin dhtr set » 

60/r» FORMRK O 7 "PRRriCI h NUMBER 


6040 r.>:NRT :i.X. CURRENT : GRIOl 7 FB. 2, 4X> ' OR 1 D? = ^ F9. 2^ 4X/ ^UR1U3 
4X. GRll'4 =. .. FB. X. FILM.t 4X,. Fll «= 

r>X. FlIM ^ R.4>:. T-!IM4 • 7 KO. i?, 4;:.. time u ,F9. 3, 4X.. ^DIST 

^ 4 ^F9. 

60>0 FfiRHm^'iX. 'PHRTICIF STOPS FnRF FHM RT -4KP 
6060 FiTRWm<5X.* PRRIjriF STOPS BEEliRI: bRIO t^T Sll P 2:;> 


STOPS E^PKORE SUPPRESSOR RT STEP S ?>:. T3/.^ 

FORMAT V .X, '(HRROF - '',F9. 3, ' MASS = ;2. • Fll. 

5M JD :• ','U2/> 

V.iufi F0RMHT<4X, ' SIGNAL RFl CM IHkFSHOl.l) V> 

-.•■i.lO FORMAT <5;;, TOG 0 - F9. ? T-X, 'LOG M *■• ',E9. 3.'.’ 
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SUBROUTINE IMS <MTfW NO« ETIIC^ OUT> IG IPNUH^ 

C • 

.C REVISED" 15 DEC iS75 
C 

" DtNENStON DNT^2# lBBe> 

COHHOll VCI#Vd^VC4^VC5.yCB#TCl^Tt:2.TCl#TC4^TC5>TC6#VBLl^VBL2^X^SLF 
INTEGER SNBB^OUT 
VRl--a0 

6 V 

VR^^e 
ViU12>=B . 6 
8 

SNRBM 
1VRS«2 
NBCOa 

vci«=a8 "" 

VC1=8, 0 
VC4«0.0 

vc5==^a 0 

VC8==0. 0 

TCi< 454545-I5-I5E4 
TC2==. 350140056E4 
TC3-. 494S845522E5 
TC4=^ 10r.2ft67642E3 
TC5=-\ 21 9r%22615E5 
TC6^ . 72G21894 
J^-^8 

TJHE-C>ftTft<l, ND> 

L=NI>+1 

DO 2020 

TF<TIHF. GE. ETIf^>60 TO 2925 
Tjr^: TIHE+2H-6 
KOiJNT=KOUNT+l 
DflTB<l. I>^TIME 
2020 COr>ITINliF 
2025 L- t4D-Ki 

IF<KOUrTf. LT. L>GO TO 203© 

DO 203tO I:=^L,K0UNT 
DfiTft<2.. 0^0. 0 
2O:<0 C0^4TINllF 

- IF<K WF lPNij|1>G0 TO 2O80 
IF<OUT. GT. 5O0>GO TO 2000 
KH^1TF<0, :ilO0> 

20o^?^ me ? 

2Jri0 IKir4C. FO. KOUNT+J>GO TO 2OO0 
T2r--DrtTFK1. IWr:> 

VRI :0 ^RIR<2> JNC-1> 

VR» DRTR<2/ I HO 
INC IHC+1. 

(. 

C CiRLiVM RTF JHKn SI OFF <SI..P) 

C 

VRI 2-VHI .1 
X- 12-Tl 
SI F- V/X 

JF':SI.P. I T. a OKiCi TO 2:1if:0 

I*. 
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CALCULATE CONDITION i VALUES 

Sd CALL CONDI <VA32# VRl# VR2^ VA9^ VA12> 

CALL CVOLT <VR32^ VR±> VR2/ VR9> VR12> 

KONCM=i 

GO TO 2600 

POSITIVE SLOPE - CALCULATE SNITCH STATUS 


2160 VR202^VR2 
VR112^VR12 
VR132=^VR3t2 

CALL CONDI <VR±32. VRiei> VR102. VR109, VR112> 
CALL COND2<VR202^ IL C> 

RVR1*=VR101-VR1 
RVR2‘--VR202-VR2 
IF<RVR1. LE. RVR2>G0 TO 2178 
SNf®^SfWB“l 
50 TO 2180 
2170 VR2=^VR1 

2180 fiO TO <2181, 2iai> 2184, 2186>, SWAB 

>4 ct4 crnp i fii 

2184 IF<VRi2. EQ. 8. 8>G0 TO 2298 
VR282-VR2 

CALL C0r«>2 <VR202, VR289, VK212> 

CAI.L C0t«>3 <VR209, VR212> 

RVR2=VR202~VR2 

RVR12=^VR212-VR12 

1F<RVR2. GT. RVR12>60 TO 229© 

60 TO 2200 

2186 1F<VR12. EQ- 8. ©>G0 TO 2298 
VR112=^VR12 
VR122=VR22 

CAIX CONDI <VR122, VR181, VR102, VR109, VR112> 
CAI.L CCW4D2 <VR209, VR31*: :' 

RVR2~VR102-VR2 
RVRil 2=VR2:l 2-VR12 
IF<RVR2. GT. RVR12>G0 TO 2290 
GO TO 2158 


CAl-CULATF VRI UF.S IM»F.R C0r€>ITI0N 2 


2288 CALL C0ND2 <Vk2> VR9, VR12> 

CRI J. C0r4f»l <VR32* VRi, A,R.C> 

CRI I ( VOLT VR32, VRI, VR2, VK9. VRI 2> 
KOr^l^ ‘2 
GO *f(i 2680 

2250 IF<jNr. EO. KOUNT-^J>GO TO 2888 
T2“-.DRTri<l, ir4C> 

Tj=^DRTH*:i. 

VALl^ DRTR<2, INC-t) 

VAl.2“-r*ATR<2> INO 


CRI I. CONf^J VR1R1, VRilH?, VRISS, VR.I.IJ' • 

Cfil , I C ONI '2 '• VR^ Hi,', M/ B > 

TBCVk:1Ci.t . l.F VR;;tV> SNRH- SHHn-H 
(lO 10 

/ 

T CRMiiiRrr- vHUii-r*. unh-r coNomoH / 


INC- IHC+:1 

V-VRl.i?-VrtL;t 

X'-TSJ-U. 

SLP-V/X 

VR?&?.~WR2 

Vkl .t?'- VR;l? 

VkJ \>R?ri 
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2298 SHRB-SHRB-2 

2368 C8LL C0NDKVR32# VRjb VR2> B> 

CBLL C0ND3<VR9. VR12> 

CBLL CVOLT <VR32. VRjU VR2> VR9^ VR12> 

K0KD=3 
GO TO 2668 

2358 IF<INC. Ea KOUNT-M>GO TO 2888 
T±i^i)RTft<l> IHC-1> 

T2*^DBTR<1^ 1MC> 

ViRLl=DftTft<2, IHC-l) 

VRL2=^DBTft<2^ INC> 

INC>^ir4C^l 

V*VftL2-V«Ll 

Xs=T2-“Tl 

SLP=V/X 

1F<SLF. QE. 8. 0>GO TO 2378 

VR282=VR2 

VR112«VR12 

VRA32»=VR32 

CBLL C0N>1 <VldL32/ VR101 . VR102^ 

CBLL C0ND2<VR282> C> 

RVRl-VR101-VRi 
RVR2^ VR282-VR2 
IF<RVR1. LE. RVR2>G0 TO 236© 

60 TO 2378 
2360 VR2=VR1 

VR4,32:=-VR32 

2378 GO TO <2388^ 2375> 2380, 2375 >> SNRB 
2375 CRIX COr«^ <VR132, VR181, V*R102, VR189, VR112> 
CRLL COt4D3<B, VR312> 

IF<VR101. LE. VR312> 60 TO 2378 
GO TO 23G© 

2378 SWRB-SI48B+2 
60 TO 2150 
2388 VR202- VR2 

CfilX C0rC»2<VR202.. R 8> 

CftLL C0ND3<8> VR312> 

IF<VR2©2. LE. VR312>60 TO 2385 
GO TO 240© ' 

2385 SWRB=SWftB-^2 
GO TO 2280 
2390 SHriB-SWRB-2 
C 

C CBLCULRTE VRLUES UNDER CONDITION 4 
C 

240© CRLL C0NDKVR32.» VR1> R/ B, C> 

CRLI- C0ND2<VR2.. R. B> 

CRLL C0Nr>“<VR9. VR:12> 

CRl.L CVOLT *:VR32. VR:l. VR2. VR9. VR12> 

KUND-4 
GO TO 260© 

240© IF<1NC . F.Q. KOUNT+:i>GO TO 2800 
Ti DRTR<1. I?4i>:l> 

12 DRIR<:b 1N»: > 

VRI..:l- DRrri<2. TN<>:1> 

VRI 2M tRTR*:?. I NO 
INC lNC+:t 
V^VRL2“VRI-:1 

SLP-V/X 
VR202^ VR2 
VP J J 2- VR J 2 
VR.I .<2“-.VR32 

CRLI CONDi < VRA32, VRICH, VR:1.02> VR109, VR13 2> 
Chi I ( i:iND2<VR2n2.. fh B> 
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CHLL C0ra>S<fhVK3iSt> 
tFKVmML L6. VR202>QO TO 2S00 
1P<VR202. LE. VR312>Q0 TO 2475 
00 TO 24M 
2475 SM fO»- SMf>B»2 
00 TO 2200 
2500 SMW D»SUnP »l 

IF<VR101. LE. VR312>G0 TO 2550 
(jO to 2300 
2550 S HnO«5MnD »2 
00 TO 2150 

2600 IF<IPNUN. ea 0>Q0 TO 2601 • 

1F<K. NE. 1FNUH>G0 TO 2603 

2601 IF<J-0UT>2.^3« 2602, 2602 

2602 

MRITE<6. 3O00>T2, VR32, VRl, VR2, VR12, VR9, SHRB 
GO TO 26^ 

2603 J^J*± 

2605 IF<VR9. LE. :l£- 01)G0 TO 2620 
00 TO <26x0, 2T‘50>, 1VR9 
2610 PT2=T2 
1VR9«2 

TPHfl=PT2H>Tl 
NUM=TPHR/ 4E-04 
NPHfl=NPHft+NUM 
GO TO 2750 

2620 T1HID»ETIHE-. lE-2 , 

IFCT2 LT. TiniD>GO TO 2750 
GO TO <2630, 2640>, IVRO 
2630 lF<iriTC. LT. K0UNT+1>G0 TO 2750 
PT2»=T2 
TPHfl--PT2-PTl 
NUH»TPHft/. 4E-04 
NPHft=^NPHft+WJH 
GO TO 2800 
2640 PT1«T2 
lVRSi=l 

IF<NftCt:. EQ. 0>NPHft--=l 
NflCC'-NfiCC+l 

2po GO TO <2100, 2250, 2350, 2450>, K0»4> 

2800 CONTINUE 

WRI TE<6, 3700>NPH8, NftCC 
3000 FORMflK?}:, 6<E12. 6, 6X>, 5X, I2> 

'• “><' »=<- »x. 'v 

3700 0'', 5X, 'PHft = IS, 10X, ^flCCUMULRTOR COUNT » •*, I3> 


ORIGINAL PAGE IS 
OS' POOR QUALnY 
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SUBROUTIME CONDl<VKI2# VIdL. Vl^ 

HEVl^ 2€ SEPT 1S7S 

COmON VCI#VC1>VC4>VC5.VC8^TC1>TC2^TC3^TC4#TCS.TC6^V<«^VAU^^X/SLF 
REfIL HR12#HR32 
INTEGER VS1AB/VS2RB 
VSlB»-104iVR32 

VR32=< VNLl-VCl >^EkP<-TCl^X>-i-SLP%<l-EXP<-TCl*>X> >/TCl 

VS1B=-10’<^32 

VS1N&4 

XE<aBS<VSUR>. LT. 5. 8>Q0 TO 288 

vsift»siGN<5. vsia> 

V51RB»VSiRB^l 

280 ir<^»S<VSlB>. LT. 5. 8>G0 TO 488 
VS1B»S1GN<5. 8, VS1B> 

VSlRB^VSlilB^a 

480 GO TO <500# 608^ 680/ V^IRB 

500 VRl^::. 35780487E-2^<VCI-VBL1>^<TC1*«€XP<-TC1^X>-TC2’*€KP<-TC2>»X>>-. 957 
188487E-2^a.P^<EXP<-TC2<^X>-EXP<-TCl>i'X>>-VCl^EXP<-TC2’i^X) 

VR2^VR± 

IF<VR2. LT. -5. 8>VR2=-5. 0 

IF<VR2. GT. 5. 8>VR2»^5. 8 

VS2R=VR12 

VR12=^VR2 

VS2B-VR12 

IF<RBS<VS28>. LT. 5. 8>G0 TO 53:0 
VS20:^SIGN^5 8/ VS2A> 

VS2AB=^VS2AB-H 

530 IF<AB5<VS2B>. LT. 5. 8>G0 TO 558 
VS2B=SI6N<5. 8, VS2B> 

VS2AB:^VS2RB+2 

558 GO TO <568,57O.570/570>, VS2RB 

568 VR9-<VC8*3. 3E--8-VC4%6. 8E-6>^EXP<-TC4>t^X>/6. 833E-6-^8. 99517^<18=*'<VCI- 
2VBL1>^<4. 4568?3S'65*EXP<-TCl^X>-3. 45762909B’^EXP<-TC2%X>+7. 35133E-4% 
3EXP<-TC4’iX> >~jl8’^SLP^<. 9885j|.66724E~83>^< 1-EXP <-TCl>t^X> >-. 987498S784E- 
483=^< J -EXP^-TC^’fcX) >+. 6982 198058E-05*< 1-EXP <-TC4^X> > >-VCl^<3. 5814005 
56’f'EXP<-TC2^X>-. 1052867643*EXP<-TC4^X>>/3. 396il3796> 

GO TO 900 

570 MR J 2=^ < VS2B- VS2A > /X 

S\jm^<VCS^2. 3F.-8-VC4’»=6. 8F-6+VS2R^6. 8F-6>4:EXP<-TC4^X>/6. 833E-6 
«iUM2=9. 452E-:s^i IR.12^ < Jl -EXP<-TC4^X> > 

VR9‘ SUM1+5LIM2 
GO TO 900 

600 MR32“<VS1B-VS1P>/X 

VRl^ <V5.lft-VCjl > ♦EXP<-TC2*X>+MR'?2^<i-EXP<-Tr2’tX^ *'/lC2 
VR2= VPl 

IF<VR2. LT. -5. 0>VR2=-5. 0 

1F<VR2. GT. 5. 0>VR2=5. 8 

VS2fi=^VR12 

VRj 2=^VR2 

V2.2P-VR12 

VS2HB‘-1 

IF<RBS<VS2ft>. LT. 5. 8>G0 TO 630 
VS2P S I CiN < 5. 0> VS2R > 

VS2RB VS2RB-^1 

630 IF*:HBr“V:vS2P>. LT. 5. 0>CiO TO 658 
V?2B= MGNC5. 0. VF.2B:> 

V“.2HR=A'^S.:*RP-»-2 

650 GO TO <66.0. rVrO. 570.. 57C>. VS2RB 
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660 SUH1«1. 463486024E5«<3. 3E-S*i«VC8-6. 8E-6*i«VC4>*i€XP<-TC4«icX> 

SUK4»2. 9363214TOE--e4«><VSlfr-Wl>^<TC2>NEXP<-TC»iiX>-T(>l*MEXP^ 

Sm^2. 930321408E-04«i«m32«i‘<E>a><-TC4*i‘X>--EXF<-TC2*>X>> 
VR9«SUH3-i>SUH4<^SUI« 

GO TO 900 
900 qOUTlNUE 
RETURN 
EN> 

SUBROUTINE C0ND2<VR2> VR9^ VR12> 

C 

C REVISED 26 SEPT 1975 
C 

COHMON VCI>VCl^VC4.VC5.VC8,TCl^TC2.TC3,TC4>TC5,TC6,VfiLl#VflL2>X.SLP 
RE«- m±2.m32 
INTEGER VS2RB 

VR2=*5»<l--E>a=><-TC3<«X> >^VR2»iCXP<-TC3<«X> 

IF<VR2. LT. -5. 0>VR2*--5. 8 
1F<VR2. GT. 5. 0>VR2«5. 0 
150 GO TO <200. 300/ 300/ 300>/ VS2f» 

200 VR9*^<VC8»i‘3, 3E-8-VC4^6* 8E-6>»#tEXP<-TC4»i‘X>/6. 833E-6-0. 99517*<VR2*<TC4 
ln«EXP<-TC4%X>-TC3^EXP<-TC3*X>>/4. 937516846E4-f5. 010661914^<EXP<^TC3^ 
2X>-EXP<-TC4*X> > > 

GO TO 900 

300 HR12«<VS2B-VS2ft>/X 

SUM1»<VC8^3. 3E-8-VC4^6, 8E-6-^VS2ft>t:6;. 8E-«>^EXP<-TC4%X>/6. 833E-*6 
SWi2==9. 452E-3^MR124*<1-EXP<-TC4*X>> 

VR9*^SW^11+SUM2 
9Ml CWTIfsttJE 
RETIME 
END 

SUBROUTINE C0r^lD3<VR9/ VR12> 

COHHON VCI/ VCl/ VC4/ VC5/ VC8/ TCI. TC2/ TC3/ TC4/ TC5/ TC6/ VRLl/ VRL2/ X. SLP 

SUMR*<TC5^EXP<“TC5=t‘X>-TC6^EXP<-TC6*X) >/2. 195J 53128E4 

SU^®=<EXP<-TC6♦X>-FXP<-TC5♦X> >/2. 195153128E4 

VR12« < VC4+VC8 > ♦Slff1R+. 2180074 122E5^VC4^Sl^tB 

VR9^VC8*SUMR+2, 228163993E7>^<VC8^3. 3E-8-VC4’*6. 8E-6>*SUMB 

RETURN 

END 

SUBROUTINE CV0LT<VR32/ VRl. VR2/ VR9/ VR12) 

C 

C REVISED 24 SEPT 1975 
C 

COMMON VCI/ VCl/ VC4/ VC5/ VC8/ TCI/ TC2/ TC3/ TC4/ TC5/ TC6/ VRLl/ VRL2/ X 

VC^^VRL2-VR32 

VS1B“-10*VRI<2 

IF<VS1B. LT. -5. 0>VS1B=r-r.. 0 

IF<VS1B. CiT. 5. 0>VSlF-5. 8 

VCJ>VS1B-VR:1 

VC5:^VR2 

VC4-VR12-VR9 

VC8^‘‘VR9 

RETURN 

END 
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A.3 PSQ72IIIT MTA SELECTION PROGRM 
A.3.1 SiaMiy 

The sensor and electronics nodel prograas produce sets of data for 
particles on particular paths. These data sets Include all types of events 
in random order fOr a particular sensor and patii. The data selection progna 
was prepared to allow selection of all particles giving a particular response 
or combination of responses. The range of particles obtained can tien be cor- 
related with the lunar data for that response or combination of responses , 
thereby giving Important data for the formation of hypotheses regarding particle 
sources and transport tteory. 

The program selects the responses to be analyzed by reftrendng^ a code 
Inserted on an Input card. The output can be selected as either a printed 
listing or a CalComp plot. 

A. 3.2 Description 

The type and nunber of selections are read from a data card. This card 
defines the type of event to be selected, the velocity of the particles of In- 
terest, idiether or not the data are to be plotted, the data set reference number 
of the data to be analyzed, and the number of data sets to be recorded per list/ 
plot. 

If a list Is desired, the headings are written out; If a plot Is re- 
quired, a card Is read which defines the size of the axes and scales. The data 
required by the plotting routine to set up the axes and titles are then produced. 
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The data set to be analyzed Is read fron tape a record at a tine, and 
each record Is analyzed for confomnce with the characteristics selected on the 
Input card and eltiier plotted, listed, or rejected. When all records from that 
data set are analyzed, a check Is Instituted to determine If more than one event 
type Is to be plotted or listed on the one output medium or whether more analyses 
are to be performed. 

The selections available are as follows. 

The variable KIND, of dimension 8, selects the options by setting a 1 
In the respective array member corresponding to the Item number below: 

1. All niA events listed or plotted. 

2. Coincident film and collector grid events. 

3. Film (Hfily events. 

4. Multiple accumulator events. 

5. Multiple, adjacent, film events. 

6. Multiple, adjacent, grid events. 

7. Multiple, nonadjacent, film events. 

8. Multiple, nonadjacent, grid events. 

The desired sensor and the east sensor shielded film are selected by data 
set reference nunber. Particular velocities or all velocities are selected by 
the velocity parameter on the Input cards. If all velocities are required, VEL 
Is set to zero. The plotted output symbol Is related to IK, vdiich indicates the 
selection code. IK is a combination of the codes listed in KIND, i.e., 1 to 8 
for single plots or 24, say, for coincident, multiple accumulator events. 


A-79 



BSR 4234 


A. 3.3 Method of Use 

Three Input data cards are required If a plot of the data Is requested; 
If a printed list Is requested, the third card nist be ountted. 

The first card data requirements are: 

Coluiwn Requirement 

1-8 KIND; Place 1 In the positions corresponding to the desired options. 

9-15 VEL: the velocity of the desired selections. If all velocities are 
required, leave columns blank (Format F7.2). 

16 LOP: Insert a 1 If a plot is desired, o^erwlse leave blank. 

17-18 IDSR: Input data set reference number (Format 12); from JCL card. 

19-22 IK: Code Indicating type of selection for titles and plot symbols, 
e.g., 1 through 8 for single selections or 24 for coincident multiple 
accumulator, etc. (Format 14). 

23-24 NDSPP: Number of data sets to be recorded/plotted. If more than one 
data set or selection Is to be recorded on the same list or plot, 
another card identical to card 1 Is required with columns 23 and 24 
blank. 

The second card requires an alpha-numeric title In the first 28 columns. 
This title Is used In both the plotted and printed outputs. 

The third card is Identical to card two in program P5072S6F. A data 
tape Is required which carries the results from program P5072SGF. 

The program will repeat for each additional set of cards. 
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A.3.4 Flow Charts and Progran Listings 

A flow chart of the program Is shown In Figure Ar3 and a program listing 
Is provided on pages A-83 through A-84. 
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Figure A-3 P5072INT Data Selection 
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P50271NT 


F 


C 


c 


\ 


DUCNSION TITt.E<7>/ KIND<8> 

REat. M 

180 RERD<5. lie. EN0»500>K1ND> VEL« LCN»/ ZDSR> IK. NDSPP 
110 F0RMRT<811#Fr. 2. 11, 12, 14, 12> 

RERD<0, 120>TITLE 
NR1TE<6, 121>TITLE 

120 F0RMftT<7fl4> 

121 F0RMftT<2X, 7fl4> 

NR1TE<6, 270)IDSR, VEL, IK, NDSPP 

IF<LOP. EQ. 0>GO TO IS© 

DEFINE PLOT AXES AND TITLES 

READ<0, 120>QRG, QMN, QFR, QCON, NRG, NtIN, HFR, NCON, AXLEQ, AXLEM 
120 F0RMAT<2<2F4. 2, E8. 2, F10. 8>, 2F4. 2> 

NRITE<6, 180>QRQ, ©MN, QFR, QCON, NRG, lOMN, NFR, NCON, AXLEQ, AXLEH 
180 FORMAT<20X, 2<2F4. 2, E8. 2, F10. 8>, 2F4. 2> 

CALL PLOTS< IBUF, 1000, 0> 

CALL PLOTCl. 0, 1. 0, -2> 

CVC»2. 20259*WCON 
DTV=1XCVC 

CALL LQAXIS<0. 0, 0. O, 5H-LN M, -5, AXLEM, 0. O, NFR, DTV> 

CVC«2. 3025&*QCON 
DTV=1/CVC 

CALL LGAXIS<0. 0, 0. 0, 5H-..N Q, 5, AXLEQ, S0. , QFR, DTV> 

VRX-RXl ^ 

CftLL SVMPO!.<R. 5, VRX/ , 14. IGHVELOCITV <M/S> . 0. e. 15) 

IF<VEL. EO a 0>GO TO 115 

CftLL NUMRbR<9S»9. . 999. . . 14. VEL. 0. 0. -1) 

GO TO 125 

115 CftLL SVMB01 a 999. . 999. . . 14. 5H ftU. . 0. 8. 5) 

125 CftLL SYMBOL <91:^9. .999. .. 14. 20H DftTft SELECT CODE .0.0.20) 

VK=IK 

rftLL NUMBKR<999. . 999. . . 14. VK. O. O. -1) 

'/ftX--VftX-0. 25 

CftLL SYMBOL <0. 5. YftX. . 14. TITLF. fl. O. 28) 

GO TO 20C 
190 MRITE<6. 140) 

140 F0RMRT<22X. 'VELOCITY^ IIX. ‘ CHARGE i:<X. ^MftSS*'. 18X. ^GRID 7X. *^FI 
2LM ID'.BX. *" PHft'.BX. ftCC*V/> 


SELECT OHTft REOLURED BY INPUT CODES 
200 READ*: IDSR. ENr>-*-:<.50>*|:VEL. Q. M. IGTOT. IFTOT. NPHR. NRCC 
IF < VEL. EO. 0. 0>GO TO 210 
1F<1-VEL. ME. VFDGO TO 200 
21 IF<KIN['K-i \ FO. 1>G0 TO 286 
IF<KINIK2 *. r*r. 0>G0 TO .•*20 


IF<IGTO‘L FO. 0>GO TO 200 


:20 

IF 


FO. 

0>GO TO 

2j:-0 


IF 

OGTUI. GT. O: 

ro 2 t: 

JG 

!70 

IF 


. Ft* 

0.:*G0 10 

240 


IF 

<Nftcc i;i 

. 2>( 

'■*0 TO 2'0»' 

'1 

‘40 

IF 

<KTMDv5> 

Ft!. 

0>GO T(* 

250 


IF 

»‘JFTni F 

(■' :l 

‘GO TO 

»G 


IF 

a FT 01 F 

(!. *"•! 

♦GC* TO :r. 

lO 


IF 

< IFTOT. F 

C* V' ! 

A^0 10 

*0 


IF 

OFTMT. [ 

t! i; 

’ •(io 10 ; 

>■ iH 


IF 

OF TOT F 

!.♦ t ^ 

i>uo III ; 

.i 1 


IF 

*:iFTcrr. f 

0 

.•..'0(1 TO 

•'! il i 
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GO 1 fl .^‘00 

.750 IF <L.IMI l:i*» 0>GO TO 2G0 

iF'.‘. IGTOl Fr.‘. / >G0 TO 2Fr» 
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1F<I0T0T. 1:0. 6>00 TO 260 
1F<1QT0T. ECL 7>00 TO 2€0 
IF<I0T0T. EO. 12>G0 TO 2E0 
1F<IQT0T. EO. 14>Q0 TO 260 
ir<lGTOT. EO. 15>Q0 TO 260 
GO TO 200 

260 lF<Kim<7>, EO. 0>Q0 TO 270 
1F<IFT0T. EQ. 5>G0 TO 270 
1F<IFT0T. EQ 9>Q0 TO 270 
IF<IFTOT. EQ. 10>G0 TO 270 
GO TO 200 

270 IF<KIND<8>. EQ. 0>Q0 TO 280 
IF<IGTOT. EO. 5>G0 TO 280 
IF<IQTOT. EO, 9>Q0 TO 280 
IF<IGTOT. EQ. 10>G0 TO 280 
GO TO 200 
280 CONTINUE 

IF<LOP. EQ. 0>6O TO 30O 
C PREPftRR TO PLOT PHR DRTB 

phr-nphr 

QVRL«fiLOG<Q> 

QSV=^QVRL/<-QFR> 

RAND«RLOQ<QSV> 

PBNDQ^^RfiND’t'QCON 

WVfiL=^RLOG<ri> 

WSV=^WVAL/<-WFR> 

RRHD--FiLOG<WSV> 

RRNDM^^RftND^WCON 
C PLOT PHR DRTR 
INTEQ«IK 

IF<IK. GT. e>INTEQ-0 

CBLL SYMBOL <:RRNDM> RRNDQ> . 04. INTEQ. 0. O. -1> 

CRLL NUMBFR<999. . 999. . . 07. PHR. O. ©. -1> 

GO TO 200 
C LIST DRTR 

300 WRITE<6.. 160>tVEL. Q. M. IGTOT. IFTOT. NFHR. NRCC 
160 F0RMRT<21X.* 3<tl0. 4. 8X>. 14. iOX. 14. IIX. 13. 9X. 13> 

GO TO 200 
350 RFWINU IDSR 
J«J+1 

IF<J. EQ. N0SPP>GO TO 400 

RFRD<5. 110.. Fr4D==600>Kir4D. VEL. LOP. IDSR. IK 

NR I TEC 6. 370 > IDSR. VEL. IK 

GO TO 200 

370 F0RMRTC5X. • DRTR SET 13. 5X. 'VEL ^ F7. 2. 5X. ^SELECTION '.2I4> 
400 C.0NTn4UE 

IFCLOP. EO. 0>GO TO 450 

CRLL SYMBOL <20. . 8. O. . 14. . 8. 0. > 

CRLL SYMROl <20. . 10. . . 14. 11. 8. 8. -2> 

CRLL PI 0T<22. . -1. -3> 

458 GO TO :IO0 
500 CONTINUE 

IF' LC»P. EQ. 0>GO TO 510 
CRLL PLOT<0. . 8. . 999 > 

518 CONTINUE 
on TO 70 i 
608 NRITFCf.. 170> 

170 FCiRHRI <21 X. ‘'INSUFFICIENT DRTR CRRDS' ) 

700 f.CiNlINUE 
END 
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APPENDIX B 

SUMMARY STATUS MD 
PROPOSED TASKS 
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CHARACTERIZE PARTICLES RELATIVE TO OBSERVED RESPONSES 
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BENDIX PROJECT ENGINEER 
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SUPPORT LUNAR SURFACE OPERATIONS 
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